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' Recent theoretical work on solid-state proposals for the implementation of quantum computation 

^ and quantum information processing is reviewed. The differences and similarities between micro- 

^ , scopic and macroscopic qubits are highlighted and exemplified by the spin qubit proposal on one 

1 side and the superconducting qubits on the other. Before explaining the spin and supercondcut- 
^ ' ing qubits in detail, some general concepts that are relevant for both types of solid-state qubits 

Q ' are reviewed. The controlled production of entanglement in solid-state devices, the transport of 
_ carriers of entanglement, and entanglement detection will be discussed in the final part of this 
( , review. 
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I. INTRODUCTION 

The capabilities of information processing devices are derived from their physical properties; in Landauer's words, 
"Iriforinatioii is physical" (Landaucr, 1991). The pioneers of quantum information processing recognized that if a 
device was quantum mechanical, then it could have computational powers exceeding those of a classical machine. A 
sign for the superiority of quantum hardware is that typical simulations of quantum systems on classical computers 
appear to be computationally hard. 

This article is intended to give an overview of the theory of solid-state quantum information processing. For a 
general introduction to quantum computation (QC) and quantum information, we refer the reader to (Nielsen and 
Chuang, 2000). Although the distinction between different quantum devices is probably less fundamental than that 
between quantum and classical devices, Landauer's motto can also be applied here. In other words, the specific 
physical properties of the quantum hardware do matter. Two rather different categories of this hardware are those 
involving atomic systems, e.g., 

• atoms in an ion trap, 

• atoms in an optical lattice, 

• ensemble of nuclear spins in a liquid, 
and those involving solid-state systems, e.g., 

• spins of electrons in semiconductor quantum dots, 

• nuclear spins of donor atoms in a semiconductor. 
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• superconducting microcircuits containing Josephson junctions. 

This list is by no means complete; an informative collection of various proposals can be found in (Braunstein and Lo, 
2000). While there have so far been more successful demonstrations involving atomic systems in the laboratory, many 
solid-state systems are scalable, i.e., one can fabricate systems with many quantum bits (qubits) using essentially the 
same fabrication technique that is proposed or used for a single qubit. 

A. What actually has to be achieved? DiVincenzo's criteria 

For the following discussion of attempts to implement a quantum computer (or parts of it) in solid-state systems, 
it may be useful to review what actually has to be achieved. An excellent summary of the criteria for the physical 
implementation of quantum computation are DiVincenzo's following "five requirements" (DiVincenzo, 1997, 2000). 

1. A scalable physical system with well characterized qubits 

A quantum bit, or qubit, is a suitable quantum-mechanical two-state system (see item LA. 3 for more about what 
it means for the qubit to be quantum mechanical). A pure state of the two-state system then takes the form 

|^)=a|0)+/3|l), (1) 

where the amplitudes a and /? are complex numbers such that |ap + = 1. The states |0) and |1) form an 
orthonormal basis of the Hilbert space 7^2 = span{|0), |1)} of the qubit. A good example of a quantum two-state 
system is the spin 1/2 of an electron, where | t) = |0) and | [) = |1). The Hilbert space of the entire system then 
needs to be a tensor product of a large number n of such two-state systems, 

7^ = T^f " = Hz O • • • O . (2) 

n factors 

An excellent tutorial on the physical meaning of the tensor product in Eq. (2) and the difference between classical and 
quantum bits can be found in (Mermin, 2003). A system is scalable if it can be realized (in principle) for arbitrary 
n. Some early atomic qubit realizations are not (easily) scalable, and one of the biggest motivations for studying 
solid-state qubits is the hope that they will be scalable like conventional solid-state integrated circuits. A collection 
of identical particles, e.g., the Fermi sea of electrons in a metal, typically does not represent a set of well characterized 
qubits. The qubits need to be "labeled" in order to make them distinguishable, e.g., in an arrangement where single 
electrons sit on localized sites (quantum dots, donor levels of impurity atoms) and can be addressed, e.g., as "spin of 
the i-th dot" . 



2. The ability to initialize the state of the qubits 

Before a quantum computation is started, a fresh register of qubits, e.g., in the state 

|^) = |O)®" = |O)0|O)®---®|O), (3) 

^ V ' 

n factors 

is required. This requirement looks more innocent than it actually is, since it is not always easy to create such 
low-entropy states, e.g., if the temperature is not sufficiently low. Suppling a quantum computer with fresh "zeros" 
is also essential for quantum error correction, where the entropy that accumulates due to decoherence is pumped out 
of the quantum memory (Nielsen and Chuang, 2000). For this purpose, it also matters how fast the fresh "zeros" can 
be supplied. 

3. Long relevant decoherence times, much longer than the gate operation times 

A decoherence time characterizes how long it takes until the quantum phase coherence of a system (e.g., a qubit) is 
lost due to its interaction with the environment. Frequently used figures of merit are the so-called energy-relaxation 

time Ti and the decoherence time T2 of a single qubit (the notation originates from the NMR literature). To illustrate 
the meaning of these two quantities, let us assume for the moment that Ti ^ T2 (this need not be the case; Ti and T2 
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can also be of the same order of magnitude). A pure state of a qubit, Eq. (1), has degraded to an incoherent mixture 
after a time of the order of T2, described by the density matrix 

p=|ano)(0| + |/3ni)(l|. (4) 

An elementary introduction into the meaning of the density operator p can be found in quantum mechanics textbooks 
or in (Nielsen and Chuang, 2000). Note that, in the case where |0) and |1) are eigenstates with different energies, the 
processes involved in the dccohcrcnce of the qubit so far did not involve any energy exchange with the environment. 
Nevertheless, this loss of the phase information is sufHcient to disrupt a quantum computation. After a time of the order 
of Ti, energy relaxation has taken place and the system is in the thermal equilibrium state p = exp{—H/kBT), 
with the partition sum Z = Tr cxp(— iJ/fcsT), the qubit Hamiltonian the temperature T, and Boltzmann's 
constant kB- The requirement for quantum computation is that T2 ^ Top where Top denotes the time to perform a 
typical operation from the universal set (see Eq. (6) below). In order to achieve quantum computations of arbitrary 
length with the help of quantum error correction, it is required that the error probability per gate e w Top/l2 is 
below its threshold value ethres for fault-tolerant quantum computation (Nielsen and Chuang, 2000). The number 
ethres depends on the error-correcting codes used and the type of errors they have to protect against. 



4. A universal set of quantum gates 

It is required that there is a set S of unitary operators, called gates or quantum gates, operating on a bounded 
number of qubits at a time, from which all unitary operators U on any number of qubits can be composed by applying 
them in series, 



U = UkUk-i---U2Ui, 



(5) 



where Ui £ S. It has been shown that there are universal sets consisting of quantum gates that operate only on one 
or two qubits (DiVincenzo, 1995), e.g., the union of one suitable two-qubit gate U(2) with the set of all operations on 
a single qubit. 



S = {U(2)}^SU{2). 



(6) 



Examples of suitable two-quits gates U[2) are the CNOT gate, also known as quantum- XOR or simply XOR (Barenco 
et al., 1995a), with the following matrix representation in H2 <8i 'H2, 



Uxor = 
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1 



(7) 



or the square-root of SWAP (Loss and DiVincenzo, 1998), 



5 = 



/I 










1-i 

\ 








(8) 



It should be added here that there are ways of achieving unitary gates by performing non-unitary operations on a 
larger Hilbert space. There have been several proposals for doing universal quantum computation by performing von 
Neumann measurements on a subset of the qubits of a large entangled state or cluster state (Leung, 2004; Nielsen, 
2003; Raussendorf and Briegel, 2001). Another example where measurements are used to generate unitary gates is 
that of free-electron quantum computation (Beenakker et al., 2004). 



5. The ability to measure specific single qubits 

At the end of a computation, the qubits (or, at least, a subset of them) need to be measured individually in some 
fixed basis, e.g., the computational basis given by the states |0) and |1). The observable that is measured in this 
procedure is the Pauli matrix 
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B. Microscopic vs. macroscopic solid-state qubits 

The existing and proposed solid-state qubits can roughly be grouped into two categories. The qubits of the first 
category, which we will label microscopic, are similar to the atomic qubits in the sense that they are based on 
quantum objects on the atomic scale whose states |0) and |1) are distinguishable only by measuring a microscopic 
observable, such as an angular momentum on the order of Planck's constant or a magnetic dipole moment of the 
order of one Bohr magneton, ^b- Electron and nuclear spin qubits, as well as the orbital state of an electron in a 
semiconductor quantum dot, fall under this category. The second category of qubits we call macroscopic, for their 
distinguishability under measurement of a macroscopic observable, such as a current carried by a large number of 
electrons, the magnetic field induced by such a current, or the position of an electron charge in a system with two 
macroscopically distinguishable sites. The typical examples in this category are the superconducting qubits (with 
exceptions). 

C. Scope of this review article 

This is not intended to be a comprehensive review of all theoretical work that has been done in the field of solid- 
state quantum computation. Besides the discussion of some general concepts that apply for a broad range of possible 
implementations in Sec. II, we concentrate on qubits based on the electron spin (Sec. Ill) in semiconductor structures 
(quantum dots) and on superconducting circuits (Sec. IV), representing an example of a microscopic and a macroscopic 
qubit. 

Other solid-state proposals for quantum computation include quantum Hall systems (Privman et ai, 1998; Yang 
et al., 2002), anyons in fractional quantum Hall systems (Kitaev, 2003), the nuclear spin of donors in a semiconductor 
(Kane, 1998), electron charge degrees of freedom in quantum dots (Barenco et al., 1995b; Brum and Hawrylak, 1997; 
Landauer, 1996; Tanamoto, 2000; Zanardi and Rossi, 1998), "flying" electron spin qubits in surface acoustic waves 
(Barnes et al, 2000) or ballistic quantum wires (Popescu and lonicioiu, 2004), ferroelectrically coupled quantum dots 
(Levy, 2001), excitons (Biolatti et al., 2000; Chen et al, 2001; Troiani et al, 2000), SiGe quantum dots (Vrijen et al, 
2000), paramagnetic impurities in semiconductor quantum wells (Bao et al, 2003), Si-based solid-state NMR (Ladd 
et al, 2002), and electrons on the surface of liquid He (Platzman and Dykman, 1999). 

II. GENERAL CONCEPTS 

A. The Loss-DiVincenzo proposal 

The underlying idea of this proposal is that the spins of single electrons confined in semiconductor quantum dots 
(e.g., in a two-dimensional semiconductor heterostructure) are to be used as qubits (Loss and DiVincenzo, 1998). The 
required coupling between the qubits in this case is provided by the tunneling between adjacent quantum dots, giving 
rise to a nearest-neighbor exchange coupling. The resulting spin Hamiltonian is that of the Heisenberg model, 



where denotes the spin operator of the electron in the i-th quantum dot and the exchange energy between spins 
i and j. 

It has to be noted, however, that this proposal for exchange-based QC extends far beyond electron spins in quantum 
dots. Subsequent proposals for QC, using the nuclear spins of donor atoms buried in a silicon substrate (Kane, 1998), 
or using electron spins in SiGe quantum dots (Vrijen et al, 2000), electrons trapped by surface acoustic waves (Barnes 
et al, 2000), and spins of paramagnetic impurities (Bao et al, 2003), rely on the same type of interaction. 

In Eq. (10), we have also taken into account the Zeeman coupling to an external magnetic field B which may 
be spatially varying. It may also be that the Lande g-factor gi is also be site-dependent in some semiconductor 
heterostructures. The Bohr magnetic moment is denoted hy hb- Structures with two coupled quantum dots where 
the electron number could be controlled one- by-one down to a single electron per dot have recently been demonstrated 
in GaAs-AlGaAs heterostructures (Elzerman et al, 2003), see Fig. 6. 

In the "idle" phase, i.e., when no quantum gates are performed on the register, the exchange coupling would be 
switched off Jij = between all dots i and j. In order to perform an elementary two-qubit gate between dots i and j, 
the exchange coupling between dots i and j is temporarily switched on, while leaving the other exchange couplings off. 




(10) 



{id) 




back gates 



magnetized or 
high-g layer 



■ heterostructure 
quantum well 



FIG. 1 Schematics of a quantum-dot array for quantum computing according to (Loss and DiVincenzo, 1998). Quantum dots 
(dashed circles) arc defined in a two-dimensional semiconductor heterostructure with metal gates (shown schematically in grey) 
and host one (excess) electron (e) with a spin 1/2 each. By controlling the gate voltages, the coupling of adjacent quantum 
dots is switched on and off for quantum gate operations. 



(a) 



(b) 



FIG. 2 Circuit notation of two exam- 
ples of two-qubit gates that are universal 
for quantum computation when combined 
with single-qubit gates, (a) The 'square- 
root-of-swap' (S) gate, (b) the XOR gate. 



Several non-overlapping pairs of qubits can be coupled simultaneously in this way. A pulse Jy {t) with the property 

i j J,j{t')dt' = I (mod27r) (11) 

generates the above-mentioned square-root of SWAP gate (up to an unimportant global phase factor e~"^/^ which we 
omit below), 



S ~ exp j dt'H{t')^ = exp (i|Si • S^) . 



The quantum gate S can be combined with single-spin rotations 

to produce a controlled phase flip (CPF) (Loss and DiVincenzo, 1998), 

UcPF = e-'^e'^^'e-'^^'^Se''^^'S, 

which, up to a basis change, equals the quantum XOR gate: 

Uxor = VUcpfV\ 
V = exp(-i7r5|/2). 



(12) 

(13) 
(14) 



(15) 
(16) 



The effect of an inhomogcncous external magnetic field on the exchange interaction and the robustness of the 
procedure described here are discussed in (de Sousa et at, 2001). 



B. QC with anisotropic couplings 

1. Ising and transverse (XY) coupling 



Both for photon- mediated spin-spin coupling in a semiconductor microcavity (Imamoglu et al, 1999) and for 
inductively coupled superconducting qubits (Makhlin et at, 1999, 2001), which will both be further discussed in 
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Sees. III.E and IV below, the coupling takes an anisotropic form instead of being described by the isotropic Heisenberg 
Hamiltonian Eq. (10). In both cases, the form of the couphng turns out to be that of the XY (transverse) spin 
Hamiltonian, 



/O \ 

10 
10 

Vo 0/ 



(17) 



where we chose the basis of the two interacting qubits for the matrix representation of Hxy- 

It is known that any generic two-qubit Hamiltonian gives rise to a universal set of gates when combined with 
single-qubit operations. In two notable cases of anisotropic spin couplings, the Ising and the XY interactions, it is 
known how the CPF and XOR gates can be constructed. In the case of a system described by the Ising Hamiltonian 
Hi = JS'f 5*1 and a homogeneous magnetic field in z direction, there is a particularly simple realization of the CPF 
gate with constant parameters, namely Ucpf = exp(i7r(l — 25^ — 25*1 + 4S'fS'|)/4) (Loss and DiVincenzo, 1998). 

For the transverse spin-spin coupling of Eq. (17), we find that a useful two-qubit gate, such as the conditional- 
phase-flip (CPF) operation, can be carried out by combining Hxy with one-bit rotations. The unitary evolution 
operator generated by the Hamiltonian of Eq. (17) is 



Uxy{cI>) = Texp 



dtH 



XY 



exp[icf>{StS^ + SfS]) 



(18) 



where <p = J dtJ{t). The CPF gate (C/cpf) can be realized by the sequence of operators (Burkard et al, 1999b; 
Imamoglu et al, 1999) 



t/cPF = e^^/'^e*'^"' /^e*''"^ '^^ I^Uxy (7r/2)e"''^ Uxy {-k /2)e''"'l'^e'-'"'y''^ 



(19) 



where a denotes the vector Pauli operator, where S = cr/2, and = (1,1,— 1)/\/3 and Wj = (— 1, 1, l)/\/3. The 
XOR gate can be realized by combining the CPF operation with single-qubit rotations as in Eqs. (15) and (16). 

While it is impossible to generate the CNOT gate with a single use of the XY Hamiltonian (Burkard et al., 1999b), 
it is possible to generate a different universal quantum gate with the XY interaction in a single pulse; the CNOT + 
SWAP (CNS) gate [/cns = t^swApC^xoR, is generated as (Schuch and Siewert, 2003) 



J/cNS = HiUxY{'^)e 



2, 



where Hi is the Hadamard gate 



H=Xf' \ 



(20) 



(21) 



applied to qubit i. 

Gate errors due to unwanted inhomogeneous magnetic fields during an otherwise isotropic coupling, effectively 
creating an anisotropy, have been studied and quantified in (Hu and Das Sarma, 2003). 



2. Anisotropy due to the spin-orbit coupling 

The exchange interaction, Eq. (10), between electron spins in tunnel- coupled sites (such as quantum dots) can 

acquire anisotropic terms due to spin-orbit coupling during tunneling between the sites (Kavokin, 2001). Surprisingly, 
it turns out that the first-order effect of the spin-orbit coupling during quantum gate operations can be eliminated 



G 



G+ 



FIG. 3 A circuit representation for 
the conditional phase flip (CPF), 
Eq. (14). The single qubit rotations 



arc G 



and 



= e"''^ . The CPF is related to 
the XOR gate Eq. (7) by the basis 
transformation Eq. (15). 
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by using tinic-symmctric pulse shapes for the couphng between the spins (Bonesteel et al, 2001). A related, but 
independent, result shows that the spin-orbit effects exactly cancel in the gate sequence on the right hand side 
of Eq. (14) required to produce the quantum XOR gate, provided that the pulse form for the spin-orbit and the 
exchange couplings are identical (Burkard and Loss, 2002). The XOR gate being universal when complemented with 
single-qubit operations, this result implies that the spin-orbit coupling can be dealt with in any quantum computation. 
In any real implementation, there will be some (small) discrepancy between the pulse shapes for the exchange and 
the spin-orbit coupling; however, one can choose two pulse shapes which are very similar. It was shown that the 
cancellation still holds to a very good approximation in such a case, i.e. the effect of the spin-orbit coupling will still 
be strongly suppressed (Burkard and Loss, 2002). There is also an effect of dipolar interactions between adjacent 
spins, providing another anisotropic coupling; this coupling can also be treated as an anisotropic contribution to 
Eq. (10) and therefore cancels out in the gate sequence Eq. (14) for the same reasons as the spin-orbit interaction. 

The spin-orbit coupling for a conduction-band electron is given by the following Hamiltonian (Gantmakher and 
Levinson, 1987), being linear in the 2D momentum operator pj, i = x,y ([100] orientation of the 2D plane), 

H,o= Pij'^iPh (22) 

where the constants depend on the strength of the confinement in z-direction and are of the order (1-^3)- 10^ cm/s 
for GaAs hetcrostructures. Combining the isotropic Heisenberg coupling (10) with the anisotropic exchange between 
two localized spins Si and S2 one obtains the Hamiltonian (Burkard and Loss, 2002) 

H{t) = J{t){Si-S2+A{t)), (23) 

where the anisotropic part is given by the expression (Kavokin, 2001), 

A{t) = m ■ (Si X S2) + imm ■ si)im ■ sa), (24) 

and Pi = Y2j Ai(V-'i|*?'jlV'2) is the spin-orbit field, the groimd state in site (dot) i = 1, 2, and 7 » O{0^). As was 
discussed in Sec. 11. A, for ^ = 0, the quantum XOR gate can be obtained by applying H{t) twice, together with 
single-spin rotations, see Eqs. (14) and (16). Moreover, ii A = Q, then H{t) commutes with itself at different times 
and the time-ordered exponential 



U{ip) = Texp / H{t) dt (25) 




is a function of the integrated interaction strength only, = /I^^/g J{t)dt. In particular, U{ip = 7r/2) = Uli^ = S \s 
the "square-root of swap" gate. 

The interesting situation is of course 0. If in this case, /3 and 7 (and thus A) are time-independent, then H{t) 
still commutes with itself at different times and one can find a fixed coordinate system in which /3 is parallel to the 
z axis. In this basis, the anisotropic term Eq. (24) can be expressed as 

A = 0{SfS^-SfS^) + 6StSI, (26) 

with S = 7/?2. In the singlet-triplet basis with basis vectors {|T+) = |tT), \S} = (| Ti) " UT))/\/2, \To} = (|ti) + 
1 11))/^/^, \T-) = I II)} the gate sequence Eq. (14), including the anisotropy Eq. (24), yields the unitary operation 

Ug = diag(ze-^^(i+^\ 1, 1, -ie-*^(i+*)), (27) 

where diag(xi, . . . ,2:4) denotes the diagonal matrix with diagonal entries Xi, . . . ,X4. Note that the pulse strength 
if and the spin-orbit parameters only enter Ug in the = ±1 subspaces. Moreover, the terms linear in f3 have 
canceled out exactly in Ug. With the choice (p — tt/2{1 + 6), one obtains the conditional phase flip gate Ug = Ucpf = 
diag(l, 1, 1, — 1), being equivalent to the XOR up to the basis change, Eq. (16). Therefore, the anisoptropic terms 
A = const, in the spin Hamiltonian cancel exactly in the gate sequence Eq. (15) for the quantum XOR. 

We briefly discuss what happens if, as can be expected in real systems, the anisotropic terms in the Hamiltonian H 
are not exactly proportional to J{t), i.e. if A(t) is time-dependent. Generally, both (3 and 7 depend on time. In this 
more general case, we cannot exactly eliminate the effect of the anisotropy because of the time-ordering in Eq. (25) 
and since the Hamiltonian cannot be expected to commute with itself at different times, [H{t)^H{t')] 7^ 0. The 
estimated gate errors e = ||J7g — Ucpf\\^ due to the anisotropy in the case where A{t) is only weakly time-dependent 
are e < where we use A/3(f) = /3(t) — (3q and 

A = (|^|/3o/2) max |(J(t)/Jo)(/3(0//3o - 1)|, (28) 
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where Jo stands for the average exchange couphng, Jq = f/rg ^ 0. It can be shown (Burkard and Loss, 2002) that for 
tunnel-coupled quantum dots, it is possible to choose a weakly time-dependent A by using Eq. (52) for the exchange 
coupling and the result 

b{d, q) = \J{d, q)f3{d, q)\ = ho^de^{-2qd^), (29) 

where 6o = a/a%, a% = ^Jh/muJQ, and where a is a constant depending on the spin-orbit parameter (for a 5 nm 
wide [100] GaAs quantum well a w 2 mcVnm), q = uj/luq. The minimal value of the quantum dot confinement 
energy cu is denoted by wq- A possible model for the switching process is the use of a time-dependent confinement 
strength q{t) = ij{t)/u)o = cosh^{at/Tg) (where alpha is a number of order 1, e.g., a = 4). This pulse shape has 
favorable adiabatic properties (Burkard et al., 1999a; Schliemann et al., 2001b), as detailed in Sec. II. E. and leads to 
a pulsed exchange interaction J(t) — J{d,q{t)) and spin-orbit field b{t) = b{d,q{t)), where —Ta/2 < t < Ts/2. The 
resulting error was estimates in (Burkard and Loss, 2002) as A « 7 • 10~^, leading to gate errors occurring at a rate 
e w 4A^ « 2 • 10~^ being around the currently known threshold for fault tolerant quantum computation (Preskill, 
1998a,b). The error e can be further reduced by performing the gates more slowly, with a long period of constant A 
between the rise and fall of the pulses. 

C. Universal QC with the exchange coupling 

In some situations, a local controllable field or g-factor gi in the Hamiltonian Eq. (10) and thus the single-qubit 
operations SU{2) in the universal set Eq. (6) may be more costly to implement than the tunable exchange coupling 
generating the spin-spin coupling U(2) (note, however, that there exist all-electric switching schemes using g-factor 
modulation, sec Sec. III.D). A scheme has been developed in which the Heiscnberg interaction alone suffices to exactly 
implement any quantum computer circuit, at a price of a factor of three in additional qubits and about a factor of 
ten in additional two-qubit operations. However, the Heisenberg interaction by itself is not a universal gate (Barenco 
et al., 1995a), in the sense that it cannot generate any arbitrary unitary transformation on a collection of spin- 1/2 
qubits. This is why in Eq. (6), the Heisenberg interaction needs to be combined with some other means of applying 
independent one-qubit gates. The Heisenberg interaction alone docs not give a universal quantum gate because it has 
too much symmetry: it commutes with the operators S"^ and S^, where the total spin is defined as 

n 

S = ^S,, (30) 
and therefore it can only rotate among states with the same S, Sz quantum numbers. 

1. Encoding 

The exchange coupling is thus not universal in the full Hilbert space; but, by working exclusively in one symmetry 
sector of the Hilbert space with fixed S, Sz quantum numbers, the exchange coupling can be made universal. This 
restriction is achieved by defining coded qubit states, ones for which the spin quantum numbers always remain the 
same (Bacon et al., 2000; Kempe et al., 2001; Viola et al., 2000). The smallest number of spins 1/2 for which two 
orthogonal states with identical S, Sz exist is three. The space of three-spin states with spin quantum numbers 
S = 1/2, Sz = +1/2 is two-dimensional and will serve to represent our coded qubit. An explicit choice for the basis 
states of this qubit are 

|0l) = \S)\ T), (31) 
lU) = V273|T+)U)-\/l73|To)|T), (32) 

where 15) = a/1/2(| t^) — | J,t)) is the singlet state of spins 1 and 2 (see Fig. 4a) of the three-spin block, and 
|T+) = I IT) and |To) = ^/l/2{\ tJ,) -I- | 11)) are triplet states of these two spins. While in principle this solves the 
problem of exchange-only QC, in practice we would like to know what the overhead in terms of qubits (for coding) and 
gates (for operating on encoded qubits with the exchange interaction) will be, and how a universal set of operations 
on the encoded qubits can be achieved (DiVincenzo et al, 2000). It has also been found that the anisotropic XY 
interaction (17) alone is sufficient for quantum computation (Kempe and Whaley, 2002), a result which was later 
generalized to large class of anisotropic exchange Hamiltonians (Vala and Whaley, 2002). An encoding involving two 
spins per qubit has also been demonstrated for universal quantum logic starting from locally alternating g-factors 
(Levy, 2002) and from a homogeneous magnetic field combined with anisotropic exchange interactions (Wu and Lidar, 
2002a,b; Wu et al., 2004). 
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FIG. 4 Possible layouts of spin-1/2 de- 
vices, a) One-dimensional layout. We con- 
sider two different assumptions about how 
the exchange interactions can be turned 
on and off in this layout: 1) At any given 
time each spin can be exchange-coupled to 
at most one other spin (we refer to this 
as "serial operation" in the text), 2) All 
exchange interactions can be turned on 
simultaneously between any neighboring 
pair of spins in the line shown ( "ID parallel 
operation"), b) Possible two-dimensional 
layout with interactions in a rectangular 
array. We imagine that any exchange in- 
teraction can be turned on between neigh- 
boring spins in this array ( "2D parallel op- 
eration"). Of course other arrangements 
are possible, but these should be represen- 
tative of the constraints that will be faced 
in actual device layouts. 

2. One-qubit gates 

A one-qubit gate on a single three-spin block is performed as follows. The Hamiltonian H12 generates a rotation 
U12 — exp{i/h J J Si ■ S2dt) which is just a z-axis rotation (in Bloch-sphere notation) on the coded qubit, while 
^23 produces a rotation about an axis in the x-z plane, at an angle of 120° from the z-axis. Since simultaneous 
application of H12 and i?23 can generate a rotation around the x-axis, three steps of ID parallel operation (defined in 
Fig. 4a) permit any one-qubit rotation, using the classic Euler-angle construction. In serial operation, it can be found 
numerically that four steps are always adequate when only nearest-neighbor interactions are possible (e.g. the sequence 
H12-H23-H12-H23 shown in Fig. 5a, with suitable interaction strengths), while three steps suffice if interactions can 
be turned on between any pair of spins (e.g. H12-H23-H13, see Fig. 5b). 

3. Two-qubit gates 

The implementation of two-qubit gates for universal QC with the exchange interaction on two three-spin code 
blocks is less intuitive that the corresponding task for one-qubit gates. Much of the difficulty of these searches arises 
from the fact that while the four basis states |0l,1l)|0l,1l) have total spin quantum numbers S = 1, Sz = +1, the 
complete space with these quantum numbers for six spins has nine states, and exchanges involving these spins perform 
rotations in this full nine-dimensional space. Numerical searches for the implementation of two-qubit gates using a 
simple minimization algorithm (DiVincenzo et al., 2000) aided by the two-qubit gates invariants (Makhlin, 2002) have 
resulted in a sequence for an encoded CNOT operation that is depicted in Fig. 5. The solution shown in Fig. 5c 
appears to be optimal for serial operation and happens to involve only nearest neighbors in the ID arrangement of 
Fig. 4a. There also are (apparently) optimal numerical solutions for parallel operation mode. For the ID layout of 
Fig. 4a, the simplest solution found involves 8 clock cycles with just 8*4 different interaction-time parameters {H12 
can always be zero in this implementation). For the 2D parallel mode of Fig. 4b, a solution was found using just 7 
clock cycles (7*7 interaction times). 

In the present scheme, quantum computation would proceed as follows. In the beginning, all the computational 
qubits would be set to the |0l) state which is easily obtained using the exchange interaction: if a strong H12 is 
turned on in each coded block and the temperature made lower than the strength J of the interaction, these two 
spins will equilibrate to their ground state, which is the singlet state. The third spin in the block should be in the 
I 1) state, which can be achieved by also placing the entire system in a moderately strong magnetic field B, such that 
ksT << gfisB < J (it can be shown that in a slightly more general scheme involving both the Sz = +1 and Sz = — i 
subspaces, the last step can be omitted). After the computation, with the one- and two-qubit gates implemented 
according to the schemes mentioned above, the final qubit measurement, we note that determining whether the spins 
1 and 2 of the block are in a singlet or a triplet suffices to perfectly distinguish (DiVincenzo et ai, 1999) |0l) from 
|1l) (again, the state of the third spin does not enter). 
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FIG. 5 Circuits for implementing single-qubit and two-qubit 
rotations using serial operations, a) Single-qubit rotations by 
nearest-neighbor interactions. Four exchanges (double-headed 
arrows) with variable time parameters r; are always enough to 
perform any such rotation, one of the two possible layouts is 
shown, b) Non-nearest neighbor interactions. Only three inter- 
actions are needed, one of the possible layouts is shown, c) Cir- 
cuit of 19 interactions that produce a cNOT between two coded 
qubits (up to one-qubit gates before and after). The durations 
of each interaction are given in units such that for t = 1/2 the 
rotation Uij — exp{iJtSi ■ Sj/h) is a SWAP, interchanging the 
quantum states of the two spins i, j. The U parameters are 
not independent, they are related to the tiS as indicated. The 
uncertainty of the final digits of these times are indicated in 
parentheses. With these uncertainties, the absolute inaccuracy 
of the matrix elements of the two-qubit gate rotations achieved 
is no greater than 6 x 10^^. Further fine tuning of these time 
parameters would give the cNOT to any desired accuracy. In 
a practical implementation, the exchange couplings J{t) would 
be turned on and off smoothly; then the time values given here 
provide a specification for the integrated value f J{t)dt. The 
functional form of J{t) is irrelevant, but its integral must be 
controlled to the precision indicated. 



4. Protection against errors 

Codes of the type of Eqs. (31) and (32) have first been introduced as a computational basis in decoherence-free 
subspaces, i.e., subspaces of a Hilbert space which are protected against errors with a certain type of symmetry (Bacon 
et ai, 2000; Kempe et al., 2001; Lidar et ai, 1999, 1998; Zanardi and Rasetti, 1997). Moreover, it has been suggested 
that the logical subspace may be energetically separated from the remaining Hilbert space and thus protected against 
errors in a system where a certain combination of exchange couphngs always remains switched on (Bacon et al, 2001; 
Weinstein and Hellberg, 2004). 

5. Related ideas 

Encoded qubits of a different kind, so-called spin-cluster qubits (Meier et ai, 2003a, b), have been proposed in 
order to relax the requirements for control on the single-spin level while inheriting the favorable single-spin properties 
such as long decoherence time and fast gate operating time. Spin cluster qubits arc finite spin chains with Heisen- 
berg or anisotropic (XY and Ising-like) antiferromagnetic exchange interaction that can have uniform or nonuniform 
interaction constants. 

The use of many-electron QDs for exchange-based quantum computations has been analyzed in (Vorojtsov et al., 
2004). A particular implementation of three-spin QDs encoding one qubit has been put forward in (Kyriakidis and 
Penney, 2004). 

D. Optimization of quantum circuits 

A quantum gate operating on n qubits can be represented as a 2" x 2" unitary matrix. Any quantum computation 
or algorithm can be split up into a series of elementary gate operations drawn from a universal set involving only one or 
two qubits, as in Eq. (6). This is the quantum circuit representation of quantum algorithms (or, unitary operations). 
For a simple example, see Fig. 3 for a circuit representations of CPF in terms of sqrt-of-SWAP gates. However, 
quantum circuits arc in general not the most efficient way of implementing a quantum computation (Burkard et al, 
1999b). There are a number of related but different approaches using, e.g., genetic algorithms and chirped Gaussian 
pulses (Sanders et ai, 1999) or control theory (Khaneja et ai, 2001). 

If one is interested in optimizing the switching time for a desired unitary [/, with a given Hamiltonian, e.g., the 
spin Hamiltonian Eq. (10), one can depart from the circuit representation of the unitary by allowing arbitrary time 
dependent parameters p{t) in the Hamiltonian. In the case of the spin Hamiltonian (10), we have p = (J, Bi, B2, . . .). 
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Wc will only demonstrate this optimization in the case of a simple two-qubit unitary, the XOR (CNOT) gate. The 
optimization method can in principle be applied to unitaries of any size; note, however, that the optimization as an 
arbitrary classical computational task is typically a hard computation in itself. 



1. Serial pulse mode 

We first restrict ourselves to a special class of parameter functions p{t), in which at every time t, only one component 

of p{t) is non-zero. If wc further restrict ourselves to parameter functions in which the duration of the ./-pulses with 
p — (J, 0, 0, . . .) are 7r/2 pulses generating the sqrt-of-swap S, then we are back to the circuit model with the universal 
set Eq. (6) and U(^2) — S. In this case we can optimize circuits, e.g., to have as few instances of S as possible. 
E.g., it turns out that the use of two S for a CPF as in the sequence Eq. (14), and therefore, for XOR, is optimal. 
Such minimal requirements for the implementation of a unitary U can be proven by analyzing the set ViU) of 

product states ||^) e W = \^) = |^i) • • • (g) |</>m); |</'i) S which are mapped back onto product states by 

U (Burkard at al, 1999b). An alternative method for determining whether a Hamiltonian generates a gate in a single 
pulse involves the invariants under addition of single-qubit gates (Makhlin, 2002). 

2. Parallel pulse mode 

In the case where several or all parameters p can be changed simultaneously, we expect that a given quantum gate, 
say XOR, can be performed faster than by changing only one parameter at a time as in the serial pulse mode. The 
unitary time evolution operator after time t is the following functional in p, 



Ut [p(r)] = T exp ( 1 ^ iJ (p(r)) dr ] , (33) 



where T denotes the time-ordering. For a given quantum gate Ug, the integral equation Ut\p{T)] = Ug has to be solved 

for the functions p{t). An optim,al solution is given by a set of bounded hmctions |pi(T)| < Mi requiring minimal time 
t for a fixed set of bounds Mj. In order to simplify the problem, one can restrict the problem to piecewise-constant 
functions, 

f/^(pW,...,pW;</.) = e'^[7Ar(pW)...C/2(p(2))c/^(p'(i)), 

?7fc(p('=)) = exp{ttHip^''y)}. (34) 

For each of the N time intervals, one has the freedom to choose a new set of parameters p'-'^^ = (J, Bi,B2). The 
discretizcd problem can now be treated both analytically and numerically (Burkard et al., 1999b). 

One finds analytically that CPF can be implemented in a single step by fixing N = 1, i.e., all parameters in Eq. (10) 
simultaneously non-zero but constant, 

UcpF = exp [itH{J, Bi, B2)] , (35) 

The parameters are (in units of 2'jTh/t), 

J = k — n — 2m—- , 4'=~^{n + -^) 
Bi = ^(0,0,n+^ + Vfc2- J2) , B2 = ^(0,0,n+^- J2), (36) 

where n, and m are arbitrary integers, and k is an integer satisfying 2\k\ > \n + 2m + ■^l- In the specific case where 
all constraints are equal to M, we find that the solution forfc = l,m = n = 0, 

J=l BJ = i(l + V3), i3f = i(l-V3), <^ = -| (37) 
has the shortest switching time, 



(38) 
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less than half the time which is used for the serial pulse quantum circuit Eq. (14), tcPF,s = 1-5 • 2nH/M. Numerically, 
one finds that XOR requires at least N = 2 steps, 



f/xOR 



e"fe"f exp 



exp 



with the parameter values (in units of 2-Kh/t) 



k 


jik) 




^2x 








^2z 


1 


0.187 


-0.025 


0.464 


0.205 


0.195 


-0.420 


0.395 


2 


0.617 


-0.220 


0.345 


-0.384 


0.244 


0.353 


0.108 



(39) 



(40) 



and the global phase cj) = —0.8481 • tt, where the time t has to be chosen such that none of the parameters exceeds the 
bound M. The total switching time for equal boimds is in this case txOR.p = (0.4643+0. 6170)27r7i/M = 1.0813-27r?i/-^, 
compared to ^xOR.s = 2 • 2TTh/M for the serial switching. 

3. Anisotropic systems 

Parallel switching is also possible with the XY dynamics Eq. (17). It can be shown that Ucpf requires two pulses, 

UcpF = e'''U2Ui, (41) 



where Uk = exp f27riiJxy,B ( J^*"' , B^J'^ , B^*^) ) 



1,2. 



Note that all magnetic fields can be chosen homogeneous (B-j^' 



(fc) 



B 



(k) 



{By = 0). Here we give one possible realization which is found numerically (^ = — 37r/4): 



B^*^)) and perpendicular to the ?/-axis 



k 






Bi'^ 


1 


0.7500 


0.7906 


0.5728 


2 


0.5000 


0.0000 


0.2500 



(42) 



The total switching time for CPF, assuming equal bounds Mj = Mb = M for J and B, is ^ = 1.291 • 2nh/M, 
compared to ^cpf s = 2.167 • 2nh/M for the serial pulse sequence defined in Eq. (19). 

In order to produce the XOR gate Eq. (7) we can implement the basis change Eq. (15) using the single-qubit 
rotation V. This procedure requires a total of four steps for the XOR gate. Another way of achieving XOR is the 
following sequence which we found numerically and which takes only three steps: 



Uxor = exp(3i7r/4)C/3J72f/i, 



(43) 



with the following parameters: 



k 


jik) 


Bix 


B2x 




B2y 


Biz 


B2z 


1 


1.802 


0.615 


2.045 


0.020 


0.316 


0.794 


0.130 


2 


3.344 


0.348 


0.718 


0.259 


0.493 


1.583 


1.062 


3 


1.903 


1.193 


0.705 


0.413 


-0.305 


0.589 


0.604 



(44) 



The total switching time of txoR.p = 7.29- 27r?i/M (compared to 2.67 ■2'jTh/M using CPF and a basis change) indicates 
that Eq. (44) is not an optimal solution. 



E. Adiabaticity 

Quantum gates are generated by controlling the parameters in the Hamiltonian Eq. (10), Jij{t) and Bj(t) (or 
gi{t)), as a function of time. E.g., the exchange coupling J depends on time via some physically controlled quantity, 

such as an electric gate voltage v{t), i.e., J{t) = J{v{t)) (similarly for the effective g-factor g{t)). According to 
Eq. (11), only the time integral J{v(t))dt needs to assume a certain value (modulo 2it) in order to generate the 
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correct quantum gate and the pulse form of v{t) does not matter. However, the exchange interaction J{t) needs to be 
switched adiabatically in order to avoid unwanted excitations in the system. The adiabaticity condition is (Burkard 
et ai, 2000a, 1999a, b) \v/v\ ^ Se/h, where Se is the energy scale on which excitations may occur. Here, 6e denotes 
the energy-level separation of a single dot, i.e., the smaller of either the single-electron level spacing or the on-site 
Coulomb energy U required to add a second electron to a dot. A rectangular pulse leads to excitation of higher 
levels, whereas an adiabatic pulse with amplitude vq is e.g. given by v{t) — vq sech(t/Ai) where At controls the width 
of the pulse. We need to use a switching time Tg > At, such that v{t = Ts/2)/vo becomes vanishingly small. We 
then have \v/v\ = \ta,nh{t / At)\ / At < 1/At, so we need 1/At <C 6e/h for adiabatic switching. The Fourier transform 
v{lij) = Atoo7rsech(7ra;At) has the same shape as v{t) but width 2/7rAt. In particular, v{uj) decays exponentially in 
the frequency uj, whereas it decays only with l/oj for a rectangular pulse. 

Adiabatic switching of the exchange coupling in two coupled quantum dots and the error probability for different 
pulse forms have been studied numerically in (Schliemann et al, 2001b). Furthermore, corrections to the fully 
adiabatic result have been investigated (Requist et ai, 2004). 

III. ELECTRON SPINS 

Being a natural two-level system, the spin 1/2 of the electron represents an ideal candidate for a qubit. On the 
one hand, the electron spin is typically quite well isolated from charge degrees of freedom (not completely, though, 
due to, e.g., the spin-orbit coupling). In some situations, electron spin dccoherence times in solids appear to be 
relatively long, exceeding microseconds (Awschalom and Kikkawa, 1999; Kikkawa and Awschalom, 1998; Kikkawa 
et ai, 1997). On the other hand, single spins in solid-state structures are not readily available and controllable. 
However, large experimental efforts are currently made to isolate and control single spins in solid-state structures. 
The spin-based proposals for quantum information processing which will be discussed below are all based on artificial 
nano- or micrometer-scale semiconductor structures, such as quantum dots (QDs) or microcavities. 

A. Quantum Dots 

In (Loss and DiVinccnzo, 1998), a quantum register is proposed in which single electrons are trapped in quantum 
dots (QDs) that are arranged in an array or lattice in a semiconductor structure, e.g., as in Fig. 1. Electrically defined 
QDs in two-dimensional semiconductor heterostructures (typically, GaAs) are well-studied objects (Kouwenhoven 
et al., 2001) in which charge transport has attracted much attention (Averin and Nazarov, 1992; Kouwenhoven et ai, 
1997a; van der Wiel et ai, 2003). In recent years, the controlled storage of a single electron — and thus a spin 1/2 or 
qubit — in a QD has been achieved (Ciorga et ai, 2000; Elzerman et ai, 2003). Structures in which two QDs, each 
containing a well-controlled number of electrons (down to a single electron), are adjacent and tunnel-coupled, have 
been fabricated and studied (Elzerman et ai, 2003). In Fig. 6, we show an electron micrograph of a structure of the 
type that was used in (Elzerman et ai, 2003). The tunneling of electrons between the two dots is predicted to give 
rise to the spin exchange coupling JSi • S2 in Eq. (10). In the next section, we are going to outline a theory of this 
spin exchange mechanism. 




FIG. 6 Electron micrograph of a structure comprising two 
QDs, defined by metal electrodes (bright structures) on the 
surface of a GaAs/AlGaAs heterostructure (Courtesy of J. 
Elzerman, TU Delft). The charge on the dots is controlled 
in steps of single electron charges, down to one electron per 
dot, by tuning the voltage applied to the plunger gates Pl,r 
and is monitored by measuring the conductance of (i.e., the 
currents Jqpc through) the quantum point contacts (QPCs) 
Q-R and Q-L. Conductance spectroscopy was performed by 
measuring the current /dot (Elzerman et al, 2003). 
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FIG. 7 Two coupled QDs with one var- 
lence electron per dot. Each electron is 
confined to the xy plane. The spins of 
the electrons in dots 1 and 2 arc de- 
noted by Si and S2. The magnetic 
field B is perpendicular to the plane, 
i.e. along the z axis, and the elec- 
tric field E is in-plane and along the x 
axis. The quartic potential is given in 
Eq. (48) and is used to model the cou- 
pling of two harmonic wells centered at 
(±a, 0, 0). The exchange coupling J be- 
tween the spins is a function of B, E, 
and the inter-dot distance 2a. 



B. Exchange in laterally coupled QDs 

Due to the Coulomb interaction and the Pauli exclusion principle, the ground state of two coupled electron sites 
(atoms, QDs) in the absence of a magnetic field is a spin singlet (a highly entangled spin state), while the spin triplet 
states (one of them entangled) arc typically separated by some energy gap J. This energy gap is called exchange 
coupling, as it arises from virtual electron exchange between the two sites due to the interaction. The virtual electron 
exchanges are allowed for opposite spins (spin singlet, 5 = 0) but forbidden by the Pauli principle for parallel spins 
(spin triplet, S = 1), therefore the energy of the singlet is lowered by the interaction. 

We now introduce a model for the two laterally coupled QDs containing one (conduction band) electron each 
(Burkard et ai, 1999a). The two-dot system is shown schematically in Fig. 7. It is essential that the electrons are 
allowed to tunnel between the dots, and that the total wave function of the coupled system must be antisymmetric. It 
is this fact which introduces correlations between the spins via the charge (orbital) degrees of freedom. The electronic 
Hamiltonian in the effective-mass approximation for the coupled system is then given by 

H=J2 Kr^, P^) + C + Hz= iJorb + Hz, (45) 

1=1,2 

where the single-particle Hamiltonian, 

h{ri. Pi) = ^{pi- iM^i)) ^ + eXiE + V{ri), (46) 
describes the electron dynamics confined to the xy-p\a,ne and 

C = r, (47) 

K\ri — r2| 

represents the Coulomb interaction (unscreened in this case where the dot diameter is small or comparable to the 
screening length). The electrons have an effective mass m {m = 0.067 rrie in GaAs) and carry a spin-1/2 Sj. The 
dielectric constant in GaAs is k = 13.1. We allow for a magnetic field B = {0,0, B) applied along the 2;-axis and 
which couples to the electron charge via the vector potential A(r) = -|-(— y, a;, 0). We also allow for an electric field 
E applied in-plane along the x-direction, i.e. along the line connecting the centers of the dots. The coupling of the 
dots (which includes tunneling) can be modeled by a quartic potential, 

V{x,y) = ^(^-^,{x^-a^r + y'), (48) 

which separates (for x around ±o) into two harmonic wells of frequency ojq, one for each dot, in the limit of large 
inter-dot distance, i.e. for 2a ^ 2aB, where a is half the distance between the centers of the dots, and ob = y^h/mcoo 
is the effective Bohr radius of a single isolated harmonic well. This choice for the potential is motivated by the 
experimental fact (Kouwenhoven et al., 1997b; Tarucha et al, 1996) that the spectrum of single dots in GaAs is well 
described by a parabolic confinement potential, e.g. with Hluq = 3meV (Kouwenhoven et al, 1997b; Tarucha et al., 
1996). We note that in this simplified model, increasing (decreasing) the inter-dot distance is physically equivalent to 
raising (lowering) the inter-dot barrier, which can be achieved experimentally by e.g. applying a gate voltage between 
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the dots (Livcrniore et ai, 1996; Waugh et ai, 1995). Thus, the effect of such gate voltages is described in this model 
simply by a change of the inter-dot distance 2a. 

The magnetic field B also couples to the electron spins via the Zeeman term Hz = 5Mb ' where g is 

the effective g- factor {g « —0.44 for GaAs), and /zb the Bohr magneton. The ratio between the Zeeman splitting 
and the relevant orbital energies is small for all -B-values of interest here; indeed, guBB/hwo ^ 0.03, for _B <c -Bo = 
{fvbJo/ UBiim/nie) ~ 3.5 T, and gi^BB/fuvi, < 0.03, for B ^ Bq, where = eB/2mc is the Larmor frequency, and 
where we used Hloq = 3 meV. Thus, we can safely ignore the Zeeman splitting when we discuss the orbital degrees of 
freedom and include it later into the effective spin Hamiltonian. 



1. The Heitler-London approach 

We consider first the Heitler-London (HL) approximation (also known as valence orbit approximation), and then 
refine this approach by including hybridization as well as double occupancy in a Hund-MuUiken approach, which will 
finally lead us to an extension of the Hubbard description. We will see, however, that the qualitative features of J as a 
function of the control parameters are already captured by the simplest HL approximation for the artificial hydrogen 
molecule described by Eq. 45. 

The HL approximation is borrowed from molecular physics. In the present case, think of a hydrogen molecule H2. 
The HL approach starts from single-dot ground-state (s wave) orbital wavefunctions ip{r) and combines them into 
the (anti-) symmetric two-particle orbital state vector 

I*.) = ^iia, (49) 

the positive (negative) sign corresponding to the spin singlet (triplet) state, and S = J (i^r(^^^(r)(^_a(r) = (2|1) 
denoting the overlap of the right and left orbitals. A non- vanishing overlap implies that the electrons tunnel between 
the dots (see also Sec. III.B.3). Here, </?-a(i") = (r|l) and </5+a(i") = (r|2) denote the one-particle orbitals centered 
at r = (=pa, 0), and \ij) — \i)\j) are two-particle product states. The exchange energy is then obtained through 
J = et — Cg = (\l/_|iJorb|*-) — (*-i-|-fforb|*-i-)- The single-dot orbitals for harmonic confinement in two dimensions in 
a perpendicular magnetic field are the Fock-Darwin states (Darwin, 1930; Fock, 1928), which are the usual harmonic 
oscillator states, magnetically compressed by a factor b = uj/luq = \/l + luI/luq, where Wl = eB/2mc denotes the 
Larmor frequency. The ground state (energy hu = bhua) centered at the origin is 



^{x,y) = ^^e—i^'^yy^\ (50) 

Shifting the single particle orbitals to (±a,0) in the presence of a magnetic field we obtain (p±a{x,y) = 
eyip{±iya/2lg)ip{x =F a,y), where the phase factor involving the magnetic length Ib = \fhcfeB is due to the gauge 
transformation A±a = B{—y, a,0)/2 ^ A = B{—y, x, 0)/2. We obtain (Burkard et ai, 1999a) 

.= _^(a2|C + *|12)-S?ML±iffil>), ,51, 
where the overlap becomes S = exp(— mwa^/ft — a?h/Al^muj). Evaluation of the matrix elements of C and W yields 



J = 



sinh(2d2(25- i)) 



zVb(e-»''\{bd^) - e'^^(^-V'')lo(rf2{6 - 1})) + ^ (l + hd^) 



(52) 



where we introduce the dimensionless distance d = a/as, and Iq is the zeroth order Bessel function. The first and 
second terms in Eq. (52) are due to the Coulomb interaction C, where the exchange term enters with a minus sign. 
The parameter c = y/TT/2(e^ / kub) / fu^o (~ 2.4, for fuvo = 3meV) is the ratio between Coulomb and confining energy. 
The last term comes from the confinement potential W. The result J{B) is plotted in Fig. 8 (dashed line). Note that 
typically | J/ftwol ^ 0.2. Also, we see that J > for S = 0, which must be the case for a two-particle system that 
is time-reversal invariant (Mattis, 1988). The most remarkable feature of J{B), however, is the change of sign from 
positive to negative a.t B = Bl, which occurs over a wide range of parameters c and a. This singlet-triplet crossing 
occurs at about Bl = 1.3 T for tuvo = 3meV (c = 2.42) and d = 0.7. The transition from antiferromagnetic (J > 0) 
to ferromagnetic (J < 0) spin-spin coupling with increasing magnetic field is caused by the long-range Coulomb 
interaction, in particular by the negative exchange term, the second term in Eq. (52). As B ^ Bq (f« 3.5T for 
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FIG. 8 Exchange energy J in units of meV plotted against 
the magnetic field B (in units of Tcsla), as obtained from the 
s-wave Heitlcr-London approximation (dashed line), Eq. (52), 
and the result from the improved sp-hybridized Heitler-London 
approximation (triangles) which is obtained numerically as ex- 
plained in the text. Note that the qualitative behavior of the 
two curves is similar, i.e. they both have zeroes, the s-wave 
approximation at Bl, and the sp-hybridized approximation at 
Bl^ , and also both curves vanish exponentially for large fields. 
Bo = {hujo/ fiB){'m/ms) denotes the crossover field to magneti- 
cally dominated confining {B ^ Bo). The curves are given for 
a confinement energy huJo = 3 meV (implying for the Coulomb 
parameter c = 2.42), and inter-dot distance a = 0.7ob. 



hujQ = 3meV), the magnetic field compresses the orbits by a factor b B/Bq ^ 1 and thereby reduces the overlap 
of the wavefunctions, « exp{—2(f{2b — 1/6)), exponentially strongly. Similarly, the overlap decays exponentially 
for large inter-dot distances, d ^ 1. Note however, that this exponential suppression is partly compensated by the 
exponentially growing exchange term (12|C|21)/S'^ oc exp(2d^(6— 1/6)). As a result, the exchange coupling J decays 
exponentially as exp(— 2d^6) for large b or d, as shown in Fig. 9b for B = (6 = 1). Thus, the exchange coupling J 
can be tuned through zero and then suppressed to zero by a magnetic field in a very efficient way. 



2. Limitations and extensions of HL 

We note that the HL approximation breaks down explicitly (i.e. J becomes negative even when B = 0) for certain 
inter-dot distances if the interaction becomes too strong (c exceeds « 2.8). 

The HL method can be improved by taking into account more than one single-dot orbital. Admixture of higher 
orbitals can be taken into account using a variational approach; the orbitals obtained in this way are termed hybridized 
orbitals, in analogy to hybridized molecular orbitals in chemistry. Some results obtained with sp-hybridized QD 
orbitals are plotted in Fig. 8. 

Another limitation of the HL approximation its restriction to quantum dots that are occupied with a single electron. 
Even with a single orbital, the Pauli principle allows for the presence of a second electron with opposite spin on a 
QD orbital. While this admixture of double occupancy is suppressed by the repulsive Coulomb interaction between 
electrons, it nevertheless plays a relevant role. 



3. The Hund-IVIulliken approach and the Hubbard Limit 

The Hund-Mulliken (HM) approximation (also known as molecular orbit approximation (Mattis, 1988)) extends 

the HL approach by including also the two doubly occupied states, which both are spin singlets (Burkard et al, 
1999a). This extends the orbital Hilbert space from two to four dimensions. First, the single particle states have to 
be orthonormalized, leading to the states $±a = (<^±a ~ 9^^a )/^l-2Sg + g^, where S again denotes the overlap of 
if-a with <^_|_(j and g = {1 — \/l — S'^)/S. Then, diagonalization of 
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in the space spanned by «'^„(ri, ra) ^ $± a(ri)$±a(r2 ), 1'i(ri,r2) = [4>+a(ri)4>_a(r2) ± 4>_a(ri)4>+a(r2)]/V2 yields 
the eigenvalues es± =2e + UH/2 + V+± ^JU^/A + e^o = 2e + Un-'2.X + V+ (singlet), and et = 2e 4- V_ (triplet), 
where the quantities are given in (Burkard et al., 1999a). The exchange energy then becomes 

J = et-e.-=y-^ + i^C/2+164. (54) 

In the standard Hubbard approach for short-range Coulomb interactions (and without B-field) (Mattis, 1988) J 
reduces to —U/2 + y/U"^ -\- 16t'^/2, where t denotes the hopping matrix element, and U the on-site repulsion. Thus, 
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FIG. 9 The exchange couphng J obtained from HM (full line), Eq. (54), and from the extended Hubbard approximation 
(dashed Une), Eq. (55). For comparison, we also plot the usual Hubbard approximation where the long-range interaction term 
V is omitted, i.e. J = At^/Un (dashed-dotted line). In (a), J is plotted as a ftinction of the magnetic field B at fixed inter-dot 
distance {d = o/ob = 0.7), and for c = 2.42, in (b) as a function of inter-dot distance d = a/as at zero field {B = 0), and again 
c = 2.42. For these parameter values, the s-wave Heitler-London J, Eq. (52), and the HM J (full line) are almost identical. 



tn and Uu are the extended hopping matrix element and the on-site repulsion, resp., renormaUzed by long-range 
Coulomb interactions. The remaining two singlet energies es+ and Esq arc separated from et and es_ by a gap of 
order Uu and arc therefore neglected for the study of low-cncrgy properties. Typically, the "Hubbard ratio" t-a/Un 
is less than 1, e.g., if d = 0.7, hcoo = 3meV, and B = 0, we obtain tu/Uu = 0.34, and it decreases with increasing B. 
Therefore, we are in an extended Hubbard limit, where J takes the form 

. = |i + V. (55) 

The first term has the form of the standard Hubbard approximation (Fradkin, 1991) but with tn and Un being 
renormalized by long-range Coulomb interactions. The second term V is new and accounts for the difference in 
Coulomb energy between the singly occupied singlet and triplet states ^j.. It is precisely this V that makes J 
negative for high magnetic fields, whereas t'^/Un > for all values of B (sec Fig. 9a). Thus, the usual Hubbard 
approximation (i.e. without V) would not give reliable results, neither for the B-dependence (Fig. 9a) nor for the 
dependence on the inter-dot distance a (Fig. 9b). Since only the singlet space has been enlarged, it is clear that we 
obtain a lower singlet energy than that from the s-wavc Heitler-London calculation, but the same triplet energy ct, 
and therefore J = et — Cg exceeds the s-wave Heitler-London result, Eq. (52). However, the on-site Coulomb repulsion 
U oc c strongly suppresses the doubly occupied states ^^.^ and already for the value of c = 2.4 (corresponding to 
tiujo = 3meV) we obtain almost perfect agreement with the s-wave Heitler-London result (Fig. 8). For large fields, 
i.e. B ':$> Bq, the suppression becomes even stronger (U oc VB) because the electron orbits become compressed 
with increasing B and two electrons on the same dot are confined to a smaller area leading to an increased Coulomb 
energy. Being a completely orbital effect, the exchange interaction between spins of course competes with the Zeeman 
coupling Hz of the spins to the magnetic field. In our case, however, the Zeeman energy Hz is small and exceeds the 
exchange energy (polarizing the spins) only in a narrow window (about 0.1 T wide) around Bl^ and again for high 
fields (B > 4T). 



4. Numerical work 

While the calculation discussed in Sec. ELB.S above take only the ground-state orbital in each QD into account, the 
HM, like the HL, approximation can be refined to include more orbital levels of the QDs. Such extended calculations 
are usually done numerically, and are very closely related to Hartree-Fock (HF) calculations. Note, however, that HF 
is not sufficient for the purpose of calculating a spin exchange coupling, since it is not capable of including entangled 
(quantum correlated) states such as the spin singlet or m = triplet. This is typically remedied by invoking the 
so-called configuration-interaction method which includes linear superpositions of HF states. Numerical studies of 
the double-dot system with one (Hu and Das Sarma, 2000) and three (Hu and Das Sarma, 2001) electrons per QD 
showed good agreement with the somewhat more crude approximations discussed above. 

At finite magnetic field, the exchange coupling Eq. (52) can be tuned through zero by changing electrostatic 
properties (QD size, distance, electric field). It has been confirmed both numerically and in actual experiment 
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FIG. 10 The exchange couphng J measured as a func- 
tion of the apphcd magnetic field B± using conductance 
spectroscopy in a two-electron dot system defined in 
a GaAs/AlGaAs heterostructure. There are signatures 
that a double has been formed although a single dot 
structure was used in the experiment (Zumbiihl et al., 
2004). The shape of the dot is not circular, but some- 
what elongated. The dot spectra appear to be consis- 
tent with a parabolic potential with harmonic energies 
huja = 1.2 meV and tu vb = 3. 3 meV, corresponding to a 
spatial elongation of ■\/u!b/'^a ~ 1-6- (Figure courtesy of 
D. M. Zumbiihl, Harvard University). 



(Kyriakidis et ai, 2002) that singlet-triplet crossings can be induced in a single QD by changes in the dot potential 
at constant magnetic field. 



5. Measurements of QD exchange 

Signatures of singlet-triplet crossings have been observed using transport spectroscopy in lateral GaAs quantum 
dot structures (Zumbiihl et al., 2004) (see Fig. 10). Although a single dot structure was used, there are signatures 
that a double dot was formed in the experiment (Engcl et al., 2004). 

These data seem to be in rather good qualitative agreement with theory (Burkard et al., 1999a), bearing in mind that 
the absolute magnitude of the exchange coupling J strongly depends on the inter-dot distance which is a free parameter 
of the theory. Similar double-dot experiments with the double-dot systems shown in Fig. 6 are in preparation. 



C. Exchange in vertically coupled QDs 

While lateral QDs are adjacent to each other in a two-dimensional electron gas, vertically coupled QDs are stacked 
on top of each other in a three-dimensional semiconductor structure. Vertical coupling occurs both in QDs etched 
out of multilayer structures and electrically gated (Austing et al., 1998) and in self-assembled QDs originating from 
the Stranski-Krastanov growth (Fafard et al., 1999; Fricke et al., 1996; Luyken et al., 1998). A system of vertically 
coupled QDs is shown schematically in Fig. 11. There are a number of works on the exchange coupling between spins 
located in this type of coupled QDs (Burkard et al, 2000c; Imamura et al, 1998, 1996, 1999: Oh et al, 1996; Tokura 
et al., 1999). This Section in devoted to the exchange coupling between spins in vertically coupled dots under the 
influence of both in-plane magnetic and electric fields, B|| and E\\, and perpendicular fields B^, E± (Burkard et al., 
2000c). Electronic spectra and charge densities for two electrons in a system of vertically tunnel-coupled QDs at 
zero magnetic field, B = 0, were calculated in (Bryant, 1993). Singlet-triplet crossings in the ground state of single 
(Wagner et al., 1992; Wojs et al., 1996) and coupled dots with two (Oh et al., 1996) to four (Imamura et al., 1998, 
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FIG. 11 (a) Sketch of the vertically coupled double 
quantum-dot system. The two dots may have different 
lateral diameters, aB+ and as-. We consider magnetic 
and electric fields applied either in-plane (By, Sy) or per- 
pendicularly {B±, E±). (b) The model potential for the 
vertical confinement is a double well, which is obtained 
by combining two harmonic wells at z = ±a. 
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1996, 1999; Tokura et ai, 1999) electrons in vertically coupled dots in the presence of a magnetic field perpendicular 
to the growth direction {B± in Fig. 11) have been predicted. In order to apply the HL and HM methods introduced 
in Sec. III.B to vertically coupled QDs, a Hamiltonian of the same form as Eq. (45) can be used, but with the 
single-particle Hamiltonian replaced by (Burkard et ai, 2000c) 



Hr, p) = ^ (p - ^A(r)) \ezE+ Vi{t) + K(r). 



(56) 



Also note that r is a three-dimensional vector here, as opposed to a two-dimensional vector in Sec. III.B. The 
potential Vi in h describes the lateral confinement, whereas Vy models the vertical double-well structure. For the 
lateral confinement the parabolic potential 

2 [ ag_(a;^ -l-y^), z <0, 

is chosen, where we have introduced the anisotropy parameters ao± determining the strength of the vertical relative 

to the lateral confinement. In the presence of a magnetic field B± perpendicular to the 2DES, the one-particle 
problem has the Fock-Darwin states (Darwin, 1930; Fock, 1928) as an exact solution. Furthermore, it has been shown 
experimentally (Fricke et ai, 1996) and theoretically (Wojs et ai, 1996) that a two-dimensional harmonic confinement 
potential is a reasonable approximation to the real confinement potential in a lens-shaped SAD. In describing the 
confinement Vy along the inter-dot axis, a (locally harmonic) double well potential of the form (see Fig. lib) 



K = ^(.^-aT, (58) 

can be used; in the limit of large inter-dot distance o » ob, the potential Vy in the vicinity of z Ki ±o becomes a 
harmonic well of frequency lUz- Here a is half the distance between the centers of the dots and aB = yjfij (mujz) is 
the vertical effective Bohr radius. For most vertically coupled dots, the vertical confinement is determined by the 
conduction band off'set between different semiconductor layers; therefore in principle a square-well potential would 
be a more accurate description of the real potential than the harmonic double well (note however, that the required 
conduction-band offsets are not always known exactly). There is no qualitative difference between the results presented 
below obtained with harmonic potentials and the corresponding results obtained using square-well potentials (Seelig, 
1999). 



2 



1. Perpendicular Magnetic Field B± 

For a magnetic field B = B± (cf. Fig. 11) and E = 0, one obtains the ground-state Fock-Darwin (Darwin, 1930; 
Fock, 1928) solution 



<pUx,y,z) = (!^)^/^^e— (-(^^+^^)+(^=F«)^)/2^ 



(59) 



corresponding to the ground-state energy e± = /k<;^(l -|- 2a±)/2. In Eq. (59) we have introduced a±(B) = 



*^o± + ^l{B)^/uj'^ = yJ(y-Q± + B^/Bq, with u)l{B) = eB/2mc the Larmor frequency and Bq = 2maOz/e the mag- 
netic field for which lOz = ujl- The parameters a±{B) describe the compression of the one-particle wave function 
perpendicular to the magnetic field. The HL with this Hamiltonian yields 



2^2 



z^e^^'^' (l_erf(dv^))-^ — 



yjl- (a+ + a- - 1)' 



: arccos(a+ -|- a_ — 1) 



+ 4 ("0+ - "0-) 



a+ — Qf- 



(l-erf(d)) + -(l + d2) 



(60) 



where erf (a;) denotes the error function. We have introduced the dimensionless parameters d = a/as iov the inter-dot 
distance, and c = ■^/Tr/2{e'^ / kqb) /hi^z for the Coulomb interaction. Note that a±, fi = 2a+a_/ (a+ -|- a_), and the 
overlap 



S = 2^^ ' 



exp(-rf ), 



(61) 
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FIG. 12 Exchange energy J (left graph) and single-electron tunneling amplitude t (right graph) as a function of the applied 
magnetic field for two vertically coupled small (height 6nm, width 12 nm) InAs (m = O.OSme, k = 14.6) self-assembled QDs 
in a ccntcr-to-ccntcr distance of 9nm (d = 1.5). The box-shaped symbols correspond to the magnetic field B± applied in z 
direction, the circle symbols to the field -B|| in x direction. The plotted results were obtained using the HM method and are 
reliable up to a field Bo ~ 15 T where higher levels start to become important. 



depend on the magnetic field B. The first term in the square brackets in Eq. (60) is an approximate evaluation of the 
direct Coulomb integral (12|C|12) for d > 0.7 and for magnetic fields B < Bq- The second term in Eq. (60) is the 

(exact) exchange Coulomb integral (12|C|21)/5^, while the last two terms stem from the potential integrals, which 
were also evaluated exactly. If the two dots have the same size, the expression for the exchange energy Eq. (60) can 
be simplified considerably. 

For two vertically coupled dots of equal size, we set ao+ = ckq- = ao in Eq. (60) and using Eq. (61), we obtain 

cy^e^"'^' (l-erf (dV2^)) - - , arccos(2a - 1) ^ (l d^) , (62) 

V V // TT .^1 — (2a — 1)2 4 



J = 



sinh(2(i2) 



where a = y^aQ + B'^/B^. As before, the first term in Eq. (62) is the direct Coulomb term, while the second term 
(appearing with a negative sign) is the exchange Coulomb term. Finally, the potential term in this case equals 
W = (3/4) (1 -|- cP) and is due to the vertical confinement only. For two dots of equal size neither the prefactor 
2S^/ (1 — 5^) nor the potential term depends on the magnetic field. Since the direct Coulomb term depends on B± 
only weakly, the field dependence of the exchange energy is mostly determined by the exchange Coulomb term. The 
exchange coupling can also be calculated using the HM method (Burkard et ai, 2000c). 

The dependence of the exchange energy J on an electric field E± applied in parallel to the magnetic field, i.e. 
perpendicular to the xy plane, withing the HL approximation, was found to be 

J(B,E^)^JiB,0) + hw,^—^-i^-^j , (63) 

where Eq = mcu'^/eaB- The growth of J is thus proportional to the square of the electric field E±, if the field is 
not too large. This result is supported by a HM calculation, yielding the same field dependence at small electric 
fields, whereas if eE±a is larger than Uh, double occupancy must be taken into account. The electric field causes 
the exchange J at a constant magnetic field B to cross through zero from J{E = 0, B) < to J > 0. This effect is 
signaled by a change in the magnetization M (Burkard et al., 2000c). 



2. In-plane magnetic field By 

In this section we consider two dots of equal size in a magnetic field which is applied along the a::-axis, i.e. 
in-plane (see Fig. 11). Since the two dots have the same size, the lateral confining potential Eq. (III.B) reduces 
to V{x,y) = mcij^«o(^^ + where the parameter ao describes the ratio between the lateral and the vertical 

confinement energy. The vertical double-dot structure is modeled using the potential Eq. (58). The situation for an 
in-plane field is a bit more complicated than for a perpendicular field, because the planar and vertical motion do 
not separate. In order to find the ground-state wave function of the one-particle Hamiltonian /i^.^,, an approximate 
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variational method can be applied (Burkard et al., 2000c) with the result 

/mujz \ I , ^1 
ba = ) ("oap)'' exp 



'P±a 



2h 



(64) 



where the parameters a{B) = a/oq + {B/BqY and /3(i?) = ^1 + {B/BqY, describing the wave- function compression 
in y and z direction, respectively, have been introduced. 
The resulting exchange coupling J in this case is 



J{B,d) = Jo{B,d) - h(j^ 



13 -a 



B 
B^ 



(65) 



where Jq denotes the expression from Eq. (51). The variation of the exchange energy J as a function of the magnetic 
field B is, through the prefactor 2S'^/{1 - S"'), determined by the overlap S{B, d) = exp[-d2(/3(B) + {B / Bof) / a[B)], 
depending exponentially on the in-plane field, while for a perpendicular field the overlap is independent of the field 
(for two dots of equal size), see Eq. (61). The HL result can again be improved by performing a molecular-orbital 
(HM) calculation of the exchange energy, which we plot in Fig. 12 (left graph, circle symbols). 



3. Electrical switching of the interaction 

Operating a coupled quantum dot as a quantum gate requires the ability to switch on and off the interaction 

between the electron spins on neighboring dots. A simple method of achieving a high-sensitivity switch for vertically 
coupled dots involves a horizontally applied electric field E\\ . The idea is to use a pair of QDs with different lateral 
sizes, e.g. a small dot on top of a large dot (00+ > Qo-. see Fig. 11). Note that only the radius in the xy plane has 
to be different, while it can be assumed that the dots have the same height. Applying an in-plane electric field 
in this case causes a shift of the single-dot orbitals by /\x± = eE^^/m.uj'^aQj^ = £'||/£'o«o±' where Eq = hujz/eaB, see 
Fig. 13. It is clear that the electron in the larger dot moves further in the (reversed) direction of the electric field 
(A.T_ > Aa;+), since its confinement potential is weaker. As a result, the mean distance between the two electrons 
changes from 2d to 2d', where 



d' = ^d'' + ^{Ax--Ax+f = ^d^+A''(^^^ , (66) 

with A = (l/aQ_ — l/aQ_|_)/2. Using Eq. (61), we find that the wavefunction overlap scales as S* oc cxp(— rf'^) cx 
exp[— ^^(iSy/i^o)^]. Due to this high sensitivity, the electric field is an ideal "switch" for the exchange coupling J 
which is (asymptotically) proportional to 5^ and thus decreases exponentially on the scale Eo/2A. Note that if the 
dots have exactly the same size, then A = and the effect vanishes. An estimate of J as a function of £"11 can be 
obtained by substituting d' from Eq. (66) into the HL result, Eq. (60). A plot of J(£^||) obtained in this way is shown 
in Fig. 13 for a specific choice of GaAs dots. Note that this procedure is not exact, since it neglects the tilt of the 
orbitals with respect to their connecting line. Exponential switching is highly desirable for quantum computation, 
because in the "off" state of the switch, fluctuations in the external control parameter (e.g. the electric field i?||) 
or charge fluctuations cause only exponentially small fluctuations in the coupling J. If this were not the case, the 
fluctuations in J would lead to uncontrolled coupling between qubits and therefore to multiple-qubit errors. Such 
correlated errors cannot be corrected by known error-correction schemes, which are designed for uncorrelated errors 
(Preskifl, 1998b). 



D. Single-qubit operations 

Single-qubit operations with the Hamiltonian Eq. (10) require a time- varying Zeeman coupling (^/i^S • B)(t) 
(Burkard et al., 1999a; Loss and DiVincenzo, 1998), which can be controlled by changing the magnetic field B or the 
g- factor g. Effective magnetic fields/g-factors can be produced by coupling the spin via exchange to a ferromagnet 
(Loss and DiVincenzo, 1998) or to polarized nuclear spins (Burkard et al, 1999a). We review here how the g-factor 
of an electron in a semiconductor heterostructure can be modulated by shifting its orbital between layers of host 
material with different g-factors (DiVincenzo et al., 1999; DiVincenzo and Loss, 1999). 

The spin-orbit coupling can leads to large deviations of the Lande g-factor (both in the positive and negative 
direction) in bulk semiconductors from the free-electron value go = 2.0023. The effective g-factors in these materials 
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FIG. 13 Controlled switching of the exchange coupling J between 
dots of different size by means of an in-plane electric field -Ey at 

zero magnetic field, 5 = 0. The coupling is "on" at i? = 0. When 
is applied, the larger dot is shifted to the right by Aa;_, whereas 
the smaller dot is shifted by Aa;+ < Aa;_, where Ax± = E^^/EoOQ-j. 
and Eo = hoJz/eaB- With _E|| increasing, J decreases exponentially, 
J Rj Ri exp[-2^^(_E||/So)^]. The parameters used for this plot 
arc hujz ~ 7mcV, d = 1, ao+ — 1/2 and ao- = 1/4, yielding 
Eo — hu>z/eaB = 0.56mV/mn and A = (qo+ — ao_)/2Qo+ao- = 
6. The coupling J decreases exponentially on the scale Eo/2A = 
0.047 mV/nm for the electric field. 
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FIG. 14 Effective g- factor Qes of electrons confined in 
a AlxGai-x As-Ga As-In^ Al^Gai-x-j/ As- Ala, Gai - ^ As 
heterostructure (a; = 0.3, y = 0.1) as a function 
of the applied electric field E in growth direction. 
The widths of the quantum well and the barriers are 
w = wb = 10 nm. The g-factors which are used for the 
materials are indicated with dashed horizontal lines. 



range from large negative to large positive numbers. In confined structures such as quantum wells, wires, and dots, 
the g-factor is modified with respect to the bulk value and sensitive to an external bias voltage (Ivchenko et al., 
1997). In the case of a layered structure, the effective g-factor of electrons can be varied by electrically shifting their 
equilibrium position from one layer (with g-factor gi) to another (with another g-factor g2 Qi)- The bulk g-factors of 
the layer materials and linear interpolations between them, have been used here as an approximation which becomes 
increasingly inaccurate as the layers become thinner (Kiselev et al., 1998). 

Let us assume that by replacing some fraction y of Ga atoms in the upper half of a AlGaAs-GaAs- AIGaAs quantum 
well by In atoms (we have used y = 0.1) we obtain the following layered heterostructure: 

Alj;Gai_a;As-GaAs-InyAla;Gai_a;_yAs-Ala,Gai_xAs, 

where x denotes the Al content in the barriers (typically around 30%). In such a structure, the effective g-factor can be 

modified by changing the vertical position of the electrons via top or back gates. If the electron is mostly in a pure GaAs 
environment, then its effective g-factor will be around the GaAs bulk value (gcaAs = —0.44) whereas if the electron is 
in the InAIGaAs region, the g-factor will be more negative due to the large negative InAs value ((/inAs = —15). The 
one-dimensional problem of one electron in such a structure has been analyzed numerically. When the effective mass 
m{z) is spatially varying, the Hamiltonian in the effective mass approximation can be written as 



' ' Viz) 



dz 2m{z) dz 



^{z) = E^{z). (67) 



This problem can be discretized in real space and subsequently diagonalized numerically (DiVincenzo et al., 1999). 
Finally, the effective g-factor is calculated by averaging the local g-factor g{z) over the electronic density in the 
ground-state (see Fig. 14), 

5eff = J dzg{z)\^iz)\\ (68) 

The option of performing single-qubit rotations by electrostatically controlling the g-factor makes all-electric control of 
a spin-based quantum computer (an array of QDs as in Fig. 1) possible and thus offers a way around the problematic 
local magnetic field implementation of single-qubit gates. Another method to circumvent single-spin operations 
completely (however, at a higher cost of gates and exchange operations) is the exchange-only architecture outlined in 
Sec. II.C. 
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FIG. 15 Left: QDs embed- 
ded inside a microdisk struc- 
ture, from (Imamoglu et ai, 
1999). Each QD is ad- 
dressed selectively by a laser 
field from a fiber-tip. The 
laser frequencies are chosen 
to select out the pair of QDs 
that will participate in gate 
operation. All dots strongly 
couple to a single cavity- 
mode. Right: Energy lev- 
els of a III-V (or II- VI) semi- 
conductor QD. It is assumed 
that confinement along the 
z-direction is strongest. 



E. Semiconductor microcavities 

Here, we present a modification of the Loss-DiVincenzo scheme (Sec. II. A) for QC based on QD electron spins. In 
contrast to the original scheme (Loss and DiVincenzo, 1998), where the spins are coupled via direct exchange, this 
couphng which is mediated through a single microcavity mode and uses laser fields to mediate coherent interactions 
between distant QD spins (Imamoglu et ai, 1999). 

The cavity scheme is shown in Fig. 15: the doped QDs are embedded in a microdisk structure with diameter 
d ~ 2 /im and thickness d ~ 0.1 /im. Experiments have shown that InAs self-assembled QDs can be embedded in 
microdisk structures with a cavity quality factor Q ~ 12000 (Gerard and Gayral, 1999). It is assumed that the QDs 
are designed such that the quantum confinement along the z-direction is the strongest. The in-plane confinement is 
also assumed to be large enough to guarantee that the electron will always be in the ground-state orbital. Because 
of the strong z-axis confinement, the lowest energy eigenstates of such a III-V or II- VI semiconductor QD consist 
of \mz = ±1/2) conduction-band states and jm^ = ±3/2) valence-band states. The QDs are doped such that each 
QD has a full valence band and a single conduction band electron: we assume that a uniform magnetic field along 
the x-direction (i?^) is applied, so the QD qubit is defined by the conduction-band states \mx = —1/2) = |i) and 
|m, = 1/2) = It) (Fig. 15, right). 



1. Single-qubit operations 

Single-bit operations are carried out in this scheme by applying two laser fields E^^xit) and EL_y(t) with Rabi 
frequencies ^l,x and ^L,y, and frequencies ujl,x and ujL,y (polarized along the x and y directions, respectively) that 
exactly satisfy the Raman-resonance condition between ||) and |t). The laser fields are turned on for a short time 
duration that satisfies a 7r/r-pulse condition, where r is any real number. The process can be best understood as a 
Raman 7r/r-pulse for the iioie in the conduction band state. The laser field polarizations should have non-parallel 
components in order to create a non-zero Raman coupling (if there is no heavy- hole light-hole mixing). These arbitrary 
single-bit rotations can naturally be carried out in parallel. In addition, the QDs that are not doped by a single electron 
never couple to the Raman fields and can safely be ignored. 



2. Two-qubit operations 

Two-qubit operations are mediated by virtual photons that are emitted to and reabsorbed from the microcavity 
field. It is assumed that the x-polarized cavity-mode with energy Wcav (S- = 1) and a laser field (assumed to be y- 
polarized) establish the Raman transition between the two conduction-band states, in close analogy with the atomic 
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cavity-QED schemes (Pellizzari et al, 1995). The Hamiltonian for a single QD is written as H = Hq + Hmt, with 

Ho= ^ w^e^e^ +a;cavaLv<^cav +<^La^aL, (69) 

(T=T,i>±3/2 

where e|, annihilate an electron with spin |, I along the x direction in the conduction band and e±3/2 annihilates 
an electron with spin ±3/2 along the z direction in the valence band, cf. Fig. 15 (right). The interaction can be 
written as 

-f^int = 9 (a+el3/2e-i/2 - at4/2^i/2 + ^-C-) > (70) 

where the operators for the circularly polarized light are expressed in terms of the x-polarized cavity mode and the 
y-polarized laser field, a± = (ccav ± iaL)/V2, and the conduction-band operators in the z basis can be expressed in 
terms of those in the x basis, e±i/2 = (ej, ± e-f)/\/2. With LO-3/2 = W3/2 = w„ and the definition = (e_3/2 — 63/2)72, 
the following result for the cavity Hamiltonian is obtained, 

-ffint = 9 (aLveJei + e|e„acav) - 19 (ot4ei - eje^ai,^ . (71) 

The valence band states are eliminated by a Schrieffer- Wolff transformation (Madelung, 1978; Schrieffer and Wolff, 
1966), iJeff = e-^He^, with 

S = (al^^ele-f - eje^acav) - i-^^ {a{elei + eje^ai) , (72) 

where Awf = — w^, — Wcav and AuJi = loi—uj^—ujl. Neglecting all terms 0{g^) and replacing e|e„ by its expectation 
value {e\ey) = 1 and by fii exp(— iw^t) one obtains the effective Hamiltonian 



HeS = WcavClLv^icav + ^ 



■'T """i 



aLv^ite-'"--* - h.c. 



, (73) 



where the sum runs over all QDs of the system, gig is the effective 2-photon coupling coefficient, (J|| = |t)(i| the 
spin projection operator for the i-th QD, and = u'^^ ^ "^l' exact two-photon-resonance condition would be 
Acj| = cj| — — u>ca,v = Aa;| = — — w}^. The derivation of i?eff assumes Aw^ Scav) 'a'ij^ ^ ksT, and 
'^IX ^ .9off > Tcav, where Fcav denotes the cavity decay rate (not included in Eq. (73)). The third and fourth terms 
of Eq. (73) describe the ac-Stark-effect caused by the cavity and laser fields, respectively. 

In order to implement a CNOT quantum gate, one would turn on laser fields lo\ and wj^ to establish near two-photon 
resonance condition for both the control (i) and the target (j) qubits, 

Aj = U)\^ - LOcav +I^L = ^3 " ^'^^^ 

If the two- photon detunings A^ are chosen large compared to the cavity linewidth and glg{t), the cavity modes can 
be eliminated with a second Schrieffer- Wolff transformation to obtain an effective two-qubit interaction Hamiltonian 
in the rotating frame (interaction picture with Hq = 

where gij{t) = 9l([{t)glfi{t)/ Ai and Aij = Aj — A^. The implementation of the conditional phase-flip (CPF) and the 
CNOT or quantum XOR gates between two spins i and j from a transversal (XY) spin coupling of the form Eq. (76) 
has been discussed in Sec. II. B. 

(2) 

The interaction Hamiltonian describes the coupling of the QD spins via the following virtual process. One of 
the QDs emits a virtual photon into the cavity while absorbing a laser photon. The cavity photon is then reabsorbed 

by the other QD while a laser photon is emitted. Due to the spin splitting in the QD spectrum, Fig. 15 (right), this 

(2) 

process is spin sensitive and leads to the spin-spin coupling H^^^ between the QDs. 
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3. Measurement 



In the cavity QED scheme, measurement of a single QD spin can be achieved by applying a laser field -Bl.j, to the 
QD to be measured, in order to realize exact two-photon resonance with the cavity mode. If the QD spin is in state 
I J,), there is no Raman couphng and no photons will be detected. If on the other hand, the spin state is |t), the 
electron will exchange energy with the cavity mode and eventually a single photon will be emitted from the cavity. 
A single photon detection capability is thus sufficient for detecting a single spin. 

4. Related proposals 

A related proposal is to use optically controlled virtual excitations of delocalized excitons as mediators of RKKY 
type spin interaction between electrons localized in neighboring QDs (Piermarocchi et al., 2002). The spin interaction 
in this case is isotropic, H = JSj • Sj . 

F. Decoherence 

The spin coherence time in semiconductors — the time over which the phase of a superposition of spin-up and spin- 
down states, Eq. (1), is well-defined —can be much longer than the charge coherence time (a few nanoseconds). In 
fact it is known from experiment that they can be orders of magnitude longer. This is of course one of the reasons for 
using spin as a qubit (Loss and DiVincenzo, 1998) rather than charge. In bulk GaAs and in CdSe quantum dots, the 
ensemble spin coherence time T2 , being a lower bound on the single-spin decoherence time T2, was measured using a 
technique called time-resolved Faraday rotation (Gupta et at, 1999; Kikkawa et al, 1997). For a detailed account of 
these experiments, we refer the reader to Chaps. 4 and 5 of (Awschalom et al., 2002). The spin relaxation time Ti in 
a single-electron QD in a GaAs heterostructure was probed via transport measurements and found to approach one 
microsecond (Hanson et al., 2004, 2003). It has been proposed to also measure the single-spin T2 in such a structure in 
a transport experiment by applying electron spin resonance (ESR) techniques (Engcl and Loss, 2001). In this scheme, 
the stationary current exhibits a resonance whose line width is determined by the single-spin decoherence time T2. 

Below, a number of decoherence mechanisms for spin in semiconductor nanostructures will be listed. It Should be 
emphasized, though, that it is usTially hard for theory to predict which mechanism is dominant. Nevertheless, the 
understanding of the underlying mechanisms for a list of possible causes can be a very valuable tool for the purpose 
of achieving long coherent operation in a future quantum device. 

1. Phonons and the spin-orbit coupling 

Phonon-assisted transitions between different discrete energy levels (or Zeeman sublevels) in GaAs quantum dots 

can cause spin flips and therefore spin decoherence (Khaetskii, 2001; Khaetskii and Nazarov, 2000, 2001). There are 
various mechanisms originating from the spin-orbit coupling which lead to such spin flip processes; the most effective 
mechanisms in 2D have to do with the broken inversion symmetry, either in the elementary crystal cell or at the 
heterointerface. The spin-orbit Hamiltonian for the electron in such a structure is given by Eq. (22). The relaxation 
rates F = T{~^ are evaluated in leading perturbation order in this coupling, with and without a magnetic field. The 
spin-orbit coupling Hgo mixes the spin-up and spin-down states of the electron and leads to a non-vanishing matrix 
element of the phonon-assisted transition between two states with opposite spins. However, one of the main findings 
of (Khaetskii, 2001; Khaetskii and Nazarov, 2000, 2001) is that the spin relaxation of the electrons localized in the 
dots differs strongly from that of delocalized electrons. It turns out that in quantum dots (in contrast to extended 
2D states), the contributions to the spin-flip rate proportional to are absent in general. This greatly reduces 
the spin-flip rates of electrons confined to dots. The finite Zeeman splitting in the energy spectrum also leads to 
contributions oc 




(77) 



where ^^.^o is the orbital energy level splitting in the QD and Tq{B) is the inelastic rate without spin flip for the 
transition between neighboring orbital levels. 
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Spin-flip transitions between Zeeman sublevels occur with a rate that is proportional to the fifth power of the 

Zeeman splitting, 



h(fiLUo) 



- f^A,. (78) 



The dimensionless constant Ap oc characterizes the strength of the effective spin-piezo-phonon coupling in the 
heterostructure and ranges from « 7 • 10~^ to « 6 • 10~^ depending on (3. To give a number, « 1.5 • 10^ s~^ for 
HuiQ = lOK and at a magnetic field i? = 1 T. 

It can be shown that under realistic conditions, a general symmetry argument leads to the conclusion that the spin 
decoherence time T2 does not have a transverse contribution (in leading order), in other words, T2 = 2Ti for spin-orbit 
(phonon) related processes (Golovach et al., 2004). 



2. Nuclear spins 

The nuclear spins of the host material can cause dccohcrcncc via spin flips that arc caused by the hypcrflnc 
interaction. A rough perturbativc estimate of this effect (Burkard et at, 1999a) suggests that the rate of such processes 
can be suppressed by either polarizing the nuclear spins or by applying an external magnetic field. The suppression 
factor is {B*/B)'^/N, where B* = AI / gji-Q is the maximal magnitude of the effective nuclear field (Ovcrhauscr field), 
N the number of nuclear spins in the vicinity of the electron, and A the hyperfine coupling constant. In GaAs, the 
nuclear spin of both Ga and As is / = 3/2. The field B denotes either the external field, or, in the absence of an 
external field, the Overhauser field B = pB* due to the nuclear spin polarization p, which can be obtained e.g. by 
optical pumping (Dobers et al, 1988) or by spin-polarized currents at the edge of a 2DEG (Dixon et al, 1997). In 
the latter case, the suppression of the spin flip rate becomes l/p^N . 

A more detailed analysis treats a single electron confined to an isolated QD under the influence of the hyperfine 
interaction with the surrounding nuclei (Khaetskii et al., 2002). It turns out that the electron spin decoherence time 
T2 is shorter than the nuclear spin relaxation time T„2 determined by the dipole-dipole interaction between nuclei, and 
therefore the problem can be considered in the absence of the nuclear dipole-dipole interaction. Since the hyperfine 
interaction depends on the position via a factor |V'(r)P where V(r) is the electron wavefunction, the value of the 
hyperfine interaction varies spatially. It turns out that this is the relevant cause of decoherence. The analysis is 
complicated by the fact that in a weak external Zeeman field (smaller than a typical fluctuating Overhauser field seen 
by the electron, ~ 100 Gauss in a GaAs QD), the perturbativc treatment of the electron spin decoherence breaks 
down and the decay of the spin precession amplitude is not exponential in time, but either described by a power law, 
l/t*^/^ (for finite Zeeman fields) or an inverse logarithm, l/(lnt)''/^ (for vanishing fields). 

The decay rate 1/72 is thus roughly given by A/hN, where A is the hyperfine interaction constant, and N is the 
number of nuclei within the dot, with N typically 10^. This time is of the order of several /zs. However, it needs to 
be stressed that there is no simple exponential decay which, strictly speaking, means that decoherence cannot simply 
be characterized by the decay times Ti and T2 in this case. The case of a fully polarized nuclear spin state was 
solved exactly in (Khaetskii et al, 2002). The amplitude of the precession which is approached after the decay, is 
of order one, while the decaying part is 1/A'', in agreement with earlier results (Burkard et al., 1999a), see above. A 
large difference between the values of T2 (decoherence time for a single dot) and (dephasing time for an ensemble 
of dots), i.e. T2 ^ T2 is found and indicates that it is desirable to have direct experimental access to single spin 
decoherence times. 

The non-Markovian dynamics of a localized electron spin interacting with an environment of nuclear spins with 
arbitrary polarization p was calculated in (Coish and Loss, 2004) from a perturbative analysis of the generalized 
master equation for the longitudinal and transverse components of the electron spin. 



IV. SUPERCONDUCTING MICRO-CIRCUITS 
A. Overview 

Roughly speaking, three prototypes of superconducting (SC) qubits are studied experimentally. We only briefly 
review them here, and refer the reader to (Makhlin et al, 2001) for a comprehensive review. The charge qubit (Averin, 
1998; Makhlin et al, 1999; Nakamura et al, 1999; Pashkin et al, 2003; Shnirman et al, 1997; Vion et al, 2002), 
operating in the regime Ec ^ Ej, and the flux qubit (Chiorescu et al, 2003; Mooij et al, 1999; Orlando et al, 1999; 
van der Wal et al, 2000), operating in the regime Ej ^ Eq, are distinguished by their Josephson junctions' relative 
magnitude of charging energy Ec and Josephson energy Ej. A third type, the phase qubit (loffe et al, 1999; Martinis 
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et al, 2002), operates in the same regime as the flux qubit, but is represented purely in the SC phase and is not 
associated with any magnetic flux or circulating current. The Josephson phase qiibit consists of a single Josephson 
junction (Martinis et ai, 2002). In flux qubits, the quantum state of the SC phase differences across the Josephson 
junctions in the circuit contain the quantum information, i.e., the state of the qubit. A micrograph of the circuit for 
a SC flux qubit studied in (Chiorescu et al, 2003) is shown in Fig. 16. In charge qubits, the quantum state of the 
charge on SC islands contains the quantum information. 

Both charge and flux qubits have been described by an approximate pseudo-spin Hamiltonian of the type (Makhlin 
et al, 2001), 

H=^<7, + (79) 

where A denotes the tunnel coupling between the two qubit states |0) and |1) (eigenstates of az) and e the bias 
(asymmetry). In Sec. IV. C, a more general model, including the full Hilbert space of a SC circuit, will be discussed. 



B. Decoherence, visibility, and leakage 

1. Decoherence 

Decoherence within the model Eq. (79) can be understood phenomenologically as follows. If the qubit is initially 

prepared in state |1). then it will undergo free Larmor oscillations with frequency v = h'^\/ + e^. Ideally, the 
probability for finding the qubit in state |1) after time t would be a perfect cosine function of t. This ideal Larmor 
precession is shown as a thin dotted line in Fig. 17. Such a Larmor precession experiment (also known as Ramsey 
fringe experiment) determines how well the qubit satisfies item I. A. 3 of DiVincenzo's five criteria. Decoherence is a 
process in which the amplitude of the oscillations decays over time, as shown by the thick solid line in Fig. 17. This 
decay is often (but not always) exponential with a characteristic decoherence time T2. 

All types of SC qubits suffer from decoherence that is caused by a several sources. Decoherence in charge qubits has 
been investigated using the spin-boson model in (Makhlin et al, 2001; Makhlin and Shnirman, 2004). In flux qubits, 
the Johnson- Nyquist noise from lossy circuit elements (e.g., current sources) has been identified as one important 
cause of decoherence (Tian et al, 2000, 2002; van der Wal et al, 2003; Wilhelm et al, 2003). A systematic theory 
of decoherence of a qubit from such dissipative elements, based on the network graph analysis (Devoret, 1997) of the 
underlying SC circuit, was developed for SC flux qubits (Burkard et al, 2004b), and applied to study the effect of 
asymmetries in a persistent-current qubit (Burkard et al, 2004a). The circuit theory for SC qubits will be discussed 
further below in Sec. IV. C. For the Josephson phase qubit (Martinis et al, 2002), decoherence due to bias noise and 
junction resonators was studied in (Martinis et al, 2003; Simmonds et al, 2004). 



2. Visibility 

A different type of imperfection that typically affects SC qubits is a limited visibility v. This effect means that the 
maximum range v of the read-out probability of the qubit being in state |1) is smaller than one. This means, e.g., 
that the probability p{0) of measuring the qubit in state |1) right after preparation in this state is less than one. In 
the case of a symmetric reduction of the visibility, the relation is p(0) = (1 + v)/2. This effect is schematically shown 
in Fig. 17. 



FIG. 16 Electron micrograph of the SC flux qubit circuit 
studied in (Chiorescu et al, 2003). The logical qubit ba- 
sis states correspond to circulating SC currents in the 
smaller loop as indicated. The bright areas are the Al 
wires; the double-layer structure from the shadow evap- 
oration deposition is clearly visible. (Figure courtesy of 
I. Chiorescu and J. E. Mooij, TU Delft). 
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FIG. 17 Theoretical Larmor precession (Ramsey 
fringe) curve with decoherence time r2 < oo and 
limited visibility t; < 1 (solid thick line) , compared 
to the ideal curve (dotted thin line). The proba- 
bility p{t) to find the qubit in state |1) is plotted 
as a function of the free evolution time t. The 
Larmor frequency of the coherent oscillations is 
denoted with v. The visibility v is the maximum 
range of p{0) whereas the decoherence time T2 is 
the time over which the oscillations are damped 
out (in the case of an exponential decay) . For this 
plot, we have chosen T2 — 20 /f and v — 70 %. 

Vt 

3. Leakage 

Limitations of the visibility are often attributed to a mechanism called leakage. Since the SC phase is a continuous 
variable as, e.g., the position of a particle, superconducting qubits (two- level systems) have to be obtained by trunca- 
tion of an infinite-dimensional Hilbert space. This truncation is only approximate for various reasons; (i) because it 
may not be possible to prepare the initial state with perfect fidelity in the lowest two states, (ii) because of erroneous 
transitions to higher levels (leakage effects) due to imperfect gate operations on the system, and (iii) because of erro- 
neous transitions to higher levels due to the unavoidable interaction of the system with the environment. Apparent 
leakage effects may occur if the read-out process is not 100% accurate. Leakage effects due to the non-adiabaticity of 
externally applied fields were studied in (Fazio et ai, 1999). Recent work (Meier and Loss, 2004) shows that leakage 
in microwave-driven Josephson phase qubits leading to a reduced visibility can occur, even if the microwave source is 
pulsed slowly. 

C. Circuit theory 

A recently developed method for deriving the Hamiltonian of SC circuits from their classical dynamics, combined 
with the theory of dissipative quantum systems, can be utilized to describe decoherence in arbitrary SC circuits 
(Burkard et ai, 2004b). 

There exists a variety of theoretical models for dissipative environments in general, and dissipative electrical circuit 
elements (impedances) in particular. A dissipative (resistive) element can be modeled as a transmission line (Werner 
and Drummond, 1991; Yurke, 1987; Yurke and Denker, 1994), i.e. an infinite set of dissipation- free elements (capacitors 
and inductors), or, alternatively, within the widely known Caldeira-Leggett model (Caldeira and Leggett, 1983; Leggett 
et ai, 1987; Weiss, 1999) as a continuum of harmonic oscillators that is coupled to the degrees of freedom of the system 
(in this case, the SC circuit). In the following, the Caldeira-Legget approach will be used. 

The systematic derivation of the dynamical equations for a general (classical) electric circuit is a well-known problem 
in electric engineering that has been tackled using the elegant and convenient network graph analysis methods (Peikari, 
1974). It has been suggested early on that these methods may also be used for a description of the dissipative quantum 
dynamics of superconducting circuits (Devoret, 1997). We shall now explain the circuit graph analysis applied to SC 
qubits, both of the flux (Sec. IV. D) and charge (Sec. IV. E) type. 

D. Flux qubits 

In this Section, the results of the circuit theory for flux qubits are presented; for a derivation, see (Burkard et aL, 
2004b). The IBM qubit (Koch et at, 2003) will be used as a first example, and then followed by other examples. The 
IBM qubit is described by the electrical circuit drawn in Fig. 18. 
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1. The network graph 

As a first step in the circuit analysis of a SC flux qubit, the network graph of the SC circuit is drawn and labeled. In 
the graph, each two-terminal element (Josephson junction, capacitor, inductor, external impedance, current source) 
is represented as a branch connecting two nodes. In Fig. 18 (left panel), the IBM qubit is represented as a network 
graph, where thick lines are used as a shorthand for resistively shunted Josephson junctions, or RSJ (see Fig. 18, 
right panel). A convention for the direction of all branches has to be chosen-in Fig. 18, the direction of branches is 
represented by an arrow. 



2. The tree of the network graph 

As a second step, a tree of the network graph needs to be specified. A tree of a graph is a set of branches connecting 
all nodes without containing any loops. Here, the tree is chosen such that it contains all capacitors, as few iuchictors 
as possible, and neither resistors (external impedances) nor current sources. The conditions under which such a choice 
can be made are discussed in (Burkard et al., 2004b). The tree of Fig. 18(right) that will be used here is shown in 
Fig. 18 (center). The branches in the tree arc called tree branches; all other branches arc called chords. Each chord is 
associated with exactly one so-called fundamental loop that is obtained when adding the chord to the tree. 



3. The loop matrices 



In a next step, the loop sub-matrices Fql, ^cz, Fcbi Fkl, Fkz, and F^b need to be found. The loop sub- 
matrices FxY have entries -|-1, —1, or 0, and hold the information about which tree branches of type X belong to 
which fundamental loop associated with the chords of type Y. E.g., for our example. 



FCL 



V 



1 o\ 

-1 1 

-ij 



(80) 



where the first column determines that the capacitor Ci (part of Ji) belongs to the large loop (associated with Li), 
capacitor C2 (part of J2) belongs to the large loop (with different orientation), while capacitor C3 (part of J3) does 
not belong to the large loop at all. Similarly, the second column of Fcl says which of the capacitors are contained in 
the small loop (associated with L3). 

The loop matrices have the purpose of systematically incorporating Kirchhoff 's laws of current and energy conser- 
vation in the circuit. 



F^'^'^I = (IIP) 1 = 0, 
pWV = (-F^|l)V = *, 



(81) 
(82) 




FIG. 18 Left: The IBM qubit is an example of a network graph. Each thick line represents a Josephson element, i.e. three 

brandies in parallel, sec right panel. Thin lines represent simple two-terminal elements, such as linear inductors (L, K), external 
impedances (Z), and current sources (Is). Center: A Josephson subgraph (thick line) consists of three branches; a Josephson 
junction (cross), a shunt capacitor (C), a shunt resistor (R), and no extra nodes. Right: A tree for the circuit shown on the 
right. A tree is a subgraph containing all nodes and no loop. Here, a tree was chosen that contains all capacitors (C), some 
inductors (K), but no current sources (7s) or external impedances (Z). 
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where the external magnetic fluxes are denoted with and where the loop sub-submatrices 

p _ j Fcj FcL FcR Fez Fob \ ^^^-^ 
\ Fkj Fkl Fkr Fkz Fxb J 

are related to the full fundamental loop and cutset matrices F^^^ and F^*^) via the grouping of the branch currents 
and voltages into a tree and a chord part, I = (Itr j Ich) and V = (Vtr, Vch)- 



4. Current-voltage relations (CVRs) 

In order to derive the equations of motion and eventually, the Hamiltonian of the SC circuit, Kirchhoff's laws, 
Eqs. (81) and (82) need to be combined with the CVRs of the various branch elements. For this purpose, the tree 
and chord currents and voltages are divided up furtlier. according to the various branch types, 1^^ = (1c,1k), and 
Ich = (Ij,Ii,,Iij,Iz,lB), and similarly for the voltages. The tree current and voltage vectors contain a capacitor (C) 
and tree inductor (K) part, whereas the chord current and voltage vectors consist of parts for chord inductors, both 
non-linear (J) and linear (L), shunt resistors (R) and other external impedances (Z), and bias current sources (B). 
The branch charges and fluxes {X = C, K, J, L, R, Z, B) are formally defined as 

lx{t) = Qx{t), (84) 
Yx{t) = ^x{t). (85) 

Using the second Josephson relation and Eq. (85), we identify the formal fiuxes associated with the Josephson junctions 
as the superconducting phase differences if across the junctions. 

If = h 

where $o = h/2e is the superconducting flux quantum. Each branch type has its own current-voltage relation (CVR); 
e.g., the Josephson junction branches follow the first Josephson relation, 

Ij = Ic sin(p, (87) 

with the critical current matrix !(., while the external impedances are described by the integral relation, 

Yz{t) = f Z{t - T)Iz{T)dT = (Z * Iz){t). (88) 

The total inductance matrix 

is used for the CVR of the chord (L) and tree (K) inductances, 

(90) 




where L and Lk are the self inductances of the chord and tree branch inductors, resp., off-diagonal elements describing 
the mutual inductances among chord inductors and tree inductors separately, and Llk is the mutual inductance matrix 
between tree and chord inductors. 



5. The Hamiltonian 

The elements described above are sufficient to determine the Hamiltonian of the dissipation-free system. 



2 



-Hs = 2^cC-^Qc + 2^ j U{^), (91) 
= - V^^cos<^i + i<^^Mo<^+ J^<^^(N*, + SIb), (92) 
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where Qc are the charges conjugate to the fluxes $j = ($o/27r)(p and C is the capacitance matrix of the circuit. 
The matrices Mq, N, and S are obtained from the inductance and loop matrices L( and F (Burkard et al., 2004b). 
The Hamiltonian Eq. (91) is quantized using the commutator relation 



TT^i'Qc-j =ih5ij. (93) 

ZTT J 

The system including dissipation can be described using the Caldeita-Legget model, 



n = Hs + HB+nsB, (94) 

-E 

a 

HsB = m-ipJ2caXa+^U{ip), (96) 



= iyj^+rriacolxl], (95) 



where Ti^ is the quantized Hamiltonian Eq. (91), Hb is the Hamiltonian describing a bath of harmonic oscillators with 
(fictitious) position and momentum operators and pa with [x^jp^] — ihSap, masses ma, and oscillator frequencies 
LUct- Finally, Hsb describes the coupling between the system and bath degrees of freedom, (p and Xa, where Ca is a 
coupling parameter and m are obtained from the inductance and loop matrices and F (Burkard et al., 2004b). 

The quantum dynamics of the entire system (qubit plus bath of oscillators) is described by the Liouville equation 
p(t) — ^i[H, p{t)] = —iCp{t) for the density matrix p. The state of the system alone is described by the reduced 
density matrix psit) = Tr^ p{t). In the Born-Markov approximation, the master equation for psif) can be written in 
the form of the Redfield equations (Redfield, 1957), 

Pnm{t) = -i0JnmPnm{t) - ^ ^mklPkl{t), (97) 

kl 

where Pnm = {nlpslm) are the matrix elements of p in the eigenbasis \n) oiTis (eigcncucrgies a;„), and cJnm = w„— 
and with the Redfield tensor, 

Rnmkl = ^irra'E ^irrfc ^nky^^ ^\rrm ~ ^\mnk ~ ^\mnki (^^) 
r r 

Rer[:i, = (m.y),^(m.y)„,J(Kfe|) jJ^j]^^^^^^ , 

I^rlmLfe = -(m- y)tm(m- <f)nk-pj du> ft^-a;nfcC0th^] . (99) 

In the two-dimensional qubit subspace, the Bloch vector p = Tr(<Tp) can be introduced where cr = {ax,<yy,<Jz) are 
the Pauli matrices, and the Redfield equation (97) takes the form of the Bloch equation p = w x p — i?p + po, with 
ijj = (0, 0, woi)"'", where in the secular approximation, the relaxation matrix R is diagonal, R = diag(T2"^, T2"^, ^). 
The relaxation and decoherence times Ti and T2 are then given by 

^ = 4|(0|m.¥>|l)|V(u;oi)coth^, (100) 



^ = |(0|m.¥,|0)-(l|m.<^|l)|2:^ 

J-6 U! 



2kBT. (102) 

w->0 

In the semiclassical approximotion, Ti and T^j, can be related to the parameters A and e in the Hamiltonian Eq. (79), 



k = [^) I^V->r^("..)coth^, (103) 



, , |A(^-m|- 



2kBT. (104) 
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FIG. 19 Schematic of the Delft circuit, Fig. 16, where the crosses 
denote Josephson junctions. The outer loop with two junctions 
L and R forms a dc SQUID that is used to read out the qubit. 

The state of the qubit is determined by the orientation of the 
circulating current in the small loop, comprising the junctions 
1, 2, and 3, one of which has a slightly smaller critical current 
than the others. A bias current Ib can be applied as indicated 
for read-out. 



6. Leakage 

The leakage rate from a qubit state A; = 0, 1 into higher levels n = 2, 3, . . . outside the qubit space can be quantified 
from the Redfield equation Eq. (97) by the sum 

= 4^ |(n|m • ^\k)\^JM coth (105) 

J-L,k IKbJ- 



7. The Delft qubit 

A very successful qubit design is the Delft qubit (Chiorescu et ai, 2003) which is depicted in Fig. 16, and which 
will be discussed in this Section. A schematical drawing of the SC circuit for the Delft qubit is shown in Fig. 19. 
This design is intended to be immune to current fluctuations in the current bias Ib due to its symmetry properties; 
at zero dc bias, Ig = 0, and independent of the applied magnetic field, a small fluctuating current 5lB{t) caused by 
the finite impedance of the external control circuit (the current source) is divided equally into the two arms of the 
SQUID loop and no net current flows through the three-junction qubit line. Hence, in the ideal circuit (Fig. 19) the 
qubit is protected from decoherence due to current fluctuations in the bias current line. However, asymmetries in the 
SQUID loop may spoil the protection of the qubit from decoherence. In the case of an inductively coupled SQUID 
(Mooij et al, 1999; Orlando et al, 1999; van der Wal et al, 2000) neither a small geometrical asymmetry (imbalance 
of self- and mutual inductances in the SQUID loop), nor thejunction (critical current) asymmetry of typically a few 
percent, would suffice to cause a relevant amount of decoherence at zero bias current /b = (Burkard et al, 2004b). 
What turns out to be important here is that the circuit (Fig. 16) contains another asymmetry, caused by its double 
layer structure, being an artifact of the fabrication method used to produce SC circuits with aluminum/ aluminum 
oxide Josephson junctions, the so-called shadow evaporation technique. Junctions produced with this technique will 
always connect the top layer with the bottom layer, see Fig. 20. 

Thus, while circuits like Fig. 19 can be produced with this technique, strictly speaking, loops will always contain an 
even number of junctions. In order to analyze the implications of the double layer structure for the circuit in Fig. 19, 
the circuit is drawn again in Fig. 21(a) with separate upper and lower layers. Note that each piece of the upper layer 
is connected with the underlying piece of the lower layer via an "unintentional" Josephson junction. 

These extra junctions typically have large areas and therefore large critical ciirrents; thus, their Josephson energy 
can often be neglected. In order to study the lowest-order effect of the double layer structure, one can neglect all 
unintentional junctions in this sense and arrive at the circuit Fig. 21(b). It should be emphasized that the resulting 
circuit is distinct from the 'ideal' circuit Fig. 19 which does not reflect the double-layer structure. In the real circuit. 
Fig. 21(b), the symmetry between the two arms of the dc SQUID is broken, and thus it can be expected that bias 
current fluctuations cause decoherence of the qubit at zero dc bias current, is = 0. 

Starting from the circuit graph of the Delft qubit, the circuit theory can be used to find the Hamiltonian of the 
circuit, which can subsequently be analyzed numerically. The double-well minima </3g and was found for a range of 
bias currents and applied external fiux. The states localized at ip^ and cp-^ are encoding the logical |0) and |1) states of 
the qubit. Two special lines in the plane spanned by the bias currents and applied external fiux can now be determined, 
see Fig. 23. (i) The fine /*(/b) on which a symmetric double well is predicted, e = U{(pf^) — U {(fi) = 0. On this line, 
the dephasing time diverges, (ii) The line on which m • A(p = 0, where Aip = y'oVi the vector joining the two 
minima of the potential. On this line, the environment is decoupled from the system, and both the relaxation and 




FIG. 20 Schematics of Josephson junctions produced by the 
shadow evaporation technique, connecting the upper with the 
lower aluminum layer. Shaded regions represent the aluminum 
oxide. 
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FIG. 21 (a) Double layer structure. Dashed blue lines represent 
the lower, solid red lines the upper SC layer, and crosses indi- 
cate Josephson junctions. The thick crosses are the intended 
junctions, while the thin crosses are the unintended distributed 
junctions due to the double-layer structure, (b) Simplest circuit 
model of the double layer structure. The symmetry between 
the upper and lower arms of the SQUID has been broken by 
the qubit line comprising three junctions. Thick black lines de- 
note pieces of the SC in which the upper and lower layer are 
connected by large area junctions. 
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FIG. 22 External circuit attached to the qubit (Fig. 19) that 
allows the application of a bias current Ib for qubit read-out. 
The inductance Lgh and capacitance Csh form the shell circuit, 
and Z{u!) is the total impedance of the current source (is). The 
case where a voltage source is used to generate a current can be 
reduced to this using Norton's theorem. 



the decoherence times diverge, ri.2,0 00. The curve J*{Ib) agrees qualitatively with the experimentally measured 
symmetry line (Bertet et al., 2004), but it underestimates the magnitude of the variation in flux /' as a function of 
Ib ■ The point where the symmetric and the decoupling lines intersect coincides with the maximum of the symmetric 
line, as can be understood from the following argument. Taking the total derivative with respect to Is of the relation 
e = U{iPq; f*{lB), Ib) — U{ipi; f*{lB),lB) = on the symmetric hue, and using that ipgi are extremal points of U, 
we obtain n • A(p df* /OIb + (27r/$o)m • = for some constant vector n. Therefore, m • A(p = (decoupling line) 
and n • Aip ^ implies dj* /diB = 0. 

The relaxation and decoherence times Ti and T2 on the symmetric line have been evaluated and are plotted (Fig. 23, 
right) where e = 0, and therefore, E = A. The divergence in on the symmetric line is cut off by higher-order 
effects, whereas the divergence of Ti on the decoupling line is cut off by residual impedances, e.g., due to the junctions' 
quasiparticle resistance (Burkard et ai, 2004a). A peak in the relaxation and decoherence times where predicted from 




FIG. 23 Left: Decoupling (red solid) and symmetric (blue dashed) curves in the {Ib,/') plane, where Ib is the applied bias 
current and /' — 2n^'^/^o is the dimensionless externally applied magnetic flux threading the SQUID loop. Both curves are 
obtained from the numerical minimization of the potential Eq. (92). The decoupling line is determined using the condition 
m ■ Aip — 0, whereas the symmetric line follows from the condiction e = 0. Right: Theoretical relaxation, pure dephasing, and 
decoherence times Ti, T^, and T2 as a function of applied bias current Ib, along the symmetric line (Fig. 23, left). 
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FIG. 24 Left: Plasma frequency Usi-i as a function of the applied bias current Is- The variation is due to the change the 
effective in Josephson inductances as Ib is varied. Right: Rabi frequency of the coupling between the qubit and the plasmon 
mode. The coupling disappears at the crossing with the decoupling line (Fig. 23), i.e., when m ■ Aip — 0. 

theory can be observed experimentally (Bertet et ai, 2004). 

The symmetry breaking due to the double layer structure has another, very interesting, consequence. It causes a 
coupling between the qubit and an external harmonic oscillator, the plasmon mode formed by LC resonator in the 
SQUID (Fig. 22). This coupling can be observed as resonances in the microwave spectrum of the system. Moreover, 
it can be used to entangle the qubit with another degree of freedom (Chiorescu et ai, 2004). Here, the inductance Lgh 
and capacitance Cgh of the "shell" circuit (plasmon mode) are responsible for this resonator. This coupling has been 
studied quantitatively in the framework of the circuit theory. From the full Hamiltonian Ti-s, a two-level Hamiltonian 
in the well-known Jaynes-Cummings form can be derived, 

n^Aa, + ea, + huj.h (^b^b + ^ + ^'^z{b + b^), (106) 

with the coupling parameter (Rabi frequency) 

Note that this coupling vanishes along the decoupling line (Fig. 23, left) and also rapidly with the increase of Lsh- 
The Rabi frequency at Ib ~ is predicted to be A « 210 MHz. The Rabi frequency as a function of the bias current 
Ib, together with the plasma frequency, is plotted in Fig. 24. Strong coupling between a charge qubit and a quantum 
electromagnetic cavity formed by a SC transmission line has been observed in (Wallraff et ai, 2004). 

E. Charge qubits 

In analogy to the circuit theory for flux qubits, a general circuit theory for charge qubits will be outlined and 
illustrated with examples in this Section (Burkard, 2004). As in the case of the circuit theory for flux qubits, we are 
not restricted to a Hilbert space of the SC device which is a priori truncated to two levels only. The role of the self 
and mutual inductances in SC charge qubits have been previously studied (You et ai, 2001), in particular as a means 
of coupling two SC charge qubits (Makhlin et ai, 2001; Makhlin and Shnirman, 2004). Here, the self and mutual 
inductances in the underlying SC circuit are fully taken into account. 

1. Graph theory 

Note that the circuit theory for charge qubits is dual to that for flux qubits in the sense that the roles of chord and 
tree branches are interchanged, 

It, = (Ij,lL,Iy,Iz), Ich-(Ic.,Ic,lK), (108) 

Vtr-(Vj,VL,Vy,Vz), Veh-(Vc,,Vc,VK). (109) 




m • Acp. 



(107) 



36 



The loop matrix F now acquires the block form, 



F = 



/ 1 Fjc FjK \ 

Flc Flk 

Fvc FvK 

V Fzc FzK ) 



(110) 



The Hamiltonian has the form 



with the potential 



= 2 (Q - Cv.V)^C-i (Q - CyV) + [/(*), 



(111) 



(112) 



where the Joscphson and inductor flux variables are combined in $ = ((^, $l), with the vector operator of conjugate 
charges denoted by Q. The coupling Hamiltonian in a Caldeira-Leggett description H = Tis + + TisB now takes 
the form 



1-183= C ■^m • Q ^CcCCc = m • C 



(113) 



where C is the total capacitance matric of the circuit. The resulting Redfield equation takes the form Eq. (97) and 
Eq. (98), but with 



Rer 



(+) 

Imnk 



imr(+) 

Imnk 



i(m . Q).„(m . Q)«^-/(l-n'=l)-^^7^^^ , 

1 2 /■°° J(lij) ( 
--(m- Q)i„(m- Q)„fe-P/ rfw-j 5- ('^"'^"fe ^oth 



and with m = C m. Finally, the relaxation and decoherence times in a two-level description reduce to 



1 
1 



^|(0|m.Q|l)|2j(a;oi)coth^, 



1 1 



1 

% 

In the semiclassical limit, one finds 



iKO|m.Q|0)-(l|m.Q|l)|^:^ 



= -|m-AQ| 



Ti h 



A 



J(a;oi) coth 



01 



2A;bT' 



— = -Im-AQp — 4 



2fcBr. 

0)^0 



(114) 

(115) 
(116) 

(117) 

(118) 
(119) 



The leakage rates from the logical state fc = 0, 1 to states n = 2, 3, . . . outside the computational subspace can be 
estimated as 



Tl h 



n=2 



2A:bT' 



(120) 



2. Single charge box 



We now illustrate the circuit theory for charge qubits with some examples. The first example is the voltage-biased 
charge box, shown in Fig. 25. The inductance of the leads has been neglected for simplicity (no L and K branches). 
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FIG. 25 Circuit graph of a single voltagc-biascd charge box. 
Branches represent a Josephson junction (Ej), capacitances (Cj 
and Cg), a voltage source V, and the impedance Z. The nodes 
are shown as black dots; the node connecting the junction {Ej) 
to the gate capacitance Cg represents the SC island. 




FIG. 26 A flux-controlled Josephson junction. 



The tree of the graph is given by the Josephson, voltage source, and impedance branches. For the loop matrices, we 
simply find Fjc = Fyc = Fzc = 1- With the capacitances C = Ctot = Cj + Cg and Cy = Cg, we arrive at the 
Hamiltonian, 



ns = 



{Qj+CgVf 



2Ci 



tot 



+ Ej cos (f. 



(121) 



The coupling to the environment is characterized by m = {Cg/Ctot)- As an example, we give here the relaxation and 
dephasing times, with m = |m| = Cg/Ctot, 



— = 27rm 4|(0|n|l)| — '^oicoth^^, 

1 „ 2,,„, ,,,,24Re.^{0) 2A:bT 

— = 2WK0|n|0)-(l|n|l)p— 



(122) 
(123) 



where n = Q/2e and Rq = h/e^. In the semiclassical limit, (0|n|l) w (l/2)(A/a;oi)An and (0|n|0) — (l|n|l) « 
(e/a;oi)An. With An w 1, we reproduce the results in (Makhhn et al, 2001). Typical leakage rates are of the form of 
1/Ti, with the matrix element replaced by |(0|n|/c)| and |(l|n|A;)|, where k>2 labels a state other than the two qubit 
states, and with wqi replaced by LVik (/ = 0, 1). 



3. Flux-controlled Josephson junction 

A flux-controlled Josephson junction is a SC loop with two junctions which acts as an effective Josephson 
junction with a flux-dependent Josephson energy (Makhlin et al, 1999). The circuit Fig. 26 we use to describe the the 
flux-controlled junction comprises a chord inductance {K) with inductance L. The tree consists of the two Josephson 

branches. The only relevant loop matrix is Fjx = 1 — 1 ^ . In the limit L ^ 0, and if Ej\ = Ej2, we find 
Fj/fV + = — <P2 + ^ ^ 0, which leads us to the Hamiltonian 

ns = ^-Ej{^)cosv, (124) 
where ip = 7r$/$o, C = Cji + Cj2, and Ej{^) = 2Ej cos(27r$/$o)- 
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FIG. 27 Two inductively coupled charge boxes. 



4. Inductively coupled charge boxes 

We now turn to the case of two charge boxes of the type discussed in Sec. IV. E. 2, coupled via an inductive loop 

(Makhlin et al., 1999, 2001), as shown in Fig. 27. Here, the tree consists of all Josephson, voltage source, and 
impedance branches, plus the inductive branch L, and the loop matrices are 




'LC 



(") 



(125) 



With the two capacitance matrices C = diag(Ci,C2) and Cj = diag(Cji, Cj2), we find C 



tot 



C C 



JV 



CjL — Cj^y — (Cl, 6*2)^, and Cj^ = Ci +C2. The vector m consists of the two parts mj = C and 1x1^=^01 C2 ^ • 
With Eq. (Ill) and the inverse of the total capacitance matrix, 



/ {Ci- 


h 0^)0 J2 — C| 
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( ^eff,l 
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^f:S,\2 


^-1 
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C1C2 


(Cl + C2)Cji-C2 


—C2CJI 






^eff,2 
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— ClCj2 


— C2CJ1 


CjlCj2 ) 










/ 
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where 7 = (C\ + C2)Cj\Cj2 — C\Cj2 — C^Cjx^ the Hamiltonian of the coupled system can be written as, 



(126) 



Us 



( {Qj^+C^Vjf , ^ ^ , {Ql + C^Vi+C2V2? , 

1 T^PT- + ^Ji COS <Pi ) + ^ — h 



,=1,2 V 2aff,i 

iQji + CiVi){Qj2 + C2V2) 



+ 



CeS. 



12 



E 

i=l,2 



2Ceff,L 2L 
(Qji + CiVi){QL + CyVx + (72^2) 



(127) 



Ceff, 



Li 



While the last term in Eq. (127) couples each qubit to the LC mode associated with the inductor L, and is thus 
responsible for the inductive coupling of the qubits, the second last term provides a direct capacitivc coupling between 
the qubits. In the limit Ci <C Cjj, we reproduce the results of (Makhlin et al., 2001); however, there are additional 
terms of order d/Cji, in particular the new term oc l/Cof{ ,i2 in the Hamiltonian that capacitively couples the qubits 
directly. Since the coupled system involves at least four levels (more if excited states of the LC coupling circuit or 
higher qubit levels are included), it can no longer be described by a two-level Bloch equation with parameters Ti and 
T2. We can however fix one of the qubits to be in a particular state, say |0), and then look at the "decoherence rates" 
of the other qubit. To lowest order in Ci/Cji, these rates due to the impedance Zi have the form (gj = Ci/{Ci + C2)) 



1 
1 



2^g24|(oo|ni|10)|2^5£|W^„^eoth; 



01 



Rq 

= 27r9||(00|ni|00)-(10|ni|10)| 



2kBT' 
2 4ReZi(0) 2kBT 



Q 



(128) 
(129) 
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F. Multiple sources of decoherence 

In this Section, we show that the total decoherence and relaxation rates of a quantum system in the presence of 
several decoherence sources are not necessarily the sums of the rates due to each of the mechanisms separately, and 
that the corrections to additivity (mixing terms) can have both signs (Burkard and Brito, 2004). To this end, we 
investigate the decoherence due to several sources in superconducting (SC) flux qubits; the general idea of the present 
analysis may however be applied to other systems as well. 



1. Dissipative dynamics 

As an example, the gradiometer qubit drawn in Fig. 28 with n = 3 junctions will be discussed. The gradiometer 
qubit is controlled by applying a magnetic flux $c to the small loop on the left by driving a current Ibi in a coil 
next to it, and simultaneously by applying a magnetic flux $ on one side of the gradiometer using Ib2- The classical 
equations of motion of the SC circuit have the form 



dU 
dip 



(130) 



where C is the capacitance matrix, U the potential, and M(t) the dissipation matrix. The convolution is defined as 
{f*g){t) = f {t — T)g{T)dT . The dissipation matrix in the Fourier representation can be found from circuit theory 
(Burkard et al., 2004b) and has the form 



M.{uj) = TaLz{ui) ^m^ = m(Lz(a;) +Lc 



-1 - T 



(131) 



where m and Lc denote real matrices that can be obtained from the circuit inductances. One can assume the matrices 
Z and Tiz to be diagonal because the impedances Zi are independent. 

A Caldeira-Leggett Hamiltonian Ti = Tis +'Hb +'Hsb can be constructed that reproduces the classical dissipative 
equation of motion, Eq. (130), and that is composed of parts for the system (S), given in Eq. (91), for m > 1 harmonic 
oscillator baths (B), and for the system-bath (SB) coupling. 



T^SB 




(132) 
(133) 



where = {Xai, ■ ■ ■ ,Xam), and Cq is a real nx m matrix. Defining the matrix spectral density J{u)) of the 
environment, where Sij{X.) = 5(Xy), one obtains the relation 



JM = |X]^«™«^'^«^^('^~'^«)^«'= (^) ImM(a;) = 5]l-^j('^)mj(a;)mj(a; 



(134) 



where the spectral decomposition (with the eigenvalues Jj{oj) > and the real and normalized eigenvectors mj{co)) of 
the real, positive, and symmetric matrix ImM(w) has been used. The integer m < n,nz denotes the maximal rank of 









h 
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FIG. 28 Circuit graph of the gradiometer qubit (Koch et al., 
2004), under the influence of noise from two sources Zi and Z2. 
Branches of the graph denote Josephson junctions Ji, induc- 
tances Li and Ki , current sources Isi , and external impedances 

Zi, and are connected by the nodes (black dots) of the graph. 
Inset: A resistivoly-sliunted Josephson junction (RSJ) Ji, rep- 
resented by a thick line in the circuit graph, is modeled by an 
ideal junction (cross) with critical current Id, shunt resistance 
Ri, and junction capacitance Cj. 
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ImM(ti)), i.e., m = max^ (rank [IniM(a))]). Using Eq. (134), and choosing Cdj = jaji^ii^aj), we find that Jj{uj) is the 
spectral density of the j-th bath of harmonic oscillators in the environment, Jj(oj) = (7r/2) X)Q,(7^j/'TiajWaj)<5(w— 
From this model one obtains the Redfield equation, Eqs. (97), with a Redfield tensor, Eq. (98), of the form 

R.r«, = ^r„J(l--l)»'„.-|^^, (135) 
■-r|;L, = -ipr^f£^?<^(.-.„.co.hf), 

where <p j. = (n|<^|A;). For two levels n = 0, 1, and within the secular approximation, the relaxation and decoherence 
rates Tj~ and T2 are found to be 

Tf ^ = 4<plj{woi)fo, coth (^) = \<foi-^jM\'JjM coth (^) , (136) 



y 2 J ~ A^^^oi "'jv^'uiyi ^jv^ui; ^ 2 J' 
= |lim(<^oo-V'ii)^^(¥'oo-¥'ii) = |E|mi(0)-(¥'oo-<^ii)l'^ 



(137) 



where we have used the spectral decomposition, Eq. (134), to obtain the second equality in each line. 



2. Mixing Terms 

In the case where is diagonal, or if its ofF-diagonal elements can be neglected because they axe much smaller 
than Lz{uj) for all frequencies u>, one finds, using Eq. (131), that the contributions due to different impedances Zi are 
independent, thus m = nz and M(w) = mL2(w)~^m^ = mjYhJiuj/{Zj{u}) + iivLjj), where m^- = uij is simply 
the j-th column of the matrix m and Ljj is the j-th diagonal entry of Lc. As a consequence, the total rates 1/Ti and 
l/T^ are the sums of the individual rates, 1/t[^^ and 1/t'^\ where 

1 /i/^^oVi - 121} woiCoth(/3u;oi/2) 

In general, the situation is more complicated because current fluctuations due to different impedances are mixed by 
the presence of the circuit. In the regime Lc <C \^z{^)^ can be expanded as 

L-/ = (L2(u;) + Lc)-^ = L^^ - L^'LeL^i + L^^LeL^^LcL^' - • • • . (140) 
The series Eq. (140) can be partially resummed, 



where the first term simply gives rise to the sum of the individual rates, as in Eqs. (138) and (139), while the second 
term gives rise to mixed terms in the total rates. The rates can therefore be decomposed as {X = 1,2, cf)) 

One finds for the mixing term in the relaxation rate 

4^^) 'v'JimImL-L(a;oi)m^Voi coth(^) . (143) 



rTi(mix) 

One can show that there is no mixing term in the pure dcphasing rate, i.e., i/j'^™"^) = Q, and hence y^'™'^' = 2^^™"'-'. 
In the case of two external impedances, nz = 2, Eq. (140) can be completely resummed, with the result 
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where are the matrix eleraents of Lc and where the approximation in Eq. (144) holds up to 0(Z ■^). In lowest 
order in l/^i, one finds, 

1 ^-f^Vlm /;^--|t"^ cothf^V (145) 




y^(mix) \2t:) Zi(cjoi)^2(c^oi) V 2 

with Lpi2 = (tpoi • mi)(<^oi • m.2). 

If = Zi(tjjQi) are real (pure resistances) then, as predicted above, the imaginary part of the second-order term 
in Eq. (144) vanishes, and we resort to third order, 

ImL-i = ( ^ ^ -C ) (146) 

neglecting terms in 0{R~'^). If L12 ^ Ljj, we obtain ImL^^?^ « {uj'^ L 12/ RiR2){Lii/ Ri + L22/ R2)(^x^ and 

For the gradiometer qubit (Fig. 28), we find L12 ~ Afi2Mi3M34/LiL3, Ln k, L2, L22 ~ i^4, where Lk denotes the 
self- inductance of branch Xk {X—L or K) and Mki is the mutual inductance between branches Xk and Xi, and where 
we assume My ^ L^. The ratio between the mixing the single- impedance contribution scales as 



1/n 



0) i?2 
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where we have assumed i?i ~ i?2 = R, Ln w L22 = L, and y>Q]^ • riii w (^qj^ • m2- 

The relaxation time Ti was calculated at a temperature T <C fe^oi/^s for the circuit Fig. 28, for a critical current 
Ic — 0.3 /lA for all junctions, and for the inductances Li = 30 pH, L3 — 680 pH, L2 = L4 = 12 nH, M12 ~ -y/LiX^, 
M34 ~ ^/L^LI (strong inductive coupling), M35 = 6pII, with wqi = Stt • 30 GHz, and with the impedances Zi = i?, 
Z2 = i? -f iR\ia, where R and i?im = ±10kf2 are real (i?im > corresponds to an inductive, i?im < to a capacitive 
character of Z,i). In Fig. 29, Ti was plotted with and without mixing for a fixed value of M13 = 0.5 pH and a range 

of i? = ReZ;. In the inset of Fig. 29, Ti (with mixing) and {{t[^'^)~^ + (rf V^)"^ (without mixing) are plotted for 
i? = 75 SI for a range of mutual inductances M13; for this plot, the double minima of the potential U and ip^^i were 
computed numerically for each value of M13. The plots (Fig. 29) clearly show that summing the dccohcrcncc rates 
without taking into account mixing term can both underestimate or overestimate the relaxation rate 1 /Ti , leading to 
either an over- or underestimate of the relaxation and decoherence times Ti and T2. 





FIG. 29 The relaxation rate Ti without the mixing term 
(dashed blue line), and including the mixing term for 
i?im = -1-10 kn (solid red line) and 7?im = -lOkSl (dot- 
dashed light blue line), for M13 = 0.5 pH as a function 
of KeZi. Inset: Ti for R = ReZi — 75fl for a range of 
mutual inductances Afis. 
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V. ENTANGLEMENT 

A composite system is entangled if its wavefunction \1/ cannot be expressed as a tensor product "^a (8i ^'b of 
wavefunctions "if a and "ifs for the parts A and B of the system (a more general, but similar, definition exists for 
mixed states). A variety of quantum communication scenarios (Bennett and DiVincenzo, 2000), some of which have 
been implemented successfully in quantum optics, require maximally entangled states of two qubits, also known as 
EPR pairs (Einstein et al., 1935), such as the spin singlet, 

1^) = ^ (in) -lit))- (149) 

The triplet state |To) = | Ti) + | it) is another maximally entangled state, while the other two triplet states |T+) = 1 1|) 
and |T_) = | J,J,) are not entangled. An important feature of these states is that they are non-local, in the sense 
that they violate Boll's inequalities (Bell, 1966; Mcrmin, 1993). A universal quantum computer can, by definition, 
produce arbitrary quantum states, and, in particular, entangled ones such as the singlet \S). E.g., the square root 
of swap gate S, see Eq. (8), has the ability to turn the unentangled state | ti) into a maximally entangled one, 
^\ Ti) = (I Ti) ^ *l iT))/(l ^ *); which is equivalent (up to a singlc-qubit operation) to the singlet. There may be cases 
in which only certain entangled states arc required for quantum communication, while quantum computation itself 
does not need to be performed. In this case, a physical device dedicated to the task of producing entangled states of 
some sort may be sufficient. We call devices of this sort entanglers (Burkard et al., 2000b) and discuss a number of 
conceivable implementations of entanglers in solid state systems below. 

Experiments with entangled photons have tested Bell's inequalities (Aspect et ai, 1982), and various quantum 
communication protocols, such as dense coding (Mattle et al., 1996) and quantum teleportation (Boschi et al., 1998; 
Bouwmeester et al., 1997). However, none of these protocols been implemented so far with massive particles (such as 
electrons). 

Unfactorizable states like Eq. (149) are very common in solid-state systems. Interacting many-particle systems 
possess very complicated and entangled ground states. Not all of these are necessarily useful for quantum information 
processing, though, because (i) it is essential that there is a physical mechanism to extract and separate a pair of 
entangled particles from the many-body system in such a way that they can be used for quantum communication, 
and (ii) for indistinguishable particles, not all states that "look entangled" really are. A measure of entanglement 
which excludes pure antisymmetrization was defined in (Schliemann et al., 2001a, b). 



A. Production of entangled electrons 

1. Superconductor-normal junctions 

A superconductor (SC) with s-wave pairing symmetry contains an entire "reservoir" of spin singlet states as in 
Eq. (149) in the form of Cooper pairs that form the SC condensate (Schrieffer, 1964). It is thus natural to think 
that such a system can act as a source of spin-entangled electrons. A proposed setup (Recher et al., 2001) is shown 
in Fig. 30. The SC is held at the chemical potential /ig, and is weakly coupled by tunnel barriers to two separate 
quantum dots Di and £>2 which are in turn weakly coupled to Fermi liquid leads Li and L2, both held at the same 
chemical potential /iti = /Lt2- The tunneling amplitudes between SC and dots, and dots and leads, are denoted by Tsd 
and Tdl- 

Andreev tunneling is a process in which two electrons (one with spin up and one with spin down) can tunnel 
coherently through a normal barrier, while at the same time, single-particle tunneling as suppressed (Hekking et al., 
1993). The setup Fig. 30 with the intermediate quantum dots is designed to force two electrons from a Cooper pair 
to tunnel coherently into separate leads rather than both into the same lead. The double occupation of a quantum 
dot is suppressed by the Coulomb blockade mechanism (Kouwenhoveu ct al, 1997c). 

The flow of entangled electrons from the SC via the dots to the leads is controlled by applying a bias voltage 
= ^g — fii > 0. The chemical potentials ei and €2 of the two quantum dots can be tuned by external gate voltages 
(Kouwenhoven et al., 1997c) such that the coherent tunneling of two electrons into different loads is at resonance, 
while coherent tunneling of two electrons into the same lead is suppressed by the on-site Coulomb repulsion C/ of a 
quantum dot. 

The relevant parameters describing the device and the desired regime of operation are discussed in (Recher et al, 
2001). It is required that the barriers of the dots are asymmetric, \Tsd\ \T£)l\, temperature is sufficiently small, 
An > ksT, and A, U, 6e > A/j, > ji, ksT, and 7; > 7s, where Se is the single-level spacing of the dots, A is the SC 
energy gap, and 7; = 2'jti'i\T£)l\'^ are the dot-lead tunnel rates. The figure of merit for the device is the ratio between 
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the desired current Jiof pairwise entangled electrons tunneling into different leads and the unwanted current I2 of 
electron pairs that end up in the same lead (Recher et ai, 2001), 



h 
h 



— 



sin(A;i?(5r) 



kpSr 



1 _ 1 1 
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where kp denotes the Fermi wavevector, 7 = 71 + 72, and ^ the SC coherence length. 

The desired current I\ decreases exponentially with increasing distance = |ri — between the tunneling points 
on the SC, the scale given by the superconducting coherence length ^. With ^ typically being on the order of /xm, 
this does not pose severe restrictions for a conventional s-wave SC. In the important case < 6r ^ ^ the suppression 
is only polynomial oc l/ikpSr)"^, with kp being the Fermi wavevector in the SC. One also observes that the effect of 
the quantum dots consists in the suppression factor (j/E)'^ for tunneling into the same lead. One therefore has to 
impose the additional condition kpSr < £/^, which can be satisfied for small dots with £/j « 100 and kp^ pa 1 A. 

As an alternative to quantum dots as a means to separate the entangled electrons from the SC, it hsa been proposed 
to use a Luttinger liquid (see Sec. V.A.2 below) or a resistive lead where the dynamical Coulomb blockade effect helps 
to separate the electron pair (Recher and Loss, 2003). 



2. Superconductor-Luttinger liquid junctions 



In the Andreev entangler (Sec. V.A.I), entangled electron pairs are separated by the Coulomb repulsion in quantum 
dots that are attached to the SC which acts as a source of entangled spin singlets. In related work (Recher and Loss, 
2002a,b) and (Bena et at, 2002), it was suggested that the strong Coulomb interactions in a one-dimensional conductor, 
forming a Luttinger liquid (Tsvclik, 2003) can play the same role. There is good experimental evidence for Luttinger 
liquid (LL) behavior in carbon nanotubes (Bockrath et ai, 1999). 

The setting discussed in (Recher and Loss, 2002a,b) consists of a conventional s-wave SC tunnel- coupled to the 
center (bulk) of two spatially separated, for all practical purposes infinitely extended, one- dimensional wires (e.g., 
carbon nanotubes) each forming a separate LL. While the Coulomb interaction within each wire is essential for the 
separation of entangled pairs into distinct wires, it is assumed that the interaction between carriers in different wires is 
negligible. In the absence of backscattering, the low energy excitations of the LL are long-wavelength charge and spin 
density oscillations propagating with velocities Up = vp/Kp for the charge and = vp for the spin (Schulz, 1990), 
where vp is the Fermi velocity and Kp < 1 due to interaction. Transfer of electrons from the SC to the LL-leads is 
described by a tunneling Hamiltonian, 



HT = toY^ ^L*«(rn)+H.c., 
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where ^'s(r„) annihilates an electron with spin s at the point r„ on the SC nearest to the LL-lead n = 1,2, and 
creates it again with same spin and amplitude to at the point a;„ in LL n. By applying a bias /i = /ig — jii between 
the SC, with chemical potential /xg, and the leads, held at the same chemical potential m, a stationary current of 
pairwise spin-entangled electrons can flow from the SC to the leads. 

As in the case of the Andreev entangler with attached quantum dots, the performance of this device can be quantified 
by the ratio between the two competing currents Ji and I2 (see Sec. V.A.I). From a T-matrix calculation (Recher 
and Loss, 2002a, b), one obtains in leading order in /x/A and at zero temperature, 

27p 
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FIG. 30 The Andreev entangler. A Cooper pair is split up into 
two entangled electrons which hop with amplitude Tsd from two 
points ri, V2 of the superconductor, SC, (distance 5r = |ri — r2|) 
onto two dots Di,2 by means of Andreev tunneling. The dots 
are coupled to normal leads Li,2 with tunneling amplitude Tdl- 
In order to maximize the ofRciency of the device, wo require 
asymmetric barriers, |Tsi3|/|rDi,| <C 1. The chemical potentials 
of the SC and leads are //; and us . 
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where F is the Gamma function, A a short distance cut-off on the order of the lattice spacing in the LL, 7 = 
27ri/5fi|toP the probabihty per spin to tunnel from the SC to the LL-leads, s^nd vi the energy DOS per spin for 
the superconductor and the LL-leads at the chemical potentials and ni, resp., and 6r the separation between the 
tunneling points on the SC. The current Ii has a characteristic non-linear dependence on the voltage (electro-chemical 
potential /z), 7i oc with an interaction dependent exponent 7p = {Kp+Kp^)/4:~l/2 > 0, which is the exponent 

for tunneling into the bulk of a single LL, i.e. p{e) ~ lep'', where p{e) is the single-particle density of states (Schulz, 
1990). In the non-interacting limit 7^ = the current is given by Ii- As in Sec. V.A.I, the coherence length ^ of 
the Cooper pairs should exceed Sr in order to obtain a finite measurable current. Note that the suppression of the 
current by l/{kF5r)'^ can be considerably reduced with the use of lower dimensional SCs (Recher and Loss, 2002a, b). 
The desired current Ii now has to be compared with unwanted current consisting of electron pairs tunneling into the 
same lead and having Sr = 0. It is found (Recher and Loss, 2002a,b) that the current I2 for tunneling into the same 
lead (1 or 2) is suppressed if < A with the result, again in leading order in /i/A, 




(153) 



where Af, is an interaction dependent constant of order one, and where 7p+ = jp, and 7p_ = 7^+ + (1 — Kp) /2 > 7^+. 
Note that in Eq. (153) the current Ii needs to be evaluated at Sr = 0. In the non-interacting limit, I2 = h = Ii is 
obtained by putting 7^ = 7^^ = 0, and Up = vp. The expression Eq. (153) for I2 shows that the unwanted injection 
of two electrons into the same lead is suppressed compared to 7i by a factor of {2n/A)^'^p+ , where 7p+ = 7p, if both 
electrons are injected into the same direction (left or right movers), or by (2/i/A)^'''p" if the two electrons travel in 
different directions. It is more likely that the two electrons move in the same direction than in opposite directions, 
because 7p_ > 7p+. The suppression of the current I2 by 1/A is a manifestation the two-particle correlations in 
the LL when the electrons tunnel into the same lead, which is similar to the Coulomb blockade effect in the case of 
tunneling into quantum dots in Sec. V.A.I. As the SC gap A becomes larger, the delay time between the arrivals 
of the two partner electrons of a Cooper pair becomes shorter, and the effect the first electron in the LL has on the 
second electron tunneling into the LL increases. Increasing the bias /i opens more available channels into which the 
electron can tunnel, and therefore the effect of the SC gap A is loss pronounced. This correlation effect disappears 
when interactions in the LL are absent, 7p = "fph = 0. Experimental systems with LL behavior are e.g. metallic carbon 
nanotubes with similar exponents as derived here (Egger and Gogolin, 1997; Kane et al., 1997). 

3. Transport through quantum dots 

Entanglers with a single quantum dot attached to leads with a very narrow bandwidth (Oliver et al., 2002) or 
with three coupled quantum dots (Saraga and Loss, 2003) have been proposed. The idea behind these proposals is 
the harness the singlet ground state of a single two-electron quantum dot by extracting the two electrons into two 
separate leads. In both proposals, the separation is enhanced due to two-particle energy conservation. A double-dot 
turnstile device with time-dependent barriers was proposed in (Hu and Das Sarma, 2004). 

4. Coulomb scattering in a 2D electron system 

Scanning probe techniques can be applied to a two-dimensional (2D) electron system formed in a semiconductor 
heterostructure in order to monitor and control the flow of electrons (Topinka et al., 2000, 2001). It has been proposed 
to generate spin-entangled pairs of electrons using this technique to control Coulomb scattering in a interacting 2D 
electron system (Saraga et al., 2004). At a scattering angle of 7r/2, it is expected that constructive two-particle 
interference leads to a enhancement of the spin-singlet scattering probability, while the triplet scattering is suppressed. 
The scattering amplitudes have been calculated within the Bethe-Salpeter equation for small rg and allow an estimate 
of the achievable current of spin-entangled electrons (Saraga et al., 2004). 

5. Entangled Electrons in a Fermi Sea 

A particularly appealing aspect of the electron spin as a carrier of quantum information is that it is attached to a 
charge, and thus it can in principle be transported in a conductor. One can therefore envision solid-state structures 
(e.g., on a microchip) where entanglement is produced in one location by one of the previously discussed methods, and 
subsequently conveyed through a wire to the location where entanglement is "used up" in some quantum information 
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protocol. While transporting qubits is quite unproblcmatic in the case of photons as the quantum information carriers 
(photons have been used in many experiments to carry quantum information, even over distances of kilometers), it is 
less trivial for electrons. When an electron is injected into a metallic wire, it is immersed into a sea of other electrons 
that (i) are indistinguishable from the injected electron, and (ii) constantly interact with all other electrons (including 
the injected one) via the Coulomb interaction. In this Section, the stability of spin entanglement in the Fermi sea 
will be discussed (Burkard et ai, 2000b; DiVincenzo and Loss, 1999). It turns out that (i) the indistinguishability of 
particles in the Fermi sea is actually not a problem for the transport of spin qubits and (ii) the Coulomb interaction 
does have some effect, which is however mitigated by the phenomenon of screening which is well-known in interacting 
Fermi liquids (Mahan, 1993). More precisely, when an electron in the orbital state q is added to a Fermi sea (lead), 
the quasiparticle weight of that state will be renormalized by < < 1, i.e. some weight 1 — to find the electron 
in the original state q will be distributed among all the other electrons. Such a rearrangement of the Fermi system 
due to the Coulomb interaction happens very quickly, on a time scale given by the inverse plasma frequency. 

In order to analyze this effect quantitatively, the entangled two-electron state injected into two distinct leads 1 and 
2 can be written in second quantized notation, 

= {al^al^ ± al^ al^ ) |^o), (154) 

where s and t stand for the singlet and triplet, jV'o) for the filled Fermi sea, and n = (qj), where q denotes the 
momentum and I the lead quantum number of an electron. As usual, the operator aj,^ creates an electron in state n 
with spin a. After their injection, the two electrons of interest are no longer distinguishable from the electrons of the 
leads, and consequently the two electrons taken out of the leads will, in general, not be the same as the ones injected. 

The time evolution of the triplet or singlet states, interacting with all other electrons in the Fermi sea, is described 
by the 2-particle Green's function G*/*(12, 34; t) = i|V'34*)) which is related to the standard 2-particle Green's 
function G(12, 34; t) by 

G*/^ (12, 34; i) = - i ^[G(12, 34; t) ± G(12, 34; t)] , (155) 

where n = (n, cr) and n = (n, — cr). Assuming that at time t = 0, a, triplet (singlet) is prepared, then the fidelity of 
transmission 

P(t) = |G*/«(12,12;t)|2 (156) 

is defined as the probability for finding a triplet (singlet) after time t. 

With Eq. (155), the problem reduces to that of evaluating the standard two-particle Green's function G(12, 34; t) = 
— {Tai{t)a2{t)alal} for a time- and spin-independent Hamiltonian, H = Hq + Y^i^j ^iji where Hq describes the free 
motion of the N electrons, and Vij is the bare Coulomb interaction between electrons i and j, (...) denotes the zero- 
temperature expectation value, and T is the time ordering operator. One can assume that the leads are sufficiently 
separated, so that the mutual Coulomb interaction can be neglected, and thus the problem of finding an explicit value 
for G(12,34;t) is simplified since the 2-particle vertex part vanishes (i.e. the Hartree-Fock approximation is exact), 

G(12, 34; t) = G(13, t)G(24, t) - G(14, i)G(23, i), (157) 

with the interacting single-particle Green's functions 

G(n,t) = -i(V'o|Ta„(t)al,|Vo) = Gi(q,t), (158) 

for each lead I = 1,2. Substituting this back into Eq. (155), one obtains 

G*/«(12,34;t) = -G(l,i)G(2,i)[^i3<524T<5i4<523], (159) 

where the upper (lower) sign refers to the spin triplet (singlet). For the initial state [t — 0), or in the absence of 
interactions, one finds G(n, t) = — i, and thus G*/* reduces to (5i3(524 T ^i4<523, and P = 1. The evaluation of the 
(time-ordered) single-particle Green's functions Gi^2 close to the Fermi surface yields the standard result (Mahan, 
1993) Gi^2(q, i) ~ —izq6{eq — ep)e~"-''*-~^''^ , which holds for < f < l/Fg, where l/Fg is the quasiparticle hfetime, 
Cq = q^ /2m the quasiparticle energy (of the added electron), and ep the Fermi energy. For a two-dimensional electron 
system (2DES), as e.g. in GaAs heterostructures, Fg oc (e^ — ep)"^ log(eg — cf) (Giuliani and Quinn, 1982) within the 
random phase approximation (RPA) , which accounts for screening and which is obtained by summing all polarization 
diagrams (Mahan, 1993). The quantity of interest here is the quasiparticle weight, zp = Ti^^iikp, lo))~^ \ ^q, 



46 



evaluated at the Fermi surface, where E].ct('Z,'^) is the retarded irreducible self-energy. For momenta q close to the 
Fermi surface and for identical leads (Gi = G^) we find |G*/''(12, 34; f)p = Zp \ (5i3(524 T '5i4'523p, for times satisfying 
Q <t< l/Tq. The fidelity therefore turns out to be 

P = 4 (160) 



The irreducible self-energy Ej-et and from it the quasiparticle weight factor in two dimensions were evaluated 
explicitly (Burkard et al., 2000b), with the final result 

^"= l + r. (1/2 + 1/.) ' ^'''^ 

in leading order of the interaction parameter Vs = l/kpciB, where Ub = e^h^ /me^ is the Bohr radius. In particular, 
in a GaAs 2DES we have as ^ 10.3 nm, and = 0.614, and thus we obtain zp = 0.665. The expansion in powers 
of Ts for the exact RPA self-energy can be summed up and evaluated numerically, with the (more accurate) result 
Zf = 0.691155 for GaAs. The fidelity of transmission of the injected singlet in this case is around P « 0.2. However, 
for large electron density (small Vg), P is closer to unity. Note the fidelity of the ("postsclcctcd") singlet pairs which 
can successfully be removed from the Fermi sea, is equal to 1, provided that (as assumed here) the spin-scattering 
effects are negligible. That this is indeed the case in GaAs 2DEGs is supported by experiments where the electron 
spin has been transported phase-coherently over distances of up to 100 /xm (Awschalom and Kikkawa, 1999; Kikkawa 
and Awschalom, 1998; Kikkawa et al, 1997). 



B. Detection of spin entanglement 

Efforts arc being made to produce spin entanglement in solid-state structures; therefore, it is only natural to ask 
how one can test for the presence of entanglement in such a setting. Here, a variety of tests for spin entanglement are 
described. This investigation touches on fundamental issues such as the non-locality of quantum mechanics, especially 
for massive particles, and genuine two-particle Aharonov-Bohm effects which are fascinating topics in their own right. 
The main idea in all of the following detection schemes is to exploit the unique relation between the symmetry of the 
orbital state and the two-electron spin state which makes it possible to detect an electron spin state via the orbital 
(charge) degrees of freedom. 



1. Coupled quantum dots 

The first setup to be considered can be used to probe the entanglement of two electrons localized in a double-dot by 
measuring a transport current and its fiuctuations, or current noise (Loss and Sukhorukov, 2000). It is assumed that 
the double-dot is weakly coupled to in- and outgoing leads (at chemical potentials 1^1,2) with tunneling amplitude 
T, where the dots are shunted in parallel. The regime of interest is (i) the Coulomb blockade regime (Kouwenhoven 
et at, 1997c) where the charge on the dots is quantized and (ii) the cotunneling regime (Averin and Nazarov, 1992; 
Konig et al, 1997), where single-electron tunneling is forbidden by energy conservation. The latter regime is defined 
hy U > l/iti ± /i2| > J > ksT, where U is the single-dot charging energy, v the lead density of states, and J 

the exchange coupling between the dots. The current in the cotunneling regime is generated by a coherent virtual 
process where one electron tunnels from a dot to, say, lead 2 and then a second electron tunnels from lead 1 to this 
dot. If the bias voltage is larger than the exchange coupling, — /i2| > J , elastic as well as inelastic cotunneling 
will occur. It will be assumed that T is small enough for the double-dot to return to its equilibrium state after each 
tunneling event. An electron can either pass through the upper or lower dot, therefore a closed loop is formed by 
these two paths. A magnetic flux then gives rise to an Aharonov-Bohm phase (p = ABe/H {A being the loop area) 
between the upper and the lower paths leading to quantum interference efi'ects. This transport setting is sensitive 
to the spin symmetry of the two-electron state on the double dot; if the two electrons on the double-dot are in the 
singlet state, then the tunneling current acquires an additional phase of tt leading to a sign reversal of the coherent 
contribution compared to that for triplets. In cotunneling current, this additional phase manifests itself in the sign 
of an interference term (Loss and Sukhorukov, 2000) 

/ = eiTiy^T^ Mi^M2 ^ ^ ^^g2) 

where the upper sign refers to the triplet states in the double-dot and the lower sign to the singlet state. The shot 
noise is Poissonian with power S{0) = — e|/|, hence it has the same dependence on the state on the double dot. 
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The shot noise has also been calculated for finite frequencies in (Loss and Sukhorukov, 2000), and it was found 
that again S{uj) cx (2 ± cosc/)), and that the odd part of S{ui) leads to slowly decaying oscillations of the noise in real 
time, S{t) oc sm{fj,t)/ fj,t, fi = {fii + /U2)/2, which can be ascribed to a charge imbalance on the double dot during an 
uncertainty time /U~^. 

Note that while the scheme described above is able to distinguish the states \S) and |To) on the dots, the three 
triplets |To), |T+), and |T_), can be further distinguished by an orientationally inhomogeneous magnetic field which 
results in a spin-Berry phase (Loss and Goldbart, 1992; Loss and Sukhorukov, 2000) that leads to left, right or no 
phase-shift in the Aharonov-Bohm oscillations of the current (noise). 



2. Coupled dots with SC leads 

A related scenario of double-dots (DD) has been considered in (Choi et ai, 2000), where two quantum dots are again 
shunted in parallel between the leads, but without any direct coupling between them. The two dots are assumed to be 
coupled via tunneling (with amplitude T) to two superconducting leads. It turns out that the s-wave SC energetically 
favors an entangled singlet-state on the dots. In addition to this, the coupling to the SC provides a mechanism for 
detecting the spin state via the Josephson current through the double dot system. In leading order oc T*, the spin 
coupling is described by a Heisenberg Hamiltonian (Choi et ai, 2000) 

H,s « J (1 + cos if) (s^ • S6 - J) , (163) 

where J « 2T^/e, e is the energy difference from the dot state to the Fermi level, and ip is the average phase difference 
across the SC-DD-SC junction. The exchange coupling between the spins can be controlled by tuning the external 
parameters T and (p, thus providing another implementation of a two-qubit quantum gate (Sec. I. A. 4) or cntanglcr 
(Sec. V.A). The spin state (singlet or triplet) on the dot can be probed if the SC leads are joined with one additional 
(ordinary) Josephson junction with coupling J' and phase difference 9 into a SQUID. The SC current Is flowing in 
this ring is given by (Choi et ai, 2000) 

Ig/Ij = I ^^^^^ ~ '^'^^^ + i^'l-J) sin^ > singlet, ^^^^^ 
1 {J' / J)sm9 , triplets, 

where Ij = 2e J/ h. The spin state of the DD is now probed by measuring the spin- and flux-dependent critical current 

Ic = maxe{|/s|} by biasing the system with a dc current I until a finite voltage V appears for |J| > Ic (the SC goes 
into the normal state). Another interesting effect is long-distance coupling between spins residing in dots separated 
by 6r which is induced by the presence of the SC. The resulting exchange coupling is 



J{6r) = J(0) 



sm{kpdr)'' ^ 



kpSr 



e-2'5'-/«. (165) 



3. Beam splitter shot noise 

Pairwise spin entanglement between electrons in two mesoscopic wires can be detected from a charge current 
measurement after transmission through an electronic beam splitter (Burkard et ai, 2000b). In this scheme, the 
singlet EPR pair Eq. (149) gives rise to an enhancement of the shot noise power ("bimching" behavior), whereas 
the triplet EPR pair \To) leads to a suppression of noise ("antibunching"). This behavior can be anticipated from 
a textbook example: the scattering theory of two identical particles in vacuum (Ballentine, 1990; Feynman et al., 
1965). There, the differential scattering cross-section a in the center-of-mass frame can be expressed in terms of the 
scattering amplitude f{6) and scattering angle 6 as 

a{e) = \f{e) ± /(tt - e)\' = \m\' + \f{n - e)f ± 2Rer{e)f{n - e), (166) 

the first two terms on the right being the contributions which would be obtained if the particles were distinguishable, 

and the third (exchange) term the contribution due to the particles' indistinguishability. This last term gives rise to 
genuine two-particle interference effects. Here, the plus (minus) sign applies to spin-1/2 particles in the singlet (triplet) 
state, described by a (anti)symmetric orbital wave function. The very same two-particle interference mechanism which 
is responsible for the enhancement (reduction) of the scattering cross section a{d) near 9 = tt/2 also leads to an increase 
(decrease) of the correlations of the particle number in the output arms of a beam splitter (Loudon, 1998). 
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Beam splitter 



FIG. 31 Setup for measuring the noise of entangled states. Uncorrelatcd electrons are fed into the entangler (see text) through 
the Fermi leads 1', 2' and are transformed into pairs of electrons in the entangled singlet (triplet) state |=f), which are injected 
into leads 1, 2 (one electron of undetermined spin state into each lead). The entanglement of the, say, spin singlet can then be 
detected in an interference experiment using a beam splitter (with no backscattering) : Since the orbital wave function of the 
singlet is symmetric, the electrons leave the scattering region preferably in the same lead (3 or 4). This correlation ("bunching") 
is revealed by an enhancement of the noise by a factor of 2 in the outgoing leads. 



For the detection of spin entanglement of electrons carried by two mesoscopic wires, we propose a non-equilibrium 
transport measurement based on the set-up shown in Fig. 31. The beam splitter ensures that the electrons leaving the 

entangler (see Sec. V.A) have an amplitude t to be interchanged (without mutual interaction) such that < \t\'^ < 1. 
In the absence of spin scattering the noise measured in the outgoing leads 3 and 4 exhibits bunching behavior for pairs 
of electrons with a symmetric orbital wave function (Hanbury Brown and Twiss, 1956), i.e., for spin singlets, while 
spin triplets induce antibunching behavior, due to their antisymmetric orbital wave function. The latter situation 
has been considered for electrons in the normal state both in theory (Biittiker, 1990, 1992; Martin and Landauer, 
1992) and in recent experiments (Henny et al, 1999; Liu et al., 1998; Oliver et al., 1999; Torres and Martin, 1999). 
The experiments have been performed in semiconductor nanostructures with geometries that arc closely related to 
the set-up proposed in Fig. 31 but without entangler. It should be stressed here that if bunching (enhancement of 
shot noise) is detected, it can be interpreted as a unique signature for entanglement of the injected electrons, since 
the maximally entangled singlet is the only state leading to this effect. The effect of interactions in the leads have 
already been assessed in Sec. V.A. 5. In order to determine the shot noise of spin-entangled electrons, the standard 
scattering matrix approach (Biittiker, 1990, 1992), extended to a situation with entanglement, is appHed. 

We start by writing down the entangled incident state |±) = |V'i2*)) where we set n = (e„,n), now using the 
electron energies e„ instead of the momentum as the orbital quantum number in Eq. (154) and where the operator 
a]j^(£) creates an incoming electron in lead a with spin a and energy e. The theory for the current correlations in 
a multiterminal conductor (Biittiker, 1990, 1992) can easily be generalized to the case of entangled scattering states, 
with the important consequence that Wick's theorem does not apply. The operator for the current carried by electrons 
in lead a of a multiterminal conductor can be written as 

^«(*) = 7^ E [4a(£)«aa(£') - bi,{e)b^„is')] exp [i{s - s')t/h], (167) 

ee'a 

where the operators baai^) for the outgoing electrons are related to the operators aaai^) for the incident electrons via 
baa-is) = X]/3 Sa/3a/3(7(e), where Sa/3 denotes the scattering matrix. The scattering matrix is assumed to be spin- and 
energy-independent. Note that for a discrete energy spectrum in the leads, one can normalize the operators aaai^) 
such that {aQ,cr(£), ajj^, (e')} = 5a„'5ap5^^i /v, where the Kronecker symbol 5^^' equals 1 if e = e' and otherwise. 
Here, v stands for the density of states in the leads. Each lead is assumed to consist of only a single quantum channel; 
the generalization to leads with several channels is straightforward but not required here. The current Eq. (167) can 
be expressed in terms of the scattering matrix as 

^«(*) = 7^EE4.(^)^^7«7.(^')e'(^-^')*/", (168) 

ee' cr f^'f 

= SafjSay - S*^i3Saj. (169) 

The correlation function between the currents Ia{t) and If3{t) in two leads q,/3 = 1, ..,4 of the BS 

S^piuj) = lim — / e"^' Re TV [61^ {t)5Ip (0)x] , (170) 

where 5Ia = la — {^a), {la) = Tr(/Qx), v is the density of states in the leads, and x is the density matrix of the 
injected electron pair. Here, only the zero- frequency correlator Saj} = S^b{^) will be of interest (the dependence on 
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X was omitted). Further evaluation with x = |±)(±| yields 



^ T , (171) 

75 



where denotes the stun over 7 = 1,2 and all 5 7^ 7, and where again the upper (lower) sign refers to triplets 
(singlets). Note that the autocorrelator Saa is the shot noise in lead a. 

The result Eq. (171) can now be applied to the set-up shown in Fig. 31 involving four leads, described by the 
single-particle scattering matrix elements, S31 = S42 = r, and S41 = S32 = t, where r and t denote the reflection 
and transmission amplitudes at the beam splitter. In the absence of backscattering, S12 = S34 = Saa = 0, the noise 
correlations for the incident state |±) are 

e2 

533 = 544 = -534 = 2— T (1 - T) (1 T , (172) 

where T = |tp denotes the transmittivity of the beam splitter. For the remaining two triplet states | ft) and | H) one 
also obtains Eq. (172) with the upper sign. The mean current in lead a is, both for singlets and triplets, \{Ia)\ = e/hv. 
The noise-to-current ratio, or Fano factor, F = Saa/ \ {Ia)\ is thus found to be 

F = 2eT{l-T){lTSe,e,), (173) 

and correspondingly for the cross correlations. Eq. (173) implies that if electrons in the singlet state |— ) with equal 
energies, ei = £2, are injected pairwise into the leads 1 and 2, then the zero frequency noise is enhanced by a factor 
of two, F = 4eT(l — T), compared to the shot noise of imcorrelated particles (Biittiker, 1990, 1992; Khlus, 1987; 
Landaucr, 1989; Lesovik, 1989; Martin and Landauer, 1992), F = 2eT(l — T). This noise enhancement is due to 
bunching of electrons in the outgoing leads, caused by the symmetric orbital wavefunction of the spin singlet |— ). The 
triplet states |-|-), | tt)) and | J.|) exhibit antibunching, i.e. a complete suppression of the noise, Saa = 0. As already 
mentioned above, the noise enhancement for the singlet |— ) is a unique signature for entanglement (no unentangled 
state exists which can lead to this phenomenon) , therefore entanglement can be observed by measuring the noise power 
of a mesoscopic conductor as shown in Fig. 31. The various triplet states 1+), | tt)) and | H) cannot be distinguished 
by the noise measurement alone; this distinction requires a measurement of the spins of the outgoing electrons, e.g. 
by inserting spin-selective timneling devices (Prinz, 1998) into leads 3 and 4. The same signature of entanglement as 
for the shot noise can also be seen in the full counting statistics of the charge transport (Taddei and Fazio, 2002). 



4. Lower bounds for entanglement 

Here, the result of the previous Section is extended by providing a quantitative lower bound for the amount E 
of spin entanglement carried by individual pairs of electrons, related to the zero- frequency current correlators when 

measured in a beam splitter setup (Fig. 32). This result (Burkard and Loss, 2003) therefore relates experimentally 
accessible quantities with a measure for entanglement, the entanglement of formation E (Bennett et al., 1996b). Having 
information about E is important since it quantifies the usefulness of a bipartite state for quantum communication. 

Starting form a general state, expressed in the singlet-triplet basis as x = F\—){—\ +Go|+)(+| ^ Gi\ii){ii\ + 

Ax, where Ax are off-diagonal terms, one can decompose the current correlators Eq. (170) as 

Sa, ^ SI, = FS%-^ + GoSiy + Y: G,SIS, (174) 

i=Ta 

where S^^J = •S'jj*^^*'. The off-diagonal terms in Ax do not enter Sai3 because the operators dla{t) conserve total 
spin. The coefficients F, Go, G|, and Gi depend on the state preparation and therefore on the entangler. The 
more information about these coefficients can be gained, the better the chance to measure the entanglement of x- 
In Sec. V.B.3 it was shown that the singlet state |— ) gives rise to enhanced shot noise (and cross-correlators) at 
zero temperature, ^33^ = — «S'34^ = 2e/T(l — T)f, with the reduced correlator / = 2, as compared to the "classical" 
Poissonian value / = 1. All triplet states arc noiseless, s''^^ = sjJJ^ = s''^^ = (a, /? = 3, 4). Therefore, both the 
auto- and cross-correlations are only due to the singlet component of the incident two-particle state, 

^33 = -^34 = F5l*-> = 2e7T(l - T)f, f = 2F. (175) 

The entanglement of a bipartite pure state j?/') G Ha'Si'H.b can be measured by the von Neumann entropy S'n(IV')) = 
— TrsPB logpB (log in base 2) of the reduced density matrix pB = TrA|V')(^|) where < S'n < 1, SNd^^i)) = 
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FIG. 32 Entanglement of formation E of the electron 
spins versus singlet fidelity F and the reduced correlator 
/ = S33/2eIT{l - T). The curve illustrates the exact 
relation for Werner states. For general states, the curve 
is a lower bound for E\ allowed values for E and / (or 
F) are in the shaded region. 
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^Nd^i)) = 1, and S'NdV')) = ^ IV-) = ® \i')B- Physically, if S'n(IV')) ^ N/M then M > N copies of 

are sufficient to perform, e.g., quantum teleportation of N qubits for N,M ^ I (similarly for other quantum 
communication protocols). Generally, for a bipartite mixed state x the entanglement of formation (Bennett et at, 
1996b) is£'(x) = ^^^{i\x^),Pi)}e£{x)T,iPiS!^i\Xi)), where £(x) = {{\Xt),Pi)\T,tPt\Xi){Xi\ = x}, i-e-, the least expected 
entanglement of any ensemble of pure states realizing x- A state with i? > (£' = 1) is (maximally) entangled, and 
neither local operations nor classical communication (LOCC) between A and B can increase E. 

For an arbitrary mixed state of two qubits Xi -E'(x) is not a function of only the singlet fidelity F — (— |x|— ). 
However, E{x) — E{F) for the Werner states (Werner, 1989) 



the unique rotationally invariant states with singlet fidelity F, where \^±) = (|TT) i |||))/\/2. The entanglement of 
formation of the Werner states is known (Bennett et at, 1996b) as E{F) = E{pf) = ^2(1/2 + \/F{l - F)) if 1/2 < 
F < 1 and E{F) = E{pf) = if < F < 1/2, with the dyadic Shannon entropy iJ2 (a;) = -a; logx - (1 - x) log(l - x). 
With Eq. (175), one can now express E{pp) in terms of the reduced correlator / (Fig. 32). 

This result can be generalized to arbitrary mixed states x of two spins (qubits) using the following trick. Any state 
X can be transformed into pp with F = by a random bipartite rotation ("twirl") (Bennett et ai, 1996a,b), 

i.e. by applying U ®U with a random U G SU{2). Entanglement can only decrease (or remain constant) under the 
twirl operation because it involves only LOCC, 



Thus, the entanglement of formation E{F) of the corresponding Werner state represents a lower bound on E(x) 
(Fig. 32). A noise power exceeding / — 2F > 1 in the BS setup Fig. 31 can therefore be interpreted as a sign of 
entanglement, E{F) > 0, between the electron spins injected into leads 1 and 2. 

The lower bounds that have been discussed so far are only useful if one is assessing a source that aims at producing 
spin spinglet entanglement. However, it is possible in principle to extend this result to arbitrary entangled states if a 
means of rotating the spin of the carriers in one of the ingoing arms of the BS is available; such spin rotators could, 
e.g., be implemented by making use of the Rashba spin orbit effect (Egues et ai, 2002). 

The relation between the shot noise / and the entanglement E has also been explored in a number of non-ideal 
situations, in the presence of decoherence, backscattering, and thermally mixed input states (Burkard and Loss, 2003). 



5. Proposed tests of Bell's inequalities 

There have been a number of proposals to test Bell's inequalities (Bell, 1966; Mermin, 1993) with spm-entangled 
electrons directly with a spin-sensitive detection scheme (Kawabata, 2001; Lesovik et ai, 2001; Maitre et ai, 2000) in 
contrast to the detection scheme in Sees. V.B.3 and V.B.4 which only involves the measurement of a charge current. 
Combinations of the Andreev entangler setup with a Bell test were studied in (Chtchelkatchev et ai, 2002; Samuelsson 
et ai, 2004a). 

Bell tests for orbital entanglement with electron charge qubits in ballistic conductors have also been proposed 
(lonicioiu et al, 2001). A scheme to generate two-particle orbital entanglement in a mesoscopic normal-superconductor 




(176) 



E{F) < E{x). 



(177) 
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system and to detect it via a violation of a Bell inequality was analyzed in (Samuelsson et ai, 2003). A violation of a 
Bell inequality due to orbital entanglement in a Hanbury Brown-Twiss setting was also proposed (Samuelsson et al., 
2004b). 

C. Production of entangled photons 

The use of entangled photons from parametric down-conversion in non-linear crystals (Kwiat et al, 1995, 1999) 

has become a routine process (Boschi et ai, 1998; Bouwmccstcr et al, 1997; Mattlc et al, 1996). However, two 
disadvantages of these sources are (i) that they are quite inefficient (between 10^ and 10^° pump photons per yield 
photon) and (ii) that they are stochastic, i.e., although the rate (mean number of entangled pairs per second) is 
known, the precise instant of the emission cannot be controlled. It would be desirable to have a source of entangled 
photons that is both deterministic and efficient. The use of electron-hole recombination in a single QD was suggested 
in (Benson et ai, 2000; Moreau et ai, 2001). Non-resonant excitation of a QD is expected to produce pairs of 
entangled photons with an efficiency (production rate/pump rate) that is about four orders of magnitude bigger than 
for parametric down-conversion (Moreau et ai, 2001). Ultraviolet entangled photons have recently been generated in 
the semiconductor CuCl in a process called resonant hyper-parametric scattering (RHPS) which involves the creation 
of a virtual biexciton state in the crystal (Edamatsu et al., 2004). Although this is still a stochastic source, RHPS is 
a very efficient process. 

The production of polarization-entangled photons using the bicxcitonic groimd state has been investigated in two 
tunnel- coupled QDs (Gywat et ai, 2002). Biexcitons are bound states of two cxcitons in a semiconductor, where each 
exciton is the bound state of a negatively charged conduction-band electron and a positively charged valence-band 
hole. Excitonic absorption in single QDs has been studied theoretically (Efros and Efros, 1982), and biexcitonic 
states in single QDs have been investigated (Banyai et ai, 1988; Bryant, 1990; Hawrylak, 1999; Hu et ai, 1990; Kiraz 
et al, 2002; Nair and Takagahara, 1996; Santori et al, 2002; Takagahara, 1989). Single excitons in coupled QDs 
have been observed in experiment (Bayer et al, 2001; Schedelbeck et al, 1997) and spin spectroscopy of excitons in 
QDs was performed using polarization-resolved magnctophotoluminenscence (Johnston-Halperin et al, 2001). When 
discussing confined excitons, one needs to distinguish two regimes (Efros and Efros, 1982): (i) The weak confinement 
limit ax ae,ah, where ax is the radius of the free exciton and ag, ah the electron and hole effective Bohr radii in 
the QD. In this regime, an exciton can be considered (as in the bulk material) as a boson in an external confinement 
potential, (ii) The strong confinem,ent limit ax ^ a^-ah, where electrons and holes are separately confined in the 
QD. In this regime, electron-hole pairs cannot be considered as bosons anymore. In bulk GaAs ax ~ 10 nm, one is 
often in an intermediate regime ax ~ ae,ah for typical QD radii. 

Starting from a strong confinement ansatz and using the Heitler-London (HL) approximation, the low-energy (spin- 
resolved) biexciton spectrum (in which the electrons and holes each form cither a spin singlet or triplet) and the 
oscillator strengths Fig. 33, being a measure for the optical transition rates, have been calculated (Gywat et al, 
2002). The HL ansatz is similar to the one used for electrons in Sec. III.B.l with the Hamiltonian 

2 

= ^ ^/la^ + iTc + iTz + ifB, (178) 

a=e,h i=l 

where hai = {pai+qa^{''^ai) / c)^ /'^^a+Va^Tai) IS the single-particle Hamiltonian for the j-th electron (a = e,qe = — e) 
or hole {a = h, qh = +e). The Coulomb interaction is included by He = {^/'^)J2{a i)Ti(i3 j) qaqi3lK.\^m — r/s^l, with 
a dielectric constant k (for bulk GaAs, k = 13.18). A magnetic field B in 2; direction leads to orbital effects via the 
vector potential (in the symmetric gauge) A = B{—y, x, 0)/2 and to the Zeeman term Hz = i QoI^bBS"^, where 
ga is the effective g-factor of the electron (hole) and /ig is the Bohr magneton. 

Two-particle wave functions for electrons and for holes separately are constructed according to the HL method, i.e. 
a symmetric (|s)" = |/ = 0)", spin singlet) and an antisymmetric (|t)" = | J = 1)", spin triplet) linear combination of 
two-particle states \DD')c, = |D)a <S) \D')a, 

\ir = A^a/(|12)a + (-1)^|21)„), (179) 

where N^i = l/-\/2(l + (— 1)-'^|S'qP) and Sa = q(1|2)q denotes the overlap between the two orbital wave functions 
|l)a and |2)q. From the electron and hole singlet and triplet, the four biexciton states |/J) = (S> can be 
formed, where / = (1) for the electron singlet (triplet) and J = (1) for the hole singlet (triplet). The biexciton 
energies 

Eij = {IJ\H\IJ) =E° + E'' + EYj + Efj, (180) 
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FIG. 33 (a) and (b): Oscillator strengths fxx,x for transitions between the biexciton states — \IJ) and a single remaining 

exciton on one QD in units of fo as a function of (a) the magnetic field B (in Tesla) at i? = and (b) the electric field E (in 
mV/jUm) at i? = 0. The plotted HL energies Ejj are Ess (solid line), Est (short-dashed), Ets (dot-dashed), and Eu (dotted), 
neglecting the Zeeman energy. The parameters were chosen for GaAs with rj = uie/iUh = 1/2. (c) Entanglement of formation 
is as a function of the emission angles 8i and 02 of the two photons. Only if both photons are emitted perpendicular to the 
plane (^i = 62 = 0), the entanglement is maximal [E — 1). If at least one of them is emitted in-plane {6i — 7r/2), then the two 
photons are not entangled (_E = 0). 



with Efj EE {IJ\Ha\IJ), can be calculated analytically within HL (Gywat et ai, 2002). 

The exciton and biexciton recombination rates are determined by the oscillator strength / which is a measure for 
the coupling between the dipole moment of the exciton states and the electromagnetic field. For a transition between 
the -|- 1 and N exciton states jA^ + 1) and |A^), the oscillator strength is defined as 

fN^l^N = -^^1^, (181) 

where mo is the bare electron mass, hujN+i,N = En+i — En, and pNkx — {N + IjekA- p|-^), where GkA is the 
unit polarization vector for a photon with momentum k and helicity A = ±1, and p is the electron momentum 
operator. The inter-band momentum matrix element for a cubic crystal symmetry is given by M„x{9) = ekA'Pct)(o') — 
Pcii(cos(6') — ctA)/2 = Pcv'm<T\{0), where 9 is the angle between k and the normal to the plane of the 2D electron system 
(assuming that the latter coincides with one of the main axes of the cubic crystal), and Ep — 2p1^/m^ {— 25.7 eV for 
GaAs). 

The orbital momentum matrix element for the recombination of an exciton to the vacuum and for transitions from 
a biexciton state \XX) to a single exciton state \X) are given in (Gywat et ai, 2002); we give here our result for a 
transition between the HL biexciton states \XX) — \IJ) with one exciton on each QD and a single exciton in the 
final state \X), 

|(/J|ekA- P\X)\ = 2M„x{e)^JWj {Ceh [{^iy+-'+SeSh] + Seh [i-iySe + i-iySh]) . (182) 

The corresponding oscillator strength / is plotted in Figs. 33a and 33b, in the approximation hujxx,x ~ Eg. The 
effect of an electric field is to spatially separate the electrons from the holes, which leads to a reduction of the oscillator 
strengths. Hence, the optical transition rate can be efficiently switched off and on, thus allowing the deterministic 
emission of one photon pair. 

The HL biexciton state |/J) can be written as a superposition of dark {Sz — ±2) and bright {Sz — ±1) exciton 
states. Upon recombination of the biexciton, the emitted photon states are (up to normalization) 

|X/,7)« 1 + 1,01)1-1,02) +(-l)^+'Vl,^i)l + l,^2), (183) 

where \(J,d) = N{d){mc^+i{0)\a+) + ma-.-i{0)\a-)) is the state of a photon emitted from the recombination of an 
electron with spin Sz = cr/2 = ±1/2 and a heavy hole with spin Sz — 3(j/2 in a direction which encloses the angle 
9 with the normal to the plane of the 2D electron and hole motion. The states of right and left circular polarization 
are denoted |cr±). The entanglement between the two photon polarizations in the state Eq. (183) can be quantified 
by the von Neumann entropy E. For \ss) or \tt) and emission of the two photons enclosing an azimuthal angle 4> — 
or TT, we obtain 



(184) 
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where Xi = cos^(6'i). Wc plot E{9i,92) in Fig. 33c. Only the emission of both photons perpendicular to the plane 
(^1 = ^2 = 0) results in maximal entanglement {E = 1). The two photons are not entangled (E = 0) if at least one 
of them is emitted in-plane {6i = 7r/2). In order to observe the proposed effect, the relaxation rate to the biexciton 
ground state needs to be larger the biexciton recombination rate. The existence of such a regime is suggested by 
experiments with low excitation densities, see e.g. (Dekel et al., 2000; Ohnesorge et al., 1996). An upper limit for the 
pair production rate is then given by {tx + txx)~^, where tx is the exciton and txx the biexciton lifetime. 

References 

Aspect, A., J. Dalibard, and G. Roger, 1982, Phys. Rev. Lett. 49, 1804. 

Austing, D. G., T. Honda, K. Muraki, Y. Tokura, and S. Tarucha, 1998, Physica B 249-251, 206. 
Averin, D. V., 1998, Solid State Commun. 105, 659. 

Averin, D. V., and Y. V. Nazarov, 1992 (Plenum Press, New York), volume 294 of NATO ASI Series B: Physics. 
Awschalom, D. D., and J. M. Kikkawa, 1999, Phys. Today 52(6), 33. 

Awschalom, D. D., D. Loss, and N. Samarth (eds.), 2002, Semiconductor Spintronics and Quantum Computation, Nanoscience 

and Technology (Springer, Berlin). 
Bacon, D., K. R. Brown, and K. B. Whaley, 2001, Phys. Rev. Lett. 87, 247902. 
Bacon, D., J. Kempe, D. A. Lidar, and K. B. Whaley, 2000, Phys. Rev. Lett. 85, 1758. 
Ballentine, L. E., 1990, Quantum Mechanics (Prentice Hall, New Jersey), p. 352. 
Banyai, L., Y. Z. Hu, M. Lindbcrg, and S. W. Koch, 1988, Phys. Rev. B . 
Bao, ,J., A. V. Bragas, J. K. Furdyna, and R. Merlin, 2003, Nature Materials 2, 175. 

Barcnco, A., C. H. Bennett, R. Cleve, D. P. DiVincenzo, N. Margolus, P. Shor, T. Sleator, J. Smolin, and H. Weinfurter, 1995a, 

Phys. Rev. A 52, 3457. 
Barenco, A., D. Deutsch, A. Ekert, and R. Josza, 1995b, Phys. Rev. Lett. 74(20), 4083. 
Barnes, C. H. W., J. M. Shilton, and A. M. Robinson, 2000, Phys. Rev. B 62, 8410. 

Bayer, M., P. Hawrylak, K. Hinzer, S. Fafard, M. Korkusinski, Z. R. Wasilewski, O. Stern, and A. Forchel, 2001, Science 291, 
451. 

Beenakker, C. W. J., D. P. DiVincenzo, C. Emary and M. Kindcrmann, 2004, Phys. Rev. Lett. 93, 020501. 
Bell, J. S., 1966, Rev. Mod. Phys. 38, 447. 

Bena, C., S. Vishveshwara, L. Balents, and M. P. A. Fisher, 2002, Phys. Rev. Lett. 89, 037901. 

Bennett, C. H., G. Brassard, S. Popescu, B. Schumacher, J. A. Smolin, and W. K. Wootters, 1996a, Phys. Rev. Lett. 76, 722. 

Bennett, C. H., and D. P. DiVincenzo, 2000, Nature 404, 247. 

Bennett, C. H., D. P. DiVincenzo, ,J. A. Smohn, and W. K. Wootters, 1996b, Phys. Rev. A 54, 3824. 
Benson, O., C. Santori, M. Pclton, and Y. Yaniamoto, 2000, Phys. Rev. Lett. 84, 2513. 

Bertet, P., I. Chiorescu, G. Burkard, K. Semba, C. J. P. M. Harmans, D. P. DiVincenzo, and J. E. Mooij, 2004, Relcixation 

and dephasing in a persistent-current qubit, manuscript in preparation. 
Biolatti, E., R. C. lotti, P. Zanardi, and F. Rossi, 2000, Phys. Rev. Lett. 85, 5647. 

Bockrath, M., D. H. Cobden, J. Lu, A. G. Rinzler, R. E. Smalley, L. Balents, and P. L. McEuen, 1999, Nature 397, 598. 

Boncstccl, N. E., D. Stcpancnko, and D. P. DiVincenzo, 2001, Phys. Rev. Lett. 87, 207901. 

Boschi, D., S. Branca, F. De Martini, L. Hardy, and S. Popescu, 1998, Phys. Rev. Lett. 80, 1121. 

Bouwmeester, D., J.-W. Pan, K. Mattlc, M. Eibl, H. Weinfurter, and A. Zcihugcr, 1997, Nature 390, 575. 

Braunstein, S., and H.-K. Lo (eds.), 2000, Experimental Proposals for Quantum Computation, volume 48 (Wiley, Berlin), 

Special Focus Issue of Fortschritte der Physik. 
Brum, ,1. A., and P. Hawrylak, 1997, Supcrlattices and Microstructures 22, 431. 
Bryant, G. W., 1990, Phys. Rev. B 41, 1243. 
Bryant, G. W., 1993, Phys. Rev. B 48, 8024. 

Burkard, G., 2004, Circuit theory for decoherence in superconducting charge qubits, cond-mat/0408588. 

Burkard, G., and F. Brito, 2004, Non-additivity of decoherence rates in superconducting qubits, cond-mat/0408001. 

Burkard, G., D. P. DiVincenzo, P. Bertet, 1. Chiorescu, , and J. E. Mooij, 2004a, Asymmetry and decoherence in a double-layer 

persistent-current qubit, cond-mat/0405273. 
Burkard, G., H.-A. Engcl, and D. Loss, 2000a, Fortsch. Physik 48, 965, cond-mat/0004182. 
Burkard, G., R. H. Koch, and D. P. DiVincenzo, 2004b, Phys. Rev. B 69, 064503. 
Burkard, G., and D. Loss, 2002, Phys. Rev. Lett. 88, 047903. 
Burkard, G., and D. Loss, 2003, Phys. Rev. Letters 91, 087903. 
Burkard, G., D. Loss, and D. P. DiVincenzo, 1999a, Phys. Rev. B 59, 2070. 
Burkard, G., D. Loss, D. P. DiVincenzo, and J. A. Smolin, 1999b, Phys. Rev. B 60, 11404. 
Burkard, G., D. Loss, and E. V. Sukhorukov, 2000b, Phys. Rev. B 61, R16303, cond-mat/9906071. 
Burkard, G., G. Soolig, and D. Loss, 2000c, Phys. Rev. B 62, 2581. 
Biittiker, M., 1990, Phys. Rev. Lett. 65, 2901. 
Biittiker, M., 1992, Phys. Rev. B 46, 12485. 

Caldeira, A. O., and A. J. Leggett, 1983, Ann. Phys. (N.Y.) 143, 374. 

Chen, P., C. Piermarocchi, and L. J. Sham, 2001, Phys. Rev. Lett. 87, 067401. 



54 



Chiorescu, I., P. Bertet, K. Semba, Y. Nakamura, C. J. P. M. Harmans, and J. E. Mooij, 2004, Nature 431, 159. 
Chiorcscu, I., Y. Nakamura, C. J. P. M. Harmans, and J. E. Mooij, 2003, Science 299, 1869. 
Choi, M.-S., C. Bruder, and D. Loss, 2000, Phys. Rev. B 62, 13569. 

Chtchelkatchev, N. M., G. Blatter, G. B. Lesovik, and T. Martin, 2002, Phys. Rev. B 66, 161320. 

Ciorga, M., A. S. Sachrajda, P. Hawrylak, C. Gould, P. Zawadzki, S. Jullian, Y. Feng, and Z. Wasilewski, 2000, Phys. Rev. B 
61, R16315. 

Coish, W. A., and D. Loss, 2004, Hyperfine interaction in a quantum dot: Non-markovian electron spin dynamics, cond- 

mat/0405676 (to appear in Phys. Rev. B). 
Darwin, C. G., 1930, Math. Proc. Cambridge Phil. Soc. 27, 86. 

Dekel, E., D. Gershoni, E. Ehrenfreund, J. M. Garcia, and P. M. Petroff, 2000, Phys. Rev. B 61, 11009. 

Devoret, M. H., 1997, in Quantum fluctuations, lecture notes of the 1995 Les Houches summer school, edited by S. Reynaud, 

E. Giacobino, and J. Zinn-Justin (Elsevier, The Netherlands), p. 351. 
DiVincenzo, D. P., 1995, Phys. Rev. A 51(2), 1015. 

DiVincenzo, D. P., 1997, Mesoscopic Electron Transport (Kluwer Academic Publishers, Dordrecht), volume 345 of NATO ASI 

Series E, p. 657, cond-mat/9612126. 
DiVincenzo, D. P., 2000, Fortschr. Phys. 48, 771, quant-ph/0002077. 

DiVincenzo, D. P., D. Bacon, J. Kempe, G. Burkard, and K. B. Whaley, 2000, Nature 408, 339. 

DiVincenzo, D. P., G. Burkard, D. Loss, and E. V. Sukhorukov, 1999, in Quantum Mesoscopic Phenomena and Mesoscopic 
Devices in Microelectronics, edited by I. O. Kulik and R. EUialtioglu (Kluwer, Dordrecht), volume 559 of NATO ASI Series 

C: Mathematical and Physical Sciences, pp. 399-428. 
DiVincenzo, D. P., and D. Loss, 1999, J. Mag. Magn. Matl. 200, 202. 

Dixon, D. C, K. R. Wald, P. L. McEuen, and M. R. Melloch, 1997, Phys. Rev. B 56, 4743. 

Dobers, M., K. v. Klitzing, J. Schneider, G. Weimann, and K. Ploog, 1988, Phys. Rev. Lett. 61, 1650. 

Edamatsu, K., G. Oohata, R. Shimizu, and T. Itoh, 2004, Nature 431, 167. 

Efros, A. L., and A. L. Efros, 1982, Sov. Phys. Semicond. 16, 772. 

Egger, R., and A. Gogohn, 1997, Phys. Rev. Lett. 79, 5082. 

Egues, J. C, G. Bmkard, and D. Loss, 2002, Phys. Rev. Lett. 89, 176401. 

Einstein, A., B. Podolsky, and N. Rosen, 1935, Phys. Rev. 47, 777. 

Elzerman, J. M., R. Hanson, J. S. Greidanus, L. H. Willems van Beveren, S. De Pranceschi, L. M. K. Vandersypen, S. Tarucha, 

and L. P. Kouwenhoven, 2003, Phys. Rev. B 67, 161308. 
Engel, H.-A., L. P. Kouwenhoven, D. Loss, and C. M. Marcus, 2004, Controlling spin qubits in quantum dots, cond-mat/0409294. 
Engcl, H.-A., and D. Loss, 2001, Phys. Rev. Lett. 86, 4648. 

Fafard, S., Z. R. Wasilewski, C. N. Allen, D. Picard, M. Spanner, J. P. McCaffrey, and P. G. Piva, 1999, Phys. Rev. B 59, 
15368. 

Fazio, R., G. M. Palma, and J. Siewert, 1999, Phys. Rev. Lett. 83, 5385. 

Feynman, R. P., R. B. Leighton, and M. Sands, 1965, The Feynman Lectures, volume 3 (Addison- Wesley, Reading). 
Fock, v., 1928, Z. Phys. 47, 446. 

Fradkiri, E., 1991, Field Theories of Condensed Matter Systems (Addison- Wesley, Reading, MA). 
Fricke, M., A. Lorke, ,L P. Kotthaus, G. Medeiros-Ribciro, and M. P. Petroff, 1996, Europhys. Lett. 36, 197. 
Gantmakher, V. F., and Y. B. Levinson, 1987, Carrier scattering in metals and semiconductors (North-Holland, Amsterdam), 
chapter 1. 

Gerard, J.-M., and B. Gayral, 1999, in QED phenomena and applications of microcavities and photonic crystals (Springer- 
Verlag, Berlin). 

Giuliani, G. F., and ,J. ,]. Quinn, 1982, Phys. Rev. B 26, 4421. 

Golovach, V. N., A. Kliaetskii, and D. Loss, 2004, Phys. Rev. Lett. 93, 016601. 

Gupta, J. A., D. D. Awsclialom, X. Peng, and A. P. Alivisatos, 1999, Phys. Rev. B 59, R10421. 

Gywat, O., G. Burkard, and D. Loss, 2002, Phys. Rev. B 65, 205329. 

Hanbury Brown, R., and R. Q. Twiss, 1956, Nature 177, 27. 

Hanson, R., J. M. Elzerman, L. H. W. van Beveren, L. M. K. Vandersypen, and L. P. Kouwenhoven, 2004, in Controlling 

Decoherence, 330. W&Heraeus Seminar, Bad Honnef, Germany. 
Hanson, R., B. Witkamp, L. M. K. Vandersypen, L. H. W. van Beveren, J. M. Elzerman, and L. P. Kouwenhoven, 2003, Phys. 

Rev. Lett. 91, 196802. 
Hawrylak, P., 1999, Phys. Rev. B 60, 5597. 

Hekking, F. W. J., L. L Glazman, K. A. Matveev, and R. L Shekhter, 1993, Phys. Rev. Lett. 70, 4138. 

Henny, M., S. Oberholzer, C. Strunk, T. Heinzel, K. EnssUn, M. Holland, and C. Schonenberger, 1999, Science 284, 296. 

Hu, X., and S. Das Sarma, 2000, Phys. Rev. A 61, 062301. 

Hu, X., and S. Das Sarma, 2001, Phys. Rev. A 64, 042312. 

Hu, X., and S. Das Sarma, 2003, Phys. Rev. A 68, 052310. 

Hu, X., and S. Das Sarma, 2004, Phys. Rev. B 69, 115312. 

Hu, Y. Z., S. W. Koch, M. Lindberg, N. Peyghambarian, E. L. Pollock, , and F. F. Abraham, 1990, Phys. Rev. Lett. 64, 1805. 
Imamoglu, A., D. D. Awschalom, G. Burkard, D. P. DiVincenzo, D. Loss, M. Sherwin, and A. Small, 1999, Phys. Rev. Lett. 
83, 4204. 

Imamura, H., H. Aoki, and P. A. Maksym, 1998, Phys. Rev. B 57, R4257. 
Imamura, H., P. A. Maksym, and H. Aoki, 1996, Phys. Rev. B 53, 12613. 



55 



Imamura, H., P. A. Maksym, and H. Aoki, 1999, Phys. Rev. B 59, 5817. 

loffe, L. B., V. B. Gcshkcnbcin, M. V. Fcigcl'man, A. L. Fauchcrc, and G. Blatter, 1999, Nature 398, 679. 

lonicioiu, R., P. Zanardi, and F. Rossi, 2001, Phys. Rev. A 63, 050101. 

Ivchenko, E. L., A. A. Kiselev, and M. Willander, 1997, Solid State Comm. 102, 375. 

Johnston-Halperin, E., D. D. Awschalom, S. A. Crooker, A. L. Efros, M. Rosen, X. Peng, and A. P. Alivisatos, 2001, Phys. 

Rev. B 63, 205309. 
Kane, B. E., 1998, Nature 393(6681), 133. 

Kane, C, L. Balonts, and M. Fisher, 1997, Phys. Rev. Lett. 79, 5086. 
Kavokin, K. V., 2001, Phys. Rev. B 64, 075305. 
Kawabata, S., 2001, ,J. Phys. Soc. Jpn. 70, 1210. 

Kempe, J., D. Bacon, D. A. Lidar, and K. B. Whaley, 2001, Phys. Rev. A 63, 042307. 
Kempe, J., and K. B. Whaley, 2002, Phys. Rev. A 65, 052330. 
Khaetskii, A. V., 2001, Physica E 10, 27. 

Khaetskii, A. V., D. Loss, and L. Glazman, 2002, Phys. Rev. Lett. 88, 186802. 
Khaetskii, A. V., and Y. V. Nazarov, 2000, Phys. Rev. B 61, 12639. 
Khaetskii, A. V., and Y. V. Nazarov, 2001, Phys. Rev. B 64, 125316. 
Khaneja, N., R. Brockett, and S. J. Glaser, 2001, Phys. Rev. A 63, 032308. 
Khlus, V. A., 1987, Zh. Eksp. Teor. Fiz. 93, 2179. 

Kikkawa, J. M., and D. D. Awschalom, 1998, Phys. Rev. Lett. 80, 4313. 

Kikkawa, J. M., L P. Smorchkova, N. Samarth, and D. D. Awschalom, 1997, Science 277, 1284. 

Kiraz, A., S. Faith, C. Bechcr, B. Gayral, W. V. Schoenfeld, P. M. Potroff, L. Zhang, E. Hu, and A. Imamoglu, 2002, Phys. 
Rev. B 65, 161303. 

Kiselev, A. A., E. L. Ivchenko, and U. Rossler, 1998, Phys. Rev. B 58, 16353. 
Kitaev, A. Y., 2003, Annals Phys. 303, 2, quant-ph/9707021. 

Koch, R., J. Kirtley, J. Rozen, J. Sun, G. Keefe, F. Milliken, C. Tsuei, and D. DiVincenzo, 2003, Bull. Am. Phys. Soc. 48, 367. 
Koch, R., J. Kirtley, J. Rozen, J. Sun, G. Keefe, F. Milliken, C. Tsuei, and D. DiVincenzo, 2004, unpublished. 
Konig, J., H. Schoeller, and G. Schon, 1997, Phys. Rev. Lett. 78, 4482. 
Kouwonhovcn, L. P., D. G. Austing, and S. Tarucha, 2001, Rep. Prog. Phys. 64, 701. 

Kouwcnhovcn, L. P., C. M. Marcus, P. L. McEuen, S. Tarucha, R. M. Westervelt, and N. S. Wingreen, 1997a (Kluwer, 

Dordrecht), volume 345 of NATO ASI Series E: Applied Sciences, pp. 105-214. 
Kouwenhoven, L. P., T. H. Oosterkamp, M. W. S. Danoesastro, M. Eto, D. G. Austing, and S. Taurcha, 1997b, Science 278, 

1788. 

Kouwonhovcn, L. P., G. Schon, and L. L. Sohn (ods.), 1997c, Mesoscopic Electron Transport, volume 345 of NATO ASI Series 

E (Kluwer Academic Publishers, Dordrecht). 
Kwiat, P. G., K. Mattle, H. Weinfurter, A. Zeilinger, A. V. Sergienko, and Y. Shih, 1995, Phys. Rev. Lett. 75, 4337. 
Kwiat, P. G., E. Waks, A. G. White, I. Appelbaum, and P. H. Eberhard, 1999, Phys. Rev. A 60, R773. 

Kyriakidis, J., and S. J. Penney, 2004, Coherent rotations of a single spin-based qubit in a single quantum dot at fixed zeeman 
energy, cond-mat/0407139. 

Kyriakidis, J., M. Pioro-Ladricrc, M. Ciorga, A. S. Sachrajda, and P. Hawrylak, 2002, Phys. Rev. B 66, 035320. 
Ladd, T. D., J. R. Goldman, F. Yamaguchi, Y. Yamamoto, E. Abe, and K. M. Itoh, 2002, Phys. Rev. Lett. 89, 017901. 
Landauer, R., 1989, Physica D 38, 226. 
Landauer, R., 1991, Physics Today 44(5), 23. 
Landauer, R., 1996, Science 272, 1914. 

Loggctt, A. ,J., S. Chakravarty, A. T. Dorscy, M. P. A. Fisher, A. Garg, and W. Zwerger, 1987, Rev. Mod. Phys. 59, 1. 

Losovik, G., T. Martin, and G. Blatter, 2001, Eur. Phys. J. B 24, 287. 

Lesovik, G. B., 1989, Pis'ma Zh. Eksp. Teor. Fiz. 49, 513. 

Leung, D. W., 2004, Int. J. Quant. Inf. 2, 33. 

Levy, J., 2001, Phys. Rev. A 64, 052306. 

Levy, J., 2002, Phys. Rev. Lett. 89, 147902. 

Lidar, D. A., D. Bacon, and K. B. Whaley, 1999, Phys. Rev. Lett. 82, 4556. 

Lidar, D. A., I. L. Chuang, and K. B. Whaley, 1998, Phys. Rev. Lett. 81, 2594. 
Liu, R. C, B. Odom, Y. Yamamoto, and S. Tarucha, 1998, Nature 391, 263. 

Livermore, C, C. H. Crouch, R. M. Westervelt, K. L. Campman, and A. C. Gossard, 1996, Science 274, 1332. 

Loss, D., and D. P. DiVincenzo, 1998, Phys. Rev. A 57(1), 120. 

Loss, D., and P. Goldbart, 1992, Phys. Rev. B 45, 13544. 

Loss, D., and E. V. Sukhorukov, 2000, Phys. Rev. Lett. 84, 1035. 

Loudon, R., 1998, Phys. Rev. A 58, 4904. 

Luyken, R. J., A. Lorke, M. Haslinger, B. T. Miller, M. Fricke, J. P. Kotthaus, G. Medeiros-Ribiero, and P. M. Petroff, 1998, 
Physica E 2, 704. 

Madelung, O., 1978, Introduction to Solid-State Theory (Springer, Berlin), p. 110. 
Mahan, G. D., 1993, Many Particle Physics (Plenum, New York), 2nd edition. 
Maitre, X., W. D. Oliver, and Y. Yamamoto, 2000, Physica E 6, 301. 
Makhhn, Y., 2002, Quant. Info. Proc. 1, 243, quant-ph/0002045. 



56 



Makhlin, Y., G. Schon, and A. Shnirman, 1999, Nature 398, 305. 
Makhlin, Y., G. Schon, and A. Shnirman, 2001, Rev. Mod. Phys 73, 357. 
Makhhn, Y., and A. Shnirman, 2004, Phys. Rev. Lett. 92, 178301. 
Martin, T., and R. Landauer, 1992, Phys. Rev. B 45, 1742. 

Martinis, J. M., S. Nam, J. Aumentado, , K. M. Lang, and C. Urbina, 2003, Phys. Rev. B 67, 094510. 

Martinis, J. M., S. Nam, J. Aumentado, and C. Urbina, 2002, Phys. Rev. Lett. 89, 117901. 

Mattis, D. C., 1988, The Theory of Magnetism I (Springer, New York), volume 17 of Springer Series in Solid-State Sciences, 
chapter 4.5. 

Mattle, K., H. Wcinfurtcr, P. G. Kwiat, and A. Zcilinger, 1996, Phys. Rev. Lett. 76, 4656. 
Meier, F., J. Levy, and D. Loss, 2003a, Phys. Rev. B 68, 134417. 
Meier, F., J. Levy, and D. Loss, 2003b, Phys. Rev. Lett. 90, 047901. 

Meier, F., and D. Loss, 2004, Reduced visibility of rabi oscillations in superconducting qubits, cond-mat/0408594. 

Mermin, N. D., 1993, Rev. Mod. Phys. 65, 803. 

Mermin, N. D., 2003, Am. J. Phys. 71, 23, quant-ph/0207118. 

Mooij, J. E., T. P. Orlando, L. Levitov, L. Tian, C. H. van der Wal, and S. Lloyd, 1999, Science 285, 1036. 
Moreau, E., I. Robert, L. Manin, V. Thierry-Mieg, J. M. Gerard, and I. Abram, 2001, Phys. Rev. Lett. 87, 183601. 
Nair, S. V., and T. Takagahara, 1996, Phys. Rev. B 55, 5153. 
Nakamura, Y., Y. A. Pashkin, and J. S. Tsai, 1999, Nature 398, 786. 
Nielsen, M. A., 2003, Phys. Lett. A 308, 96. 

Nielsen, M. A., and I. L. Chuang, 2000, Quantum Computation and Quantum Information (Cambridge University Press, 
Cambridge) . 

Oh, J. H., K. J. Chang, G. Ihm, and S. J. Lee, 1996, Phys. Rev. B 53, R13264. 

Ohnesorge, B., M. Albrecht, J. Oshinowo, A. Forchel, and Y. Arakawa, 1996, Phys. Rev. B 54, 11532. 
Oliver, W. D., J. Kim, R. C. Liu, and Y. Yamamoto, 1999, Science 284, 299. 
Oliver, W. D., F. Yamaguchi, and Y. Yamamoto, 2002, Phys. Rev. Lett. 88, 037901. 

Orlando, T. P., J. E. Mooij, L. Tian, C. H. van der Wal, L. Levitov, S. Lloyd, and J. J. Mazo, 1999, Phys. Rev. B 60, 15398. 

Pashkin, A., T. Yamarnoto, O. Astaficv, Y. Nakamura, D. V. Averin, and J. S. Tsai, 2003, Nature 421, 823. 
Peikari, B., 1974, Fundamentals of Network Analysis and Synthesis (Prcnticc-Hall, Englcwood Cliffs, NJ). 
Pellizzari, T., S. A. Gardiner, J. I. Cirac, and P. ZoUer, 1995, Phys. Rev. Lett. 75, 3788. 
Piermarocchi, C, P. Chen, L. J. Sham, and D. G. Steel, 2002, Phys. Rev. Lett. 89, 167402. 
Platzman, P. M., and M. I. Dykman, 1999, Science 284, 1967. 
Popescu, A. E., and R. lonicioiu, 2004, Phys. Rev. B 69, 245422. 

Proskill, J., 1998a, Fault-tolerant quantum computation (World Scientific, Singapore), pp. 213-269. 
Preskill, J., 1998b, Proc. R. Soc. London Ser. A 454, 385. 
Prinz, G. A., 1998, Science 282, 1660. 

Privman, V., I. D. Vagner, and G. Kventsel, 1998, Phys. Lett. A 239, 141. 
Raussendorf, R., and H. J. Briegel, 2001, Phys. Rev. Lett. 86, 5188. 

Recher, P., and D. Loss, 2002a, Journal of Superconductivity and Novel Magnetism 15, 49, cond-mat/0202102. 

Recher, P., and D. Loss, 2002b, Phys. Rev. B 65, 165327. 

Recher, P., and D. Loss, 2003, Phys. Rev. Lett. 91, 267003. 

Recher, P., E. V. Sukhorukov, and D. Loss, 2001, Phys. Rev. B 63, 165314. 

Redfield, A. G., 1957, IBM J. Res. Develop. 1, 19. 

Requist, R., J. Schliemann, A. G. Abanov, and D. Loss, 2004, Double occupancy errors in quantum computing operations: 

Corrections to adiabaticity, cond-rnat/0409096. 
Sanmelssou, P., E. V. Sukhorukov, and M. Buttiker, 2003, Phys. Rev. Lett. 91, 157002. 
Sarrmolsson, P., E. V. Sukhorukov, and M. Buttiker, 2004a, Phys. Rev. B 70, 115330. 
Sarrmclsson, P., E. V. Sukhorukov, and M. Buttiker, 2004b, Phys. Rev. Lett. 92, 026805. 
Sanders, G. D., K. W. Kim, and W. C. Holton, 1999, Phys. Rev. A 59, 1098. 
Santori, C, G. S. Solomon, M. Pelton, and Y. Yamamoto, 2002, Phys. Rev. B 65, 073310. 
Saraga, D. S., B. L. Altshuler, D. Loss, and R. M. Westervelt, 2004, Phys. Rev. Lett. 92, 246803. 
Saraga, D. S., and D. Loss, 2003, Phys. Rev. Lett. 90, 166803. 

Schcdelbeck, G., W. Wegscheider, M. Bichler, and G. Abstreiter, 1997, Science 278, 1792. 
Schliemann, ,L, ,1. 1. Cirac, M. Kus, M. Lewenstein, and D. Loss, 2001a, Phys. Rev. A 64, 022303. 
Schliemann, J., D. Loss, and A. H. MacDonald, 2001b, Phys. Rev. B 63, 085311. 
SchriefFer, J. R., 1964, Theory of Superconductivity (Benjamin/Cummings, New York). 
Schrieffer, J. R., and P. A. Wolff, 1966, Phys. Rev. 149, 491. 
Schuch, N., and J. Siewert, 2003, Phys. Rev. A 67, 032301. 
Schulz, H. J., 1990, Phys. Rev. Lett. 64, 2831. 

Seelig, G., 1999, Spin interactions in tunneling- coupled vertical quantum dots, Diploma thesis. University of Basel. 
Shnirman, A., G. Schon, and Z. Herrnon, 1997, Phys. Rev. Lett. 79, 2371. 

Simmonds, R. W., K. M. Lang, D. A. Kite, S. Nam, D. P. Pappas, and J. M. Martinis, 2004, Phys. Rev. Lett. 93, 077003. 
de Sousa, R., X. Hu, and S. Das Sarma, 2001, Phys. Rev. A 64, 042307. 

Taddei, F., and R. Fazio, 2002, Phys. Rev. B 65, 075317. 
Takagahara, T., 1989, Phys. Rev. B 39, 10206. 



57 



Tanamoto, T., 2000, Phys. Rev. A 61, 022305. 

Tarucha, S., D. G. Austing, T. Honda, R. J. van dor Hago, and L. P. Kouwcnhovcn, 1996, Phys. Rev. Lett. 77(17), 3613. 
Tian, L., L. S. Lovitov, J. E. Mooij, T. P. Orlando, C. H. van dcr Wal, and S. Lloyd, 2000, Quantum Mesoscopic Phenomena 

and Mesoscopic Devices in Microelectronics (Kluwer, Dordrecht), pp. 429-438, cond-mat/9910062. 
Tian, L., S. Lloyd, and T. P. Orlando, 2002, Phys. Rev. B 65, 144516. 
Tokura, Y., D. G. Austing, and S. Tarucha, 1999, J. Phys. C 11, 6023. 

Topinka, M. A., B. J. LeRoy, S. E. J. Shaw, E. J. Heller, R. M. Westervelt, K. D. Maranowski, and A. C. Gossard, 2000, 
Science 289, 2323. 

Topinka, M. A., B. J. Leroy, R. M. Westervelt, S. E. J. Shaw, R. Fleischmann, E. J. Heller, K. D. Maranowski, and A. C. 

Gossard, 2001, Nature 410, 183. 
Torres, J., and T. Martin, 1999, Eur. Phys. J. B 12, 319. 
Troiani, F., U. Hohenester, and E. Mohnari, 2000, Phys. Rev. B 62, R2263. 

Tsvelik, A. M., 2003, Quantum Field Theory in Condensed Matter Physics (Cambridge University Press), 2nd edition. 
Vala, J., and K. B. Whaley, 2002, Phys. Rev. A 66, 022304. 

Viola, L., E. Knill, and S. Lloyd, 2000, Dynamical generation of noiseless quantum subsystems, quant-ph/0002072. 
Vion, D., A. Aassime, A. Cottet, P. Joyez, H. Pothier, C. Urbina, D. Esteve, and M. H. Devoret, 2002, Science 296, 886. 
Vorojtsov, S., E. R. Mucciolo, and H. U. Baranger, 2004, Phys. Rev. B 69, 115329. 

Vrijen, R., E. Yablonovitch, K. Wang, H. W. Jiang, A. Balandin, V. Roychowdhury, T. Mor, and D. DiVincenzo, 2000, Phys. 
Rev. A 62, 012306. 

Wagner, M., U. Merkt, and A. V. Chaplik, 1992, Phys. Rev. B 45, 1951. 

van der Wal, C. H., A. C. J. ter Har, F. K. Wilhelm, R. N. Schoutcn, C. J. P. M. Harmans, T. P. Orlando, S. Lloyd, and J. E. 

Mooij, 2000, Science 290, 773. 
van der Wal, C. H., F. K. Wilhelm, C. J. P. M. Harmans, and J. E. Mooij, 2003, Eur. Phys. J. B 31, 111. 
Wallraff, A., D. L Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, 2004, 

Nature 431, 162. 

Waugh, F. R., M. J. Berry, D. J. Mar, R. M. Westervelt, K. L. Campman, and A. C. Gossard, 1995, Phys. Rev. Lett. 75(4), 
705. 

Weinstein, Y. S., and C. S. Hellberg, 2004, Energetic suppression of decoherence in exchange-only quantum computation, 
quant-ph/0408037. 

Weiss, U., 1999, Quantum Dissipative Systems (World Scientific, Singapore), 2nd edition. 
Werner, M. J., and P. D. Drummond, 1991, Phys. Rev. A 43, 6414. 

Werner, R. F., 1989, Phys. Rev. A 40, 4277. 

van der Wiel, W. G., S. D. Franceschi, J. M. Elzerman, T. Fujisawa, S. Tarucha, and L. P. Kouwenhoven, 2003, Rev. Mod. 
Phys. 75, 1. 

Wilhelm, F. K., M. J. Storcz, C. H. van der Wal, C. J. P. M. Harmans, and J. E. Mooij, 2003, Adv. Solid State Phys. 43, 763. 

Wojs, A., P. Hawrylak, S. Fafard, and L. Jacak, 1996, Phys. Rev. B 54, 5604. 

Wu, L.-A., and D. A. Lidar, 2002a, Phys. Rev. A 65, 042318. 

Wu, L.-A., and D. A. Lidar, 2002b, Phys. Rev. A 66, 062314. 

Wu, L.-A., D. A. Lidar, and M. Friesen, 2004, Phys. Rev. Lett. 93, 030501. 

Yang, S.-R. E., J. Schliemann, and A. H. MacDonald, 2002, Phys. Rev. B 66, 153302. 

You, J. Q., C.-H. Lam, and H. Z. Zheng, 2001, Phys. Rev. B 63, 180501(R). 

Yurke, B., 1987, J. Opt. Soc. Am. B 4, 1557. 

Yurke, B., and J. S. Denker, 1994, Phys. Rev. A 29, 1419. 

Zanardi, P., and M. Rasetti, 1997, Mod. Phys. Lett. B 11, 1085. 

Zanardi, P., and F. Rossi, 1998, Phys. Rev. Lett. 81, 4752. 

Zumbiihl, D. M., C. M. Marcus, M. P. Hanson, and A. C. Gossard, 2004, Cotunneling spectroscopy in few-electron quantum 
dots, cond-mat/0408276. 



